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STATEMENT  OF  THE  PROBLEM 

The  problem  undertaken  involved  the  determina¬ 
tion  of  the  pipeline  flow  characteristics  of  the  following 
crude  oils: 

1.  The  Saskatchewan  blend,  consisting  of 
39.0%  Fosterton  crude,  34.5%  Success 
crude,  and  26,5%  Gantaur. 

2.  The  Combined  blend.,  consisting  of  80% 
Saskatchewan  blend,  plus  20%  Redwater 
crude . 

3.  The  Fosterton  crude. 

4.  The  Pembina  crude. 

The  MacMichael  and  Alberta  concentric  cylinder 
rotational  viscometers,  and  the  Pilot  Pipeline  viscometer 
were  used  to  determine  the  relationship  of  the  apparent 
viscosity  to  the  temperature  and  rate  of  shear.  The 
comparison  and  correlation  of  these  results  was  the  main 
objective . 
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ABSTRACT 

Several  Western  Canadian  crudes  were  tested  in  a 
pilot  pipeline  and  in  two  concentric  cylinder  viscometers 
to  determine  their  apparent  viscosity-shear  rate  relation¬ 
ships.  Data  obtained  on  the  five  different  liquids  indicate 
characteristics  which  classify  them  into  three  groups. 

The  first  group  consists  of  the  Saskatchewan  blend, 

the  Combined  blend,  and  the  Fosterton  crude.  Below  40<JF  to 
o 

45  F,  the  pilot  pipeline  data  indicated  Newtonian  behavior, 

whereas  the  concentric  cylinder  data  indicated  pseudoplastic 

o  o 

behavior.  Through  the  temperature  range  of  75  F  to  80  F,the  a 
apparent  viscosity  data  obtained  by  the  pilot  pipeline  was 
25t  to  40/o  higher  than  the  concentric  cylinder  data.  At  40°F 
to  45°F,  the  data  from  both  methods  correlate,  and  at  20°F 
the  pilot  pipeline  data  are  1+0%  lower  than  the  concentric 
cylinder  data. 

Apparent  viscosity  data  were  obtained  on  three 
different  samples  of  Redwater  crude  oil.  The  data ' obtained 
for  two  samples  by  A.  Masuda  (19J,  on  the  pilot  pipeline, 
are  generally  20%  lower  than  the  concentric  cylinder  data 


. 


' 


IV 


o  o  o 

on  the  other  sample  through  40  F  to  80  F.  Below  40  F,  the 

concentric  cylinder  data  indicate  pseuodplastic  effects 

o 

with  the  result  that  the  data  at  20  F  are  up  to  75/°  higher 
than  the  pilot  pipeline  data.  With  the  Redwater  crude,  no 
pseudoplastic  effects  were  determined  on  the  pilot  pipeline 
throughout  the  temperature  range  investigated. 

The  pseudoplastic  behavior  of  the  Pembina  crude  was 
detected  by  both  methods  of  investigation.  Above  40°F,  the 
apparent  viscosity-shear  rate  curves  obtained  by  the  Alberta 
viscometer  have  a  similar  curvature,  and  correlate  approximately 
to  the  model  pipeline  curve  at  100  sec.  \  The  pilot  pipeline 
indicates  that  the  shear  rate  above  which  the  apparent  viscosity 
is  constant  for  the  Pembina  crude  is  approximately  10,000sec. 

-1 

A  large  family  of  curves  with  shear  rate  varying  from  100  sec. 

to  80,000  sec.  ''“was  determined  by  the  pilot  pipe  line. 

The  poor  correlation  on  all  the  crudes  may  be 

attributed  to  several  characteristics.  The  concentric 

cylinder  instrument  cannot  reproduce  the  high  shear  rates 

present  in  the  pilot  pipe  line.  With  liquid  flowing  in  the 

laminar  region  in  the  pilot  pipeline,  a  temperature 

gradient  increases  the  friction  factor  as  much  as  k5%, 

which  results  in  an  increased  value  for  the  apparent 

viscosity.  The  rapid  loss  of  light  ends  from  the  crude 

o 

oil  samples  above  60  F  is  serious  on  the  concentric 
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cylinder  viscometer.  The  modification  of  the  equipment  may  improve 
the  correlation  so  that  the  rotational  viscometer  can  be  used  in 
place  of  the  pilot  pipeline. 

On  the  pilot  pipeline,  the  increase  in  the  apparent 
viscosity  with  pressure  was  found  to  be  7,  4,  and  14  $/100psi. ,  for 
the  Combined  blend,  the  Fosterton  crude,  and  the  Pembina  crude 
respectively.  For  the  Combined  blend  and  the  Fosterton  crude,  the 
effect  of  pressure  decreased  with  increasing  temperature. 
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INTRODUCTION 


When  Western  Canada's  crude  oil  producing  capacity 
became  larger  than  the  local  refining  market  could  absorb, 
it  became  necessary  to  move  large  volumes  of  crude  oil  to 
distant  markets.  The  most  efficient  and  economical  method 
of  moving  these  large  volumes  overland  for  either  long  or 
short  distances,,  is  by  pipeline.  This  surplus  of  crude  oil 
has  resulted  in  the  construction  of  several  thousands  of 
miles  of  pipelines  of  all  diameter  to  gather  and  deliver 
crude  oils  varying  in  gravity  from  below  20  °API  to  above 
40  °API.. 

In  the  construction  of  this  delivery  system,  it 
has  been  necessary  to  solve  many  new  problems  connected 
with  the  gathering,  storage,  and  transportation  of  crude  oil. 
Besides  the  problem  of  wilderness  terrain  and  muskeg  in 
Western  Canada,  there  is  the  problem  of  designing  a  system 
which  will  operate  economically  under  a  large  variation  in 
temperature.  At  the  lower  temperatures,  the  value  of  the 
apparent  viscosity  increases  very  quickly  with  a  change  in 
the  temperature,  and  for  some  crude  oils  the  operating 
temperature  is  below  the  wax  point.  The  higher  values  of 
the  apparent  viscosity  not  only  increase  the  pressure  drop 
in  the  pipe  lines  and  manifold  systems,  but  the  hydraulic 
efficiency  of  the  pumping  equipment  is  reduced.  Because  the 
wax,  in  the  form  of  a  sludge,  tends  to  settle  out  at  certain 
temperatures  in  the  pipeline,  it  is  necessary  to  run  plugs 
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more  often  at  these  temperatures  to  obtain  better  pipeline 
efficiency.  It  has  been  determined  that  it  is  economical 
to  heat  the  oil  at  the  pumping  stations  on  long  pipelines 
to  decrease  the  apparent  viscosity. 

When  it  is  necessary  to  pump  the  crude  oil  in 
batches,  the  temperature  of  the  oil  in  the  stock  tanks  drops, 
so  that  large  suction  lines  to  the  pumps  are  required  to 
enable  the  centrifugal  pumps  to  operate  efficiently.  To 
prevent  damage  to  the  impellers  of  the  centrifugal  pumps,  it 
is  necessary  to  install  suction  lines  from  the  storage  tanks, 
of  sufficient  size  to  prevent  the  occurrence  of  a  vapor 
pocket.  With  new  pipelines  being  built  to  handle  crudes 
with  pseudoplastic  properties,  many  new  problems  will  be 
encountered,  especially  on  the  gathering  systems. 

The  accurate  design  of  the  pipelines  and  pumping 
equipment  requires  the  accurate  determination  of  the  apparent 
viscosity-temperature  relationship.  The  determination  of  the 
apparent  viscosity-temperature  values  at  two  or  three  tempera¬ 
tures  by  means  of  precision  capillary  viscometers  has  been 
the  most  common  method.  This  method  is  based  on  the 
generalization  determined  empirically  for  Newtonian  oils  that 
the  viscosity-temperature  relationship  is  linear  on  the  ASTM 
temperature  Saybolt  viscosity  paper.  Then  the  viscosity 
values  determined  by  the  Saybolt  or  Ostwald  capillary  visco¬ 
meters  at  the  higher  temperatures  can  be  used  to  determine 
the  viscosity  values  at  the  lower  temperatures.  This  method 
is  limited  to  the  Newtonian  paraffinic  crude  oils  above  the 
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transition  temperature,  for  these  waxy  crude  oils  exhibit 
pseiidoplastic  properties  below  the  transition  temperature 
where  the  linear  relationship  is  not  applicable. 

The  value  of  the  apparent  viscosity  for  these 
waxy  crudes  as  determined  at  the  low  shear  rates  present 
in  these  precision  capillary  viscometers  is  much  higher  than 
the  apparent  viscosity  value  obtained  under  the  high  rates 
of  shear  in  commercial  pipelines.  The  over-designing  of 
commercial  pipelines  may  result  from  the  use  of  these 
values  of  the  apparent  viscosity.  By  duplicating  either 
the  Reynolds  number  and/or  the  shear  rates  of  commercial 
applications  on  test  pipelines,  the  apparent  viscosity 
may  be  determined  under  flow  conditions  dynamically 
Similar  to  those  present  in  commercial  pipelines.  The 
value  of  the  apparent  viscosity  undergoes  the  largest 
change  over  the  low  range  of  shear  rates,  so  that  the 
reproducibility  of  the  various  precision  capillary  visco¬ 
meters  will  be  poor,  unless  correlated  with  the  rate  of 
shear.  The  apparent  viscosity-shear  rate  relationship  may 
be  determined  with  precision  over  the  more  critical  low 
range  of  shear  rate  by  the  use  of  the  rotational  concentric 
cylinder  viscometer. 

A  research  program  with  the  above  objectives 
was  initiated  by  the  Interprovincial  Pipe  Line  Co.  in  1950. 
Under  the  direction  of  Professor  G-.  W.  G-ovier,  W.  Sidjak  did 
the  preliminary  designing,  and  M.  Chmilar  constructed  and 
tested  the  model  or  pilot  pipeline,  to  determine  the  friction 
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fact or- Reynolds  number  relationship.  An  investigation  of 
several  Alberta  crudes  with  the  modified  MacMichael,  indicated 
that  most  of  these  crudes  are  pseudoplastic  at  low  rates  of 
shear.  The  pilot  pipeline  has  been  used  for  the  investi¬ 
gation  of  five  different  blends  and  crudes  at  rates  of  shear 
which  may  have  been  too  high  to  detect  any  pseudoplastic 
behavior.  To  investigate  the  crude  oil  flow  behavior,  at 
the  low  rates  of  shear  where  the  pseudoplastic  effect  is 
more  prevalent,  as  well  as  check  the  accuracy  of  the  pilot 
pipeline,  two  rotational  concentric  cylinder  viscometers  were 
used.  One  was  the  modified  MacMichael  viscometer  used 
previously  by  Sidjak  and  the  other  was  the  Alberta  visco¬ 
meter,  constructed  by  Villiers  Fisher,  under  the  direction 
of  Professor  G.  W.  Govier. 
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REVIEW  OF  THEORY 

General  Classification  and  Flow  Behavior  of  Liquids 

It  is  necessary  to  have  a  complete  understanding 
of  the  flow  properties  of  all  classes  of  liquids  in  order 
to  fully  understand  and  interpret  the  data  obtained  on  the 
crude  oils  Investigated.  The  following  is  a  brief  discuss¬ 
ion  of  the  more  important  liquids  and  the  characteristics  of 
their  flow  behavior  which  may  be  compared  to  those  observed 
in  the  crude  oils  tested. 

There  are  six  general  classifications  of  liquids, 
depending  on  the  relationship  of  the  shear  stress  to  shear 
rate  in  laminar  flow.  These  classifications  are  as  follows: 
Newtonian;  pseudoplastic;  Bingham  plastic  (or  as  they  are 
more  often  referred  to,  plastic);  thixotropic;  dilatant;  and 
rheopectic,  A  more  detailed  discussion  of  all  these  liquids 
appears  in  Appendix  A.  The  more  important  classifications 
with  regard  to  this  investigation  are  the  Newtonian,  the 
pseudoplastic,  and  the  thixotropic  liquids. 

Newtonian  liquids  are  characterized  by  a  constant 
ratio  of  the  shear  stress  to  the  shear  rate  in  laminar  flow, 
and  this  constant  is  defined  as  the  coefficient  of  viscosity, 
or  more  commonly,  viscosity.  For  pseudoplastic  liquids  the 
ratio  of  the  shear  stress  to  the  shear  rate  is  a  non-linear 
function  of  the  shear  rate,  and  this  ratio,  at  any  shear  rate 
is  defined  as  the  apparent  viscosity.  The  plot  of  the  varia- 
bles-shear  stress  and  apparent  viscosity,  over  a  range  of 
shear  rates  is  called  a  consistency  curve.  General  equations 
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have  been  derived  from  the  consistency  curve  to  relate  the 
variables  for  pseudoplastic  liquids.  The  apparent  viscosity 
of  a  thixotropic  liquid  is  a  function  of  the  shear  rate 
history  and  the  shear  rate  at  the  instant  of  measurement. 

For  both  thixotropic  plastic  and  pseudoplastic  liquids,  the 
upcurve  and  downcurve  on  a  consistency  plot  do  not  coincide, 
with  the  result  that  a  hystereis  loop  is  formed.  The  size 
of  the  hysteresis  loop  indicates  the  degree  of  thixotropy  (12). 

General  Equations  of  Flow  for  Liquids  in  Pipes 


The  pressure  drop,  AP,  experienced  by  a  Newtonian 


liquid  in  flowing  a  finite  distance  in  a  level  pipe  of 
constant  cross  section  under  steady  state  conditions,  can  be 
described  by  the  Darcy-Weisbach  equation  (27,9)  : 


AP 


(1) 


whe  re 


2 

AP  -  pressure  drop,  lb.^/in. 
f  -  friction  factor 


L  »  length  of  section,  ft. 
p  =  density  of  flowing  liquid,  lb./ft.^ 

V  =  bulk  velocity,  ft. /sec. 

D  s  pipe  diameter,  ft. 

2 

Sc  2  conversion  constant  32.2  (  lb. ft.  )/(lb.fsec.  ) 


With  the  appropriate  friction  factor  from  Charts 


(27,32),  or  equations  (6),  the  above  equation  applies  to  the 
flow  of  Newtonian  liquids  in  the  laminar,  the  transition 
and  the  turbulent  flow  regions. 
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Another  relationship  which  is  a  fundamental  lav/  of 


liquid  flow,  is  the  linear  shear  stress  distribution  from  a 
maximum  at  the  pipe  wall  to  zero  at  the  axis  (3).  This  is 
expressed  as  follows: 


Tr 

where 

\ 

H3 

11 

r  /  R 

(2) 

Tw 

- 

shear  stress  at  the  wall, 

lb, f /in, 2 

Tr 

- 

shear  stress  at  any  radius 

r,  lb,f./in,2 

r 

SB 

any  radius,  ft. 

R 

- 

radius  of  the  pipe,  ft. 

For  a  pipe,  it  may  readily  be  shown 

that  the  shear 

stress  at  the  wall  for  any  type  of  flow  is  related  to  the 
pressure  drop  and  the  dimensions  of  the  system  by  (9,10)  : 

Tw  -  APD  /  4L  (3) 

Laminar  Flow  of  Newtonian  Liquids 

From  the  relationship  of  the  shear  stress  to  the 
shear  rate  for  Newtonian  liquids  (1)  : 

T  =  ^ jjl  (dv  /  dr)  /  gc  (4) 

where 

=  coefficient  of  viscosity, or  viscosity  lb. /ft 
(dv/dr)  a  shear  rate,  rate  of  shear,  or  velocity 
gradient,  sec,”“ 

It  may  be  shown  that  the  velocity  profile  in  fully 
developed  laminar  flow  in  a  pipe,  is  a  parabola.  The  shear 
rate  at  the  wall  may  be  derived  to  obtain  (1)  : 

(dv/dr )w 


sec . 


-  8V  /  D 


(5) 
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or  substituting  in  the  expression  for  the  shear  stress  at  the 
wall  of  the  pipe,  Equation  (3),  into  Equation  (4)  t 

(dv  /  dr)w  »  ~AP  Dcj/  4L/a/  (6) 

In  the  laminar  flow  region,  it  has  been  shown  by 
the  rigorous  derivation  of  the  Hagen-Poiseuille  Law  (27)  and 
by  experimental  data  (32),  that  the  friction  factor  is  : 

f  ■  64  /  NRe  (7) 

where 


Nrq  -  Reynolds  number 

=  VDp/yU'  dimensionless 
Thus  when  the  expression  for  f  is  substituted 
into  the  Darcy-Weisbach  equation,  the  familiar  Hagen- 
Poiseuille  law  is  obtained 


AP=  32 j*  LV/g0 


(8) 


Turbulent  Flow  of  Newtonian  Liquids 

In  the  turbulent  flow  region,  the  shear  stress  is 
the  summation  of  a  viscous  agent  and  a  momentum  exchange  agent. 
The  turbulent  boundary  layer  consists  of  three  regions  of 
flow:  the  laminar  sublayer;  the  buffer  layer;  and  the 
turbulent  core.  In  the  laminar  sublayer  the  viscous  forces 
control  the  velocity  profile,  and  in  the  turbulent  core  the 
momentum  exchange  agent  controls  the  velocity  profile.  In 
the  transition  zone  both  agents  are  variably  effective. 

Because  the  velocity  profile  is  the  result  of  the  varying 
effect  of  these  two  agents  across  the  radius  of  the  pipe, 
it  is  necessary  to  derive  individual  equations  for  the 
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three  regions  of  flow.  These  equations  were  derived  by  von 
Karman  and  Prandtl  (4)  and  the  result  may  be  represented 
graphically  as  the  Universal  Velocity  Profile  (4),  Fig. (3^) 
Appendix  A.  An  examination  of  this  profile  shows  that  the 
velocity  gradient  is  high  in  the  laminar  sublayer,  but 
relatively  low  in  the  turbulent  core.  These  equations, 
however,  do  not  take  into  account  the  shear  rates  in  all 
other  directions,  but  only  the  net  shear  rates  in  the 
direction  of  streaming  (27).  The  value  of  the  shear  rate 
determined  in  turbulent  flow  is  therefore  only  completely 
comparable  with  that  in  laminar  flow  in  the  laminar  sublayer 
itself.  The  shear  rate  at  the  wall  and  therefore  within  the 
laminar  sublayer  is  given  by  equation  (6). 

The  equation  for  the  average  velocity  in  turbulent 
motion,  can  be  determined  from  the  integration  of  the 
velocity  profile  (4).  From  this  equation  von  Karman  develop*, 
ed  the  friction  factor  equation  for  fully  developed  isothermal 
turbulent  flow  of  Newtonian  liquids  (3)*  This  equation  is 
applicable  to  smooth  tubes,  from  a  Reynolds  number  of  4,000 
to  3,000,000:; 

1  /  x/T  s  2.0  log  NReVT  •  0.8  (9) 

This  equation  involves  a  trial  and  error  solution 
for  the  friction  factor.  Many  other  equations  (strictly 
empirical)  have  been  derived  which  are  more  easily  applied. 

One  empirical  equation  based  on  data  obtained  on  tubes,  is 
that  of  Blasius,  which  is  valid  from  a  Reynolds  number  of 
3,000  to  100,000  (9,10)  : 


•  -  '  '  ‘  -i  -  i 

L 

- 

: 

.  1  ■'  vi  • 

* 

. 

. 

n:  . ; 


* 

•  •  '  0  ■ 

\  ’  -  ' 


. 


.  .  ■ 

■ 

.'"'f  ‘'-i-  : 


, 

; 

.  -  . 

. 

■i  \  1^;  -i 

1&  ‘I  <.-■  •" '  -1 

(  t  '  « 


10 


f  ,  0.316  /  (NRe)i  (10) 

From  data  obtained  on  crude  oil  flowing  through 
6  to  12  inch  pipelines,  the  Heltzel  equations  were  determined 
for  two  ranges  of  Reynolds  numbers.  For  the  range  from  2,500 
to  57,600,  the  friction  factor  is  (2,14)  : 

f  =  (0.364)  f  (t)R@ )0.265  (ii) 

For  the  range  above  57,600  the  friction  factor  isj 

f  =  (0.157)  /  (NRe)0-188  (12) 

These  equations  have  found  wide  application  in  the 
design  of  commercial  pipelines,  especially  for  the  handling 
of  crude  oil  and  petroleum  products. 

Laminar  and  Turbulent  Flow  Equations  For  Pseudoplastic  Liquids 
There  has  been  no  general  equation  determined  which 
completely  describes  the  laminar  flow  behavior  of  pseudoplastic 
liquids  (12).  As  a  result,  there  are  two  general  methods  for 
applying  laboratory  data  to  flow  problems. 

One  method  is  the  integration  of  the  shear  stress- 
shear  rate  functional  relationship  over  a  limited  range  of 
shear  rates  in  a  pipe.  The  data  obtained  on  a  rotational 
concentric  cylinder  viscometer  is  used  to  obtain  a  consistency 
curve  for  the  pseudoplastic  liquid.  The  functional  relation¬ 
ship  of  the  shear  stress  T,  and  shear  rate  dv/dr,  may  be 
expressed  as  follows  (1)  : 

(dv/dr^  z  j  (aPD  /  4L)  (13) 


where 


j  denotes  "function  of" 
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Equation  (13),  usually  of  the  parabolic  type,  may 
only  apply  for  a  certain  range  shear  rates.  This  function 
is  substituted  into  an  equation  relating  the  shear  rate  at 
the  wall  to  the  shear  stress  at  the  wall  which  is  independent 
of  the  shape  of  the  consistency  curve.  An  equation,  developed 
by  Rabinowitsch  (1),  which  is  independent  of  the  consistency 
curve,  is  : 

(d  (dv)  ) 

- (dv/dr )  «  3 (dv/dr )  /  T  (  (dr)  _ )  (14) 

(  dT  ) 

Substituting  into  this  equation  the  derived 
functional  relationship  equation  (13),  and  integrating 
gives  the  shear  stress  and  shear  rate  relationship  for  any 
pipe  over  the  range  of  shear  rates  of  equation  (13). 

The  other  more  direct,  but  less  rigorous  method  of 
treatment  for  the  data  obtained  by  viscometers  or  pilot 
pipelines,  is  based  on  an  apparent  viscosity  as  defined  by 
the  Hagen-Poiseuille  equation  for  laminar  flow,  or  the  Darcy- 
Weisbach  equation  (with  an  appropriate  friction  factor)  for 
turbulent  flow. 

There  are  several  ways  in  which  an  "apparent  viscosity" 
may  be  defined  for  a  pseudoplastic  liquid.  One  definition 
for  the  apparent  viscosity  is  based  on  the  laminar  flow  of 

the  liquid  under  a  constant  rate  of  shear  over  the  cross 

* 

section,  or  at  a  point  where  the  shear  stress  and  rate  of  shear 
are  known.  Then 

“  -  Tgc  /  (dv/dr)  (15) 

tf 


where 


apparent  viscosity  lb. /ft. sec 
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The  condition  of  a  constant  rate  of  shear  is  dupli¬ 
cated  approximately  in  the  rotational  viscometer  by  making 
the  cup  and  bob  radii  very  nearly  equal.  The  apparent 
viscosity  is  then  obtained  by  the  following  relationship  for 
.  the  rotational  viscometer, 

--  K  i  /  he  N  (16) 

where 

=  apparent  viscosity,  lb. /ft. sec. 

K  =  combined  instrument  and  torsion  constant  (See 
Appendix  C  for  units) 

a  net  deflection  of  the  torsion  wire  in  degrees 
s  rotary  speed  of  the  cup  l/sec. 

-  equivalent  depth  of  immersion  in.  or  cm. 

In  the  laminar  flow  of  liquids  through  pipes,  the 
apparent  viscosity  may  be  defined  by  the  Hagen-Poiseuille  law 
(equation  8),  and  correlated  with  the  shear  rate  at  the  wall 
(equation  6).  Such  a  correlation  may  be  expressed  in  equation 
form  (1),  or  may  be  expressed  graphically  for  quick  reference 
in  the  solution  of  flow  problems  for  other  systems  (12).. 

In  turbulent  flow  the  apparent  viscosity  is 
defined  as  that  viscosity  which  satisfies  both  the  equation 
4P  =  f  L  V2  p  /  2  Dgc 

and  an  appropriate  relationship  for  the  friction  factor  (5,27), 
such  as  the  Heltzel  equation  for  the  low  Reynolds  numbers. 
Rearranged  to  express  the  apparent  viscosity  explicitly 
/=  v  D  P  /  NRe 

V  D  p/(0*364/f  )2-77 


3, 000<NRe<  57,600  (17) 
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It  has  been  determined  empirically  that  the  value 
of  the  apparent  viscosity  as  defined  above  is  very  nearly  a 
constant  for  all  flow  rates  in  turbulent  flow  (1,36). 

Accordingly,  one  value  determined  in  turbulent  flow  in  a 
pilot  pipeline  or  at  a  high  rate  of  shear  suffices  for  the 
solution  of  flow  problems  in  other  systems.  This  is  supported 
by  the  theory  of  turbulent  flow  (see  Appendix  A),  in  which 
it  is  shown  that  the  viscous  effects  are  confined  to  the  thin 
laminar  sublayer  and  part  of  the  buffer  layer,  so  that  the 
pressure  drop  is  very  nearly  independent  of  the  viscosity. 

Effect  of  a  Temperature  Gradient  on  Newtonian  Liquids  in  Pipes 

The  standard  charts  and  equations  for  the  friction 
factor-Reynolds  number  relationship  have  been  determined  for 
the  isothermal  flow  of  Newtonian  liquids.  If  there  is  a  tempera¬ 
ture  drop  from  the  flowing  liquid  to  the  pipe  wall,  the  value 
of  the  friction  factor  is  increased  (10,15).  The  change  in 
the  friction  factor  is  a  function,  of'  the  temperature  drop, 
the  effect  of  temperature  on  the  liquid  viscosity,  and  whether 
the  flow  is  laminar  or  turbulent.  The  effect  of  the  tempera¬ 
ture  gradient  is  more  serious  in  laminar  than  in  turbulent 
flow,  because  of  the  distorted  velocity  profile  (10,15)*  and 
larger  temperature  gradients  present.  The  temperature 
gradient  is  also  more  serious  with  higher  viscosity  oils  than 
lighter  oils,  as  a  result  of  the  variation  of  viscosity  with 
temperature.  If  the  isothermal  friction  factor  is  used  for  a 
flow  condition  where  the  oil  is  being  cooled,  the  value  of  the 
viscosity  or  apparent  viscosity  will  be  high.  Several 
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empirical  plots  and  equations  have  been  obtained  to  express 
the  effect  of  a  temperature  gradient  on  the  friction  factor 
(15,29),  and  the  most  applicable  equation  is  that  of  Keevil 
and  McAdams  (see  Appendix  A). 

Flow  of  Newtonian  and  Pseudoplastic  Liquids  in  the 
Concentric  Cylinder  Viscometer 

The  equation  of  flow  for  Newtonian  liquids  in  a 

rotational  concentric  cylinder  system  has  been  rigorously 

derived  (12,13,17).  The  simplified  form  of  this  equation 

which  is  used  for  computing  the  viscosity,  and  apparent 

viscosity,  for  Newtonian  and  pseudoplastic  liquids  respective- 

iy,  is 

>4  *  <t>  K  /  he  N 

where 

/  »  net  deflection  in  degrees 

for  the  Alberta  viscometer  360°A=  2TT 
for  the  MacMichael  viscometer  300°M»  2  Tr 
K  =  v^ire  constant  for  one  length,  one  diameter  of 
torsion  wire 

for  the  Alberta,  dimensions  are  (in. lb ./ft .  sec.2°A) 
for  the  MacMichael,  dimensions  are  (cm.lb./ft . sec?°M) 
The  application  of  the  Newtonian  liquid  relationships 
to  pseudoplastic  liquids,  is  not  serious  for  slightly  pseudo¬ 
plastic  liquids,  because  the  variation  of  the  shear  rate  across 
the  annulus  is  small. 

For  the  pseudoplastic  liquids,  it  is  necessary  to 
correlate  the  apparent  viscosity  with  the  shear  rate  (1,3). 
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Since  the  average  weighted  shear  rate  for  Newtonian  liquids 
is  the  value  at  the  midpoint,  this  value  has  been  used  more 
widely  than  the  value  at  the  bob  wall  (34-).  The  shear  rate 
at  any  radius  r  between  Rj  and  Rg  may  be  computed  with  the 
following  relationship, 


(dv/dr)r  -  4  Ntt(r2  r2  ^ 

r2  (R2”Ri  ) 


(18) 


where 


(dv/dr)r  »  rate  of  shear  at  radius  r  sec. 
r  =  any  radius,  in. 

R]_  *  radius  of  bob,  in. 

P-2  -  radius  of  the  cup,  in. 

It  may  be  shown  that  a  smaller  change  in  shear  rate 
across  the  annulus  will  be  obtained,  by  making  R^  very  nearly 
equal  to  R2.  Equations  have  been  developed  for  the  flow  of 
pseudoplastic  liquids  in  this  system,  and  these  will  be 
discussed  in  the  appendix. 
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REVIEW  OF  PREVIOUS  WORK 

Several  investigators  who  have  experimentally 
determined  the  rheological  properties  of  certain  paraffin 
base  crude  oils  at  temperatures  below  the  pour  point  (11,33) 
showed  that  the  flow  rates  under  constant  pressure,  first 
increased  and  then  changed  reversibly  with  change  of  pressure, 
as  is  characteristic  of  thixotropic  liquids.  Some  of  the 
crudes  in  the  Wyoming-Utah  area  have  exhibited  thixotropic 
behavior  at  temperatures  as  high  as  70°F  (2).  The  pseudo¬ 
plastic  properties  of  the  waxy  paraffin  crudes  were  believed 
to  be  a  result  of  the  vaseline  fractions  present,  and  not 
an  effect  caused  by  the  recrystallization  of  wax  (33).  After 
the  addition  of  the  vaseline  fractions  to  the  crude  exhibiting 
Newtonian  behavior,  the  oil  became  pseudoplastic.  Without  the 
vaseline  fractions,  the  crude  oil  exhibited  Newtonian  behavior 
with  wax  crystals  in  solution.. 

Work  performed  on  a  Mexican  oil  showed  that  the 
thermal  history  effected  the  value  of  the  apparent  viscosity 
for  low  temperatures.  This  work  indicated  that  the  long 
molecules,  or  the  formation  of  crystals,  which  takes  place 
at  very  cold  temperatures,  increased  the  apparent  viscosity, 
but  the  effect  disappeared  w ith  time,  if  the  temperature  was 
increased  to  a,  warmer  one  below  the  pour  point  (11).  It  has 
been  determined  for  some  waxy  paraffinic  crudes  that  the 
transition  temperature  from  non-Newtonian  to  Newtonian 
characteristics,  and  the  nature  of  the  apparent  viscosity 
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temperature  curve  are  functions  of  the  thermal  history  (5). 

The  apparent  viscosity  of  a  Newtonian  lubricating 
oil  with  varying  amounts  of  soap  added,  was  determined  for 
shear  rates  varying  from  0.1  to  above  100,000  sec,”1.  The 
results  of  these  tests  in  tubes  showed  that  the  change  in 
the  apparent  viscosity  with  an  increase  in  the  shear  rate  was 
small  above  a  shear  rate  of  1,000  sec."'1  for  the  grease  with 
an  apparent  viscosity  of  2-3  poises.  With  the  heavier  greases 
the  shear  rate  at  which  the  apparent  viscosity  was  approaching 
a  constant  value  of  4<*5  poises  was  between  100,000  to 
1,000,000  sec .“1 (3 ) . 

Flow  tests  conducted  on  asphaltic  bitumens  using 
capillaries  of  varying  lengths  indicated  that  a  lower  overall 
apparent  viscosity  was  obtained  on  the  longer  tubes  for  the 
same  rate  of  shear.  This  method  of  determining  whether  a  liquid 
was  thixotropic  has  also  been  used  by  Alves  et  al  (1)  to  show 
that  both  length  and  smaller  pipe  diameters  reduced  the 
apparent  viscosity.  Pumpability  tests  have  been  conducted  on 
the  Pembina  crude,  which  indicates  that  the  apparent  viscosity 
is  a  function  of  both  the  shear  rate  history,  and  the 
instantaneous  shear  rate  (23).  These  tests  showed  that  the 
value  of  the  apparent  viscosity,  after  passing  through  a 
centrifugal  pump  was  much  lower  than  that  determined  while 
using  a  reciprocating  pump.  There  was  no  indication  of  any 
yield  value  when  starting  the  liquid  to  flow.  The  starting 
pressure  was  high  for  a  short  interval  only,  indicating  that 
the  fluid  may  exhibit  thixotropic  breakdown  until  the  apparent 
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viscosity  attains  an  equilibrium  value  for  that  rate  of  shear 

The  friction  factors  used  for  the  design  of  large 
diameter  pipelines  have  been  obtained  from  the  Stanton  and 
Pannell  data,  determined  by  tests  on  small  diameter  tubes, 
using  air  and  water  to  obtain  a  wide  range  in  the  Reynolds 
number.  In  1929,  Heltzel  correlated  data  obtained  on  6  inch 
to  12  inch  diameter  pipelines  (the  largest  at  that  time) 
transporting  crude  oil,  with  the  data  of  the  above  investi¬ 
gators  (2),  and  determined  that  all  the  friction  factor  data 
could  be  very  nearly  represented  by  two  curves  intersecting 
at  a  Reynolds  number  of  57,600.  The  two  equations  for 
these  lines  (equations  11  and  12)  may  be  used  to  determine 
the  friction  factor  from  a  Reynolds  number  of  2,500  to  above 
57,600. 

Recently,  investigations  have  been  made  to  deter¬ 
mine  the  validity  of  using  the  friction  factor  data  obtained 
on  tubes  and  pipes  up  to  12  inches  in  diameters  as  above, 
for  16  inch  to  30  inch  pipelines.  The  data,  obtained  over  a 
wide  range  of  Reynolds  numbers,  on  pipelines  varying  from 
16  inch  to  30  inches  in  diameter  have  confirmed  the  validity 
of  both  the  Stanton  and  Pannell,  and  the  Heltzel  data.  Since 
the  curves  representing  the  Stanton  and  Pannell  data  have  an 
inflection  point,  it  has  been  expedient  to  draw  a  best  fit 
straight  line  through  the  curve  and  determine  one  equation 
which  is  applicable  from  a  Reynolds  of  6,000  to  130,000.  Thi 
method  has  been  sufficiently  accurate  for  the  design  of 
commercial  pipelines  (2). 
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As  a  result  of  the  improved  mechanical  properties 
of  steel  pipe,  it  has  been  economical  to  increase  through¬ 
puts  by  increasing  the  input  pressure.  The  large  pressures 
in  use  resulted  in  a  substantial  increase  in  the  viscosity. 

This  effect  of  pressure  was  investigated  on  a  rolling  ball 
viscometer  and  the  increase  in  viscosity  per  1000  psi.  for 
a  2 6.6°  API  crude  varied  from  10%  at  69°F  to  17%  at  33.5°F. 

For  a  crude  of  41. 4° API,  the  increase  in  viscosity  varied 
from  13 %  at  69°F  to  42^  per  1000  psi.  at  33.5°F(2). 

Since  the  research  program  was  started  at  the 
University  of  Alberta  in  1950,  there  have  been  three 
investigations  into  Alberta  and  Saskatchewan  crude  oils. 

During  1950,  Sidjak  investigated  several  Alberta  crude  oils 
with  the  modified  MacMichael  (30).  The  conclusions  of  this 
investigation  were  that  the  Redwater  and  several  other  Alberta 
crude  oils  were  pseudoplastic  below  50°F  for  the  range  of 
shear  rates  from  4  to  15  sec."\  and  the  apparent  viscosity 
was  a  function  of  the  thermal  history.  Since  this  investi¬ 
gation,  it  has  been  determined  that  the  method  of  obtaining 
the  average  oil  temperature  was  not  correct,  and  that  no 
precautions  were  taken  to  eliminate  weathering  of  the  crude 
oil  samples,  while  determining  the  effect  of  thermal  history. 

During  the  next  year,  the  model  pipeline  was 
constructed,  and  the  friction  fact or- Reynolds  number  relation¬ 
ship  determined  by  running  stove  oil  through  the  three  test 
lines.  The  friction  factors  determined  by  these  tests  were 
found  to  best  fit  the  curve  represented  by  the  Heltzel  equation(6). 
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During  the  next  year,  Masuda  used  the  pilot  pipeline  to 
investigate  several  Alberta  crude  oils,  with  the  emphasis  on 
the  Redwater  crude  (19).  This  investigation  indicated  that 
all  the  crudes  investigated  were  Newtonian  over  shear  rate 
ranges  from  a  low  of  700  to  1,500  sec.*'1  to  rates  above 
20,000  see.”’1.  The  values  of  the  apparent  viscosity 
obtained  in  laminar  flow  checked  with  the  values  obtained  in 
turbulent  flow,  which  indicates  that  the  apparent  viscosity 
as  defined  by  the  Hagen- Poiseuille  equation  may  be  correlated 
with  the  value  as  defined  by  the  Darcy-Weisbach  equation,  and 
the  Heltzel  friction  factor  relationship. 
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DESCRIPTION  AND  OPERATION  OF  THE 
EXPERIMENTAL  EQUIPMENT 

The  equipment  used  in  this  investigation  consisted 
of  two  main  units:  the  pilot  pipeline  with  temperature 
control  equipment,  and  the  Alberta  and  MacMichael  visco¬ 
meters  with  constant  temperature  baths.  The  pilot  pipeline 
was  located  and  operated  separately  from  the  rotational 
viscometers.  The  concentric  cylinder  viscometers  were 
connected  to  the  same  temperature  bath,  and  were  often 
operated  together. 

Pilot  Pipeline 

The  pilot  pipeline  consists  of  many  components 
to  control  and  measure  pressure,  temperature,  and  flow  rates. 
The  pumping  system  consists  of  two  rotary  gear  pumps  and 
two  valve  manifolds  for  controlling  pressure  and  direction 
of  flow.  To  control  temperature,  a  refrigerated  bath 
surrounds  the  pre-cool  lines.  The  flow  rates  are  controlled 
by  two  sets  of  by-pass  valves  from  the  pump,  and  a  throttling 
valve  at  the  end  of  the  test  section.  There  are  both  a 
heavy  duty  and  a  precision  pressure  gauge  at  the  end  of  each 
section.  A  line  diagram  of  the  pilot  pipeline  is  presented 
in  Fig.(l)  with  the  important  components  numbered  and  listed 
in  Table  (1). 

The  pumping  system  is  designed  to  permit  a  wide 
variation  of  the  flow  rates  and  pressures.  The  high  capacitie 
and  pressures  are  developed  by  a  12  to  16  g.p.m.  1000  psi0 
continuous  operation  Rotary  gear  pump,  driven  by  a  15  H.P. 
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Lincoln  synchronous  motor.  For  small  rates  of  flow,  a 
Pasco  rotary  gear  pump  capable  of  1  l/2  g.p.m.  and  1000  psi. 
is  driven  by  a  2  H.F.  motor,  3oth  these  pumps  are 
exhibited  in  Fig,  3  and  also  the  common  suction  lines  and 
discharge  lines.  To  control  the  pressure  delivered  to  the 
system,  one  3/4  inch  globe  and  one  l/8  inch  needle  valve 
are  used  to  by-pass  from  the  pumps  back  to  the  resevoir 
tank.  After  leaving  the  pumps,  the  oil  is  passed  through 
a  pre-cool  line  of  l/2  inch  Schedule  80  pipe  before  being 
either  by-passed  again,  or  directed  into  one  of  the  test 
lines . 

After  passing  the  pre-cool  section,  the  oil  tempera¬ 
ture  is  measured  by  a  Weston  high  pressure  -40°F  to  120°F 
bimetallic  thermometer,  and  then  directed  into  any  or  all 
of  the  test  lines.  Three  Nordstrom  plug  valves  are  used  to 
direct  the  oil  into  one  of  the  following  stainless  steel 
tubing  lines  of  nominal  size,  3/8",  l/2",  or  5/8",  This 
equipment,  located  outside  the  bath,  is  partially  insulated 
with  glass  wool  and  aluminum  vapour  proofing, 

A  test  section  consists  of  8.5  feet  for  a  calming 
length,  and  19.5  feet  as  a  test  section.  The  bath  consists 
of  two  sections,  one  for  a  calming  length,  the  other  for  the 
test  lines,  connected  by  liquid  and  vapour  equilibrium  lines 
as  in  Fig. (2).  The  upstream  pressure  taps,  and  valve  mani¬ 
fold  are  located  between  the  two  sections.  Approximately 
18  l/2  feet  of  the  test  lines  are  in  contact  with  the  liquid 
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At  the  downstream  end  of  the  test  section,  the 
high  pressure  quick  shut  Nordstrom  valves  can  be  used  to 
choke  the  flow,  thereby  increasing  the  average  pressure  in 
the  test  line.  A  short  distance  from  the  pressure  tap,  the 
oil  temperature  is  obtained  by  a  Weston  thermometer. 

For  large  flow  rates,  the  calibrated  60  gallon 
tanks  are  used,  and  for  small  rates  of  flow,  the  6  gallon 
tank  is  used.  The  valve  manifold  at  the  tanks  enables  the 
oil  to  be  metered,  recirculated,  or  withdrawn.  The  tanks 
are  sealed,  and  are  capable  of  holding  1  or  2  psi.  for  long 
periods  of  time,  thus  reducing  weathering. 

The  pressure  system  is  designed  to  measure  differ¬ 
ential  pressures  and  to  protect  the  equipment  from  high 
pressures.  A  2000  psi.  gauge  is  used  at  the  pumps  prior 
to  the  by-pass  manifold.  To  protect  both  of  the  Heise  gauges 
and  also  the  downstream  Acragage,  a  5000  psi.  and  a  2000  psi. 
gauge  respectively  are  exposed  to  line  pressure.  On  lbw  press¬ 
ure  tests,  the  upstream  pressure  is  determined  with  a  500 
division  1000  psi.  Heise  gauge,  and  the  downstream  pressure 
is  measured  with  a  60  division  60  psi.  Acragage.  For  small 
flow  rates,  a  60  inch  oil  over  mercury  manometer  is  used.  The 
pressure  gauges  are  recalibrated  after  adverse  treatment,  or 
when  warranted. 

The  temperature  of  the  sealed,  insulated  bath  of 
liquid  and  gaseous  Freon  12,  is  controlled  by  a  3  H.P. 
Frigidaire  compressor,  and  strip  heaters.  Initially,  the 
compressed  Freon  was  expanded  into  l/2"  copper  tubing  in  the 
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vapour  phase  of  the  bath.  Hov/ever,  much  of  the  efficiency 
was  lost  at  the  lower  temperature,  so  the  system  was  changed 
to  enable  the  Freon  to  expand  directly  into  the  bath.  Over 
six  hundred  pounds  of  Freon  12  are  necessary  to  cover  the 
test  lines,  the  level  of  which  is  determined  by  a  Jorgenson 
sight  glass  illustrated  in  Fig.  7.>  The  bath  is  insulated 
with  approximately  2  inches  of  glass  wool  and  aluminum  vapour 
proofing.  Fig.(l)  presents  a  schematic  sketch  of  the  pilot 
pipeline  and  Table  1  presents  a  detailed  list  of  the 
important  components. 
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SCHEMATIC  SKETCH  OF  THE  PILOT  PIPELINE 
CRUDE  OIL  LINES 
PRESSURE  GAUGE  LINES 
REFRIGERATION  LINES 
DESIGNATION  OF  THE  COMPONENTS 
I  TO  16  OIL  LINES 

17  TO  23  PRESSURE  GAUGE  LINES 
25  TO  33  REFRIGERATION  8  TEMPERATURE  LINES 
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TABLE  1 


Description  of  the  Pilot  Pipeline  Fig.  (1) 

Component 

Number  Description 

1  Two  sixty,  and  one  six.  Imperial  gallon  measuring  tanks 
that  have  been  calibrated  with  water. 

2  Strainer 


5 

4 

5 

6 

7 

8 

9 

10 

11  &  15 

12 

14 

15 

16 

17 

18 

19 

20 
21 


12  to  16  I.gpm.  Roper  Rotary  Gear  Pump,  1000  psi.  The 
volume  is  varied  by  changing  the  motor  drive  pulley. 

1  l/2  I.  gpm.  Pasco  Rotary  Gear  Pump  1000  psi. 

2000  psi.  pressure  gage  with  dampening  valve. 

Safety  valve  set  at  1200  psi. 

5/4  inch  Globe  by-pass  valve. 

l/8  inch  Needle  by-pass  valve. 

l/2  inch  Sch.  80  pre-cool  line. 

l/2  inch  Globe  by-pass  valve. 

Inlet  and  outlet  Weston  bimetallic  thermometers 
«-40°F  to  120°F  in  1°F  divisions. 

Downstream  plug  valve  used  for  throttling. 

Valve  header  to  direct  test  volume  into  the  tanks. 

Sight  glasses  and  meter  sticks  for  measuring  volumes. 

Vapor  equilibrium  pipes  and  safety  valve  set  at  5  psi. 

Upstream  pressure  taps  and  valve  assembly. 

Heise  gauge  with  1000  psi.  range  in  2  psi.  divisions. 

5000  psi.  Acragage  gauge. 

60  psi.  Acragage  gauge  in  1  psi.  divisions. 

Heise  gauge  with  1000  psi.  range  in  2  psi.  divisions. 
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22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 


60  inch  Oil  over  Mercury  manometer. 

2000  psi.  gauge  used. 

Glass  wool  insulation  and  Aluminum  vapor  proofing. 

3  H.P.  water-cooled  Frigidaire  compressor,  with 
water  cooler  and  automatic  pressure  controls. 

High  pressure  throttling  valve  to  copper  tubing 
heat  exchanger. 

High  pressure  throttling  valve  direct  to  bath,  and 
liquid  Freon  equilibrium  lines. 

Freon  vapor  discharge  lines  from  bath. 

-30  in.  Hg.  to  150  psig.  bath  pressure  gauge. 

Pressuretrol  automatic  device  for  cutting  out  strip 
heaters  when  bath  pressure  reaches  set  pressure. 

Jorgenson  sight  glass  for  finding  level  of  liquid 
freon. 

Vacuum  line  to  remove  air  from  bath. 

Cooling  water. 

Strip  heaters. 
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The  operation  of  the  pilot  pipeline  consisted  of 
modifying  the  flow  rates  and  bath  temperature  to  obtain  an 
isothermal  and  constant  oil  temperature  throughout  each 
test.  To  detect  pseudoplastic  effects,  the  shear  rate  was 
varied  as  much  as  possible  at  one  temperature.  However, 
the  oil  properties  and  mechanical  construction  prevented 
attaining  the  isothermal,  and  constant  flow  rates  desired. 
The  following  discussion  will  be  concerned  with  a  general 
description  of  the  procedure  used  for  one  crude  oil.  The 
methods  have  changed  with  experience,  and  were  modified 
with  changes  in  the  oil  properties. 

Before  pumping  the  crude  oil  into  the  system,  it 
was  completely  flushed  with  varsol.  After  as  much  varsol 
as  possible  was  removed,  the  system  was  flushed  with  a 
small  quantity  of  the  crude  to  be  tested.  This  crude  was 
then  removed,  and  the  sample  was  pumped  into  the  tanks 
from  barrels,  and  the  system  sealed.  After  completely 
circulating  the  system  for  a  period  of  four  to  six  hours, 
a  sample  was  taken  for  chain  gravimetric  specific  gravity 
determination.  The  gravity  was  determined  on  the  Westphal 
balance  every  5°F  from  80°F  to  40°F  on  two  samples,  one  of 
which  was  cooled  as  readings  were  observed,  and  the  other 
was  heated  from  40°F  as  readings  were  taken  (Appendix  D). 
For  the  heavier  crudes,  only  a  small  amount  of  weathering 
occurred  while  running  the  gravity  determination  tests. 
This  procedure  was  not  satisfactory  for  the  Pembina  crude, 
so  the  gravity  was  determined  with  a  hydrometer  at  several 
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different  temperatures.  The  gravity  determinations  for  all 
the  crudes  are  tabulated  in  Appendix  D. 

After  being  pumped  for  4  to  6  hours,  it  was 
assumed  that  the  oil  was  homogeneous,  and  that  future 
weathering  would  be  negligible.  It  was  possible  to  cool 
the  oil  to  approximately  40°F,  so  for  most  of  the  tests 
above  this  temperature,  the  bath  was  kept  only  cold  enough 
to  prevent  excess  heating  of  the  oil  while  it  was  being 
pumped  through  the  pre-cool  and  test  lines.  Because  only 
the  lower  part  of  the  tanks  were  recycled  and  cooled  while 
attaining  thermal  equilibrium,  the  temperature  rose  during 
a  test  because  there  was  warmer  oil  above  the  recycled  oil. 
The  oil  was  heated  because  of  pipe  friction  as  it  passed 
through  the  pump  and  the  throttle  valve,  and  this  increased 
the  oil  temperature  in  the  tanks  during  a  test.  When  the 
temperature  drop,  or  rise,  from  inlet  to  outlet  was  as 
small  as  possible  for  one  flow  rate,  a  run  was  started.  The 
flow  was  directed  into  a  partially  empty  tank,  and  with  the 
aid  of  a  100  watt  light  for  sighting  the  oil  level  in  the 
glass  gauge,  the  stop  watch  was  started  as  the  oil  passed  one 
of  the  centimeter  marks.  The  observations  were  recorded  as 
follows:  the  upstream  pressure,  downstream  pressure,  down¬ 
stream  temperature,  the  upstream  pressure  and  temperature, 
and  start  again  at  the  upstream  pressure.  The  minimum 
height  of  fluid  measured  was  at  least  30  centimeters  in  any 
of  the  tanks,  and  if  the  temperature  was  constant,  as  much 
as  100  centimeters  of  oil  was  collected.  If  the  temperature 
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rose  over  3°F  or  4°F,  the  test  was  either  stopped  or  used 
for  computations  and  then  possibly  compared  and  usually 
omitted.  Generally  the  tests  were  long  enough  to  observe 
eight  upstream  pressure  readings  and  four  temperatures  at 
each  end.  At  the  end  of  the  test,  the  total  length  of 
time  and  the  final  height  in  the  tank  were  recorded.  The 
temperature  of  the  oil  being  measured  in  the  tanks  was 
observed  at  the  pipe  next  to  the  tanks.  Initially  this 
temperature  was  obtained  at  the  surface  of  the  tanks  (see 
Appendix  E).  All  the  data  of  the  pilot  pipeline  is 
recorded  in  Appendix  E. 

When  the  oil  over  mercury  manometer  was  used, 
there  were  twice  as  many  readings  taken  of  the  manometer 
than  of  either  of  the  temperatures,  as  the  temperature 
changes  at  low  rates  of  flow  were  very  small,  whereas  the 
pressure  changes  gave  a  quick  indication  of  steady  state 
conditions . 

For  tests  at  temperatures  below  40°F,  it  was 
necessary  to  have  the  bath  as  cold  as  possible  and  the 
flow  rates  low  enough  to  attain  the  desired  test  tempera** 
ture.  The  method  of  obtaining  the  observations  during  a 
test  was  not  changed  during  the  low  temperature  runs. 

With  the  Pembina  crudes,  it  was  found  that  the  flow  rates 
at  low  rates  of  flow  often  were  not  constant,  because  of 
the  nature  of  this  crude  oil,  and  the  method  of  bypassing. 
The  tests  were  checked  for  this  behavior  by  plotting  oil 
levels  versus  time  and  any  large  discrepancies  could  be 
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detected.  The  weathering  of  the  crude  oils  while  making  the 
runs  was  reduced  during  the  later  tests  by  running  the  low 
temperature  runs  first,  and  the  warmer  tests  last. 

To  obtain  the  effect  of  pressure  on  the  viscosity, 
the  Nordstrom  valve  on  the  downstream  side  of  the  test  lines 
was  used  to  throttle  the  oil  flow,  and  regulate  the  down¬ 
stream  pressure.  The  high,  medium,  and  low  pressure  runs 
we re  intermixed  so  that  better  correlation  of  the  data 
could  be  obtained,  as  well  as  to  detect  any  weathering. 
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Alberta  Rotational  Cylinder  Viscometer 

This  instrument,  which  was  designed  and  constructed 
by  Villiers-Fisher,  is  described  in  detail  in  his  the  sis  (30). 
A  brief  resume  of  this  description,  and  the  changes  made, 
will  be  enclosed  here.  The  Alberta  rotational  viscometer 
consists  of  a  cylindrical  cup  gear  connected  to  a  variable 
speed  drive.  There  are  a  set  of  torsion  wires,  with  each 
one  marked,  so  that  the  wires  can  be  set  in  the  instrument 
10.5  inches  from  chuck  top  to  chuck  top.  The  depth  of 
immersion  of  the  bob  in  the  cup  is  controlled  by  sliding 
the  elevator  assembly  up  or  down  until  the  grooved  notches 
are  lined  up  with  the  elevator  mount.  The  cup  is  rotated 
from  0  to  250  rpm.  by  a  G-rahaio  variable  speed  transmission, 
driven  by  a  l/8  H.P.  motor.  A  revolution  counter,  mounted 
on  the  cup  shaft,  is  used  to  count  the  revolutions  over  a 
length  of  time,  usually  thirty  seconds  to  one  minute. 

The  temperature  of  the  circulating  stream  is  measured  at  the 
inlet  and  outlet  from  the  jacket  surrounding  the  rotating 
cup.  The  oil  can  be  forced  into,  or  removed  by  suction 
from  the  cup,  through  a  hollow  shaft  at  the  bottom  of  the 
cup.  This  machine  was  designed  and  constructed  to  close 
tolerances,  so  that  it  is  possible  to  reproduce  deflection 
readings  to  within  3/£  and  speed  readings  above  0,5  rps.  to  1%, 
The  crude  samples  used  in  both  the  rotational 
viscometers  were  obtained  from  the  tanks  of  the  pilot  pipe¬ 
line  at  the  completion  of  the  tests  for  the  apparent  viscosity 
curve,  and  stored  in  sealed  five  gallon  tins.  The  oil  was 
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TABLE  2 


Description  of  the  Alberta  Viscometer,  Fig.  (2) 

Component 

Number  Description 

1  Top  and  bottom  chucks  for  holding  torsion  wires. 

2  Torsion  wire  with  10.5"  from  top  of  chuck  to  top 
of  chuck. 


3  Elevator  assembly  with  screw  to  hold  bob  at  any 
depth  of  immersion. 

4  Dampening  cup  with  one  part  attached  to  the 
stator  and  the  other  anchored  to  the  elevator 
assembly. 

5  Compass  mounted  on  the  bob,  with  pointer  mounted 
on  one  leg  of  tripod. 

6  Bob  or  stator  with  thermocouple  resting  on  the 
bottom. 

7  Seal  between  the  cup  and  the  exterior  of  the  bath. 
The  lip  also  catches  excess  fluid  when  overflowing 
the  cup. 

8  Cup  or  rotor  with  hollow  shaft  for  filling  cup 
or  removing  liquid. 

9  Valve  for  controlling  liquid  flow  in  and  out  of 
cup . 

10  Pipe  connections  and  thermometer  at  inlet  and 
outlet  to  circulating  fluid  bath  surrounding  the 
cup. 

11  Double  bearing  assembly  to  hold  cup  rigidly  in 
place. 

12  Cast  iron  and  steel  2:1  reduction  gears. 

13  Revolution  counter,  gear  connected  to  cup  shaft. 

14  G-raham  variable  speed  (0  —  500  rpm.  )  drive  and 

synchronous  rotor. 

15  Channel  shaped  motor  support  to  reduce  the 
vibration  in  the  instrument. 
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removed  from  these  tins  hy  suction  ,  and  kept  in  enclosed 
containers  until  it  was  used  in  the  viscometers. 

A  Thomas  bath  was  used  to  keep  the  temperature 
constant  for  the  runs  from  50°F  to  80°F.  The  procedure 
used  for  the  runs  in  this  range  of  temperature  was  as  follows: 

With  the  correct  torsion  wire  mounted,  and  the 
circulation  bath  at  a  constant  temperature,  the  oil  was 
poured  in  from  the  top,  until  the  cup  overflowed.  A  zero 
reading  was  observed  and  the  temperature  of  the  bob  was 
determined  by  the  thermocouple  in  it.  It  was  found  that 
thermal  equilibrium  was  attained  more  quickly  by  rotating 
the  cup  at  0.5  to  0.9  rps.  When  the  thermocouple  in  the 
stator  indicated  thermal  equilibrium,  the  deflection  and 
rotary  speed  at  that  setting  were  observed  and  recorded. 

The  speed  was  then  reduced  in  steps  to  the  lowest  rotary 
speed,  where  the  bob  temperature  was  checked  for  thermal 
equilibrium.  At  each  step  the  speed  of  rotation  was 
determined  and  the  deflection  noted.  The  speed  was  then 
increased  by  the  following  increments  on  the  Graham 
variable  transmission;  1,2,5,9,12  or  15,  18  or  20,  22  or 
24.  The  last  two  or  three  settings  were  only  attained 
with  4"  immersion,  and  with  highly  viscous  oils.  The 
rise  in  temperature  at  the  high  speeds  was  sometimes 
detected  but  it  was  usually  so  small  that  it  could  not 
be  detected  by  the  thermocouple.  The  speed  was  reduced  in 
the  same  number  of  steps  to  the  lowest  speed.  If  there 
was  no  more  than  0#5°F  difference  between  the  initial  and 
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final  temperature,  the  run  was  considered  to  be  isothermal. 
For  the  range  of  temperatures  above  40°F  the  bath  was  kept 
at  a  constant  temperature  while  the  equipment  was  cleaned 
and  the  run  repeated  on  another  oil. 

For  the  runs  at  40°F  down  to  the  coldest  tempera¬ 
ture,  the  Precision  constant  temperature  bath  was  used.  This 
equipment  made  it  possible  to  cool  the  samples  while  the 
equipment  temperature  was  attaining  equilibrium.  Air  press¬ 
ure  was  used  to  force  the  pre-cooled  oil  from  the  sample 
bottles  into  the  instrument  through  the  bottom  of  the  cup. 
Less  time  was  required  for  thermal  equilibrium  with  this 
method  and,  after  a  test  at  one  temperature,  the  oil  was 
left  in  the  equipment  while  the  temperature  was  lowered. 

When  the  equipment  was  at  the  lower  temperature,  the  old  oil 
was  removed,  and  a  fresh  pre-cooled  sample  was  forced  into 
the  cup.  The  method  of  making  observations  was  similar  to 
the  other  runs  except  with  the  Pembina  crude,  as  it  tended 
to  set  up  in  the  cup,  and  give  high  apparent  viscosity 
values.  This  effect  was  reduced  by  starting  up  the  machine 
as  soon  as  possible  after  the  oil  was  introduced.  To  detect 
and  measure  the  thixotropy  of  the  Pembina  crude, runs  were 
made  so  that  the  total  times  required  to  go  from  0  to  a 
maximum  speed,  and  then  back  again,  were  equal, 

MacMichael  Rotational  Cylinder  Viscometer 

The  modified  MacMichael  used  in  these  experiments 
had  a  circulating  jacket  constructed  around  the  cup,  to 
enable  circulating  fluid  to  keep  the  jacket  at  a  constant 
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temperature.  It  has  been  further  modified  by  the  installa¬ 
tion  of  a  copper-constantan  thermocouple  in  the  stator  bob 
wall  to  measure  the  average  wall  temperature  at  a  depth  of 
immersion  of  3.5  cm.  This  made  it  possible  to  determine  the 
temperature  at  the  point  while  the  instrument  was  moving,  and 
also  to  detect  thermal  equilibrium. 

The  depth  of  immersion  was  kept  constant  at  3.5 
cm.  throughout  all  the  tests  conducted,  to  decrease  the 
variables,  as  it  was  quite  easy  to  change  the  torsion  wires 
when  over  ranged.  With  the  cup  full  of  oil,  the  circulating 
temperature  was  obtained  and  a  zero  deflection  reading  taken. 
The  cup  was  then  rotated  at  approximately  0.4  rps.  until 
thermal  equilibrium  was  obtained.  The  revolutions  were 
counted  by  eye  and  timed  with  a  stopwatch  for  thirty  to 
fifty  seconds,  and  the  net  deflection  was  obtained.  The 
full  speed  range  of  the  MacMichael  was  covered  in  four  to 
five  speeds  appropriately  spaced,  first  in  an  increasing, 
then  a  decreasing  speed,  to  detect  any  thixotropic  effects. 

If  any  change  in  temperature  had  occurred  during  any  run, 
the  run  was  repeated. 
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Fig.  3 

Vapor  and  liquid  equalizer  lines  at  upstream 
pressure  tap  of  pilot  pipeline.  Large 
pressure  gauge  indicates  bath  pressure. 


Fig.  4 

Pasco  gear  pump  at  left  bottom.  Roper  gear 
pump  at  center  right,  and  suction  and  dis¬ 
charge  manifold. 
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Fig.  5 

The  hy  pass  valves 
on  the  pipeline 
with  Roper  pump 
at  the  bottom, 
and  high  pressure 
lines  to  precool 
section. 


Fig.  6 

Pilot  pipeline  control 
board.  Upstream  and 
downstream  Heise  gauges 
and  low  pressure  down¬ 
stream  gauge.  Controls 
and  overload  switchboxes. 


Fig.  7 

Insulated  measuring 
tanks,  control 
panel,  and  Jerguson 
sight  glass  for 
liquid  Freon  12. 
Level  on  center 
right . 
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Fig.  8 

Alberta  Viscometer.  Note 
gear  connection  chain  for 
counter,  and  clips  for 
connecting  to  thermo¬ 
couple  wires  on  stator. 
G-raham  variable  speed 
drive  is  shown  on  right 
bottom. 


Overall  view  of  pilot  pipeline. 
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Fig.  10 

Alberta  Viscometer.  Cold  temperature 
operation.  Note  ice  on  circulating 
lines  and  on  viscometer  bath  jacket. 


Fig.  11 

From. left  to  right,  the  MacMichael,  the 

Qr+£lsloS  Scientific  Constant  Temoeratn-pc 
Bath,  and  the  Alberta  ViScomete??PwHh 
connecting  circulating  lines. 
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Results 

Tests  were  conducted  on  the  following  blends  and 
crudes,  using  the  pilot  pipeline,  the  MacMichael,  and  the 
Alberta  concentric  cylinder  rotational  viscometers. 

1.  Saskatchewan  Blend  consisting  of  39$  Fosterton 
34.5$  Success  and  26.5$  Cantaur  crudes  by  volume. 

2.  Combined  Blend  consisting  of  20$  Redwater 
and  80$  Saskatchewan  Blend  by  volume. 

3.  Fosterton  crude. 

4.  A  composite  sample  of  Pembina  crudes  for 
which  the  makeup  and  history  are  not  known. 

5.  Redwater  crude.  The  results  of  A.  Masuda  (19) 
on  the  tests  in  the  pilot  pipeline,  were  used.  (See  Table  14e). 

Pilot  Pipe  Line 

The  source,  API  gra.vity  and  the  dates  when  the 
crude  samples  were  received,  are  listed  in  Table  1-D.  The 
following  observations  were  made  while  conducting  the  tests 
on  the  pilot  pipeline  : 

1.  Upstream  pressure. 

2.  Downstream  pressure. 

3.  Volume  of  oil  collected,  and  time  required. 

4.  Temperature  of  oil  collected. 

5.  Upstream  inlet  temperature  to  the  test  line. 

6.  Downstream  outlet  temperature  from  the  test  line. 


...  J 

* 

*  '  ' 

’ 

' 

„ 

. 

* 

, 

. 


OH 


■ 

_  , 

•  ..  | 

* 


. 

. 

. 

. 


43 


Pilot  Pipe  Line 

This  data  is  recorded  in  laboratory  notebooks  and 
items  1,  2,  5,  and  6,  have  been  transcribed  to  Tables  1-E 
to  5-E  for  the  respective  crudes,  with  items  3  and  4  being 
combined  to  give  the  rate  of  flow  in  cubic  feet  per  minute. 

The  results  of  the  crude  gravity  determinations  as  a  function 
of  the  temperature  are  presented  in  Tables  2-D  and  in  Fig.(l-D). 

All  the  data  were  computed  as  explained  in  the 
Appendix  (E),  and  plotted.  Many  runs  have  been  omitted 
because  of  large  temperature  drops  from  inlet  to  outlet, 
and  changes  in  the  oil  inlet  or  outlet  temperature  during 
a  test.  In  some  tests,  the  temperature  increased  very 
quickly,  so  that  an  otherwise  good  run  would  be  in  error. 

Some  of  the  tests  on  the  Pembina  crude  had  unsteady  rates 
of  flow,  and  had  to  be  scrapped.  For  the  pilot  pipeline 
the  following  computed  results  are  tabulated  in  Appendix  E  : 

1.  Friction  factor. 

2.  Reynolds  number. 

3.  Apparent  viscosity. 

4.  Rate  of  shear  at  the  walls  of  the  pipe. 

5 .  Remarks . 

For  the  Saskatchewan  and  Combined  Blends,  and  the 
Fosterton  crude,  all  the  data,  was  obtained  in  the  laminar 
region  of  flow.  With  the  Pembina  crude,  it  was  necessary 
to  determine  both  the  laminar  and  turbulent  flow  character¬ 
istics  for  the  runs  with  Reynolds  numbers  between  1400  and 
2000,  as  determined  by  the  laminar  equation,  and  the  type  of 
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flow  is  stated  under  Remarks.  For  the  higher  viscosity 
oils  tested,  satisfactory  correlation  of  the  data  was 
obtained  by  plotting  the  apparent  viscosity  as  a  function 
of  the  temperature  for  the  range  of  shear  rates  investigated. 
With  the  Pembina  crude,  it  was  necessary  to  plot  the 
apparent  viscosity  as  a  function  of  the  shear  rate,  to 
obtain  a  better  means  of  correlation. 

To  determine  the  effect  of  pressure  on  viscosity, 
it  was  necessary  to  determine  the  viscosity  curve  at  a  low 
average  pressure,  and  then  compute  the  percentage  increase 
in  the  viscosity  for  the  total  difference  in  the  average 
pressures . 

3y  averaging  these  constants  over  limited  tempera¬ 
ture  ranges,  a  curve  may  be  plotted  to  show  the  decreasing 
effect  of  pressure  on  viscosity  with  increasing  temperature. 

MacMichael  and  Alberta  Rotational  Viscometer 

The  following  data  were  obtained  and  recorded  in 
notebooks  while  making  tests  with  the  Alberta  viscometer 
on  each  crude j 

1.  The  zero  reading  of  the  stabilized  torsion  wire 
deflection. 

2.  The  torsion  wire  gauge  and  wire  number. 

3.  The  immersion  depth  of  the  stator. 

4.  Bob  temperature. 

5.  Inlet  and  outlet  bath  temperature. 

6.  The  number  of  revolutions  the  cup  rotated  in  a 
measured  length  of  time. 
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For  the  MacMichael,  the  depth  of  immersion  was 
kept  constant  at  5.6  cm.,  and  only  one  temperature  of  the 
circulating  fluid  was  recorded.  Otherwise,  the  same 
observations  were  recorded  as  with  the  Alberta.  This  data 
was  used  to  compute  the  apparent  viscosity  and  respective 
temperature  as  outlined  by  the  sample  computations  in  Appendix 
G.  This  data  is  plotted  with  the  shear  rate  as  a  parameter 
to  detect  any  pseudoplastic  effects. 

With  the  Pembina  crude,  the  effect  of  the  shear 
rate  history  was  best  determined  with  the  rotational  visco¬ 
meters.  Because  of  this  property  it  was  possible  to  obtain 
a  range  of  viscosity  values.  The  apparent  viscosity  data 
plotted  at  one  temperature  was  the  lowest  value  obtained  for 
the  respective  shear  rate.  The  effect  of  shear  rate  history 
is  best  illustrated  by  Fig.  (28),  where  the  apparent  viscosity 
values  as  obtained  on  the  upcurve  and  downcurve  are  compared. 
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DISCUSSION  OF  THE  RESULTS 
Saskatchewan  Blend 

Pilot  Pipe  Line 

The  Saskatchewan  Blend  consisted  of  39$  Fosterton, 
34.5$  Success,  and  26.5$  Cantaur  crudes  by  volume,  with  a 
mixture  gravity  of  22.2°API.  In  Fig.  (13),  the  first  of 
the  two  curves  obtained  for  this  blend  is  defined  by  the 
tests  from  February  18  to  February  25,  1954.  On  February 
25,  an  accident  indicated  that  the  tanks  were  not  sealed, 
so  that  light  ends  were  lost  at  the  higher  temperatures. 

The  partially  weathered  crude  is  represented  by  the  tests 
from  February  26  to  March  6.  The  oil  used  in  the  rotational 
viscometer  tests  was  obtained  at  the  end  of  the  second  test. 

The  data  obtained  lie  within  10$  of  the 
smoothed  curve  for  one  temperature.  At  one  value  of  the 
apparent  viscosity,  the  temperature  varies  plus  or  minus 
2°F  from  the  smoothed  curve.  There  is  no  indication  that 
this  is  a  result  of  the  variation  in  the  shear  rate. 

All  the  tests  conducted  on  this  blend  were  in 
the  laminar  flow  region.  At  80°F  the  average  Reynolds 
number  was  600,  with  a  maximum  and  minimum  value  of  1,000 
and  200;  at  40°F  the  high  and  low  values  were  40  to  100; 
and  at  20°F,  20  to  40.  The  average  shear  rate  at  80°F  was 
8,100  sec."1;  at  40°F,  the  average  rate  of  shear  was  1,300 
sec.“-*-,  with  high  and  low  values  of  3,300  and  186  sec.”1; 
and  at  20°F,  the  average  was  640,  with  high  and  low  values 

of  1,000  and  380  sec.”1  respectively.  For  these  ranges  of 
shear  rates,  no  pseudoplastic  behavior  was  detected. 
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It  was  very  difficult  to  control  the  temperature 
of  this  oil  as  a  result  of  the  heat  generated  by  the  pump, 
throttling  valve,  and  the  friction  loss  in  the  fittings. 

The  high  viscosity  at  low  temperatures  seriously  reduced  the 
heat  transfer,  so  that  large  temperature  gradients  and  low 
flow  rates  were  necessary  to  attain  the  lower  temperatures. 
Because  of  these  conditions,  it  was  necessary  to  omit  many 
of  the  tests  as  a  result  of  a  large  temperature  drop  from 
inlet  to  outlet  and/or  the  temperature  being  erratic  and 
unstable. 

The  equation  for  non-isothermal  flow  as  advocated 
by  Keevil  and  McAdam  (15)  may  be  applied  to  this  crude. 

The  change  in  the  friction  factor  as  determined  by  this 
equation  at  40°F,  is  approximately  a  30/  increase,  so  the 
apparent  viscosity  values  obtained  may  be  30/  high;  at  30°F 
the  apparent  viscosity  values  may  be  45/  higher  than  for 
isothermal  flow. 

MacMichael  Viscometer 

Above  45°F,  the  data  obtained  with  this  instrument 
results  in  a  smooth  curve,  Fig. (13),  with  an  indication  of 
weathering  at  79°F.  Below  this  temperature,  the  data  are 
scattered,  possibly  as  a  result  of  the  temperature  gradients 
and  the  variables  such  as  room  temperature.  In  obtaining 
the  average  oil  temperature,  there  are  no  corrections  for 
variations  in  the  room  temperature,  and  the  higher  thermal 
gradients  across  the  annulus  at  low  circulating  liquid 
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This  oil  exhibits  pseudoplast ic  behavior  over  the 
full  temperature  range  in  the  MacMichael  viscometer,  with 
the  minimum  pseudoplasticity  being  exhibited  from  50°F  to 
80°F,  and  the  maximum  effect  exhibited  at  20°F. 

Alberta  Viscometer 

The  apparent  viscosity  data  obtained  with  the 
Alberta  viscometer  is  plotted  on  Fig,  (14),  and  good 
correlation  is  obtained  from  80°F  to  30°F,  below  which  the 
data  vary  as  much  as  20$  in  the  apparent  viscosity,  or 
2°F  in  the  temperature.  This  may  be  a  result  of  either 
thermal  instability  or  an  instrumental  error,  such  as  high 
heat  conduct ion^and  poor  circulation  around  the  thermometers. 

The  oil  exhibits  pseudoplasticity  below  40°F  with 
the  effect  of  shear  rate  on  the  apparent  viscosity,  getting 
larger  with  decreasing  temperatures.  Below  40°F,  it  was  not 
possible  to  attain  a  shear  rate  high  enough  that  a  constant 
value  of  the  apparent  viscosity  was  obtained. 

This  crude  did  not  exhibit  any  weathering  action 
at  the  higher  temperatures,  possibly  because  it  was  allowed 
to  weather  in  the  pilot  pipeline  tanks. 

Correlation  of  Data 

The  smoothed  data  from  the  three  methods  is 
presented  in  Fig.  (15).  The  data  of  the  Alberta  and 
MacMichael  viscometers  correlate  within  ten  percent  from 
80°F  down  to  40°F,  at  which  temperature  the  MacMichael 
exhibits  values  higher  than  or  approximately  a  temperature 
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difference  of  1°F  warmer  than  the  Alberta  data.  This  may 
be  a  result  of  the  method  used  to  measure  the  oil  tempera¬ 
ture  on  the  MacMichael.  Because  the  data  of  the  Alberta 
viscometer  is  at  a  higher  rate  of  shear  and  is  more  consist¬ 
ent,  this  curve  will  be  compared  with  the  pilot  pipeline 
data. 

At  80°F,  the  pilot  pipeline  data  is  3 0 %  higher 

than,  and  at  20°F,  this  data  is  50^  lower  than,  the  Alberta 

data.  Below  40°F  this  discrepancy  may  be  a  result  of  the 

pseudoplastic  properties  of  the  blend.  This  effect  does 

increase  with  lower  temperatures.  The  maximum  shear  rate 

-1 

in  the  Alberta  is  50  to  PO  sec.  ,  whereas  the  pipeline 
shear  rates  at  the  wall  vary  from  640  to  1300  sec.“\ 

Above  40°F,  the  increasing  discrepancy  of  the 
data  may  be  a  result  of  the  large  temperature  gradient  from 
the  pipe  to  flowing  liquid.  Thi3  large  temperature  gradient 
is  attributed  to  the  low  heat  exchange  as  a  result  of  high 
viscosities  and  laminar  flow.  By  using  the  Keevil  and 
McAdams  friction  factor  equation  (15)  for  non- isothermal 
cooling  of  a  liquid,  it  is  estimated  that  the  apparent 
viscosity  data  is  30%  high  at  80°F,  as  compared  to  the  actual 
difference  of  25%.  The  increase  in  the  apparent  viscosity 
as  a  result  of  pressure  is  estimated  to  be  7%, (Fig.  35) > 
although  this  value  varies  with  the  test  pressure  and  tempera 
ture. 

Since  the  pilot  pipeline  was  operated  under 
non*»i  so  thermal  conditions  as  a  necessity,  and  because  of 
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the  correlation  ootained  with  the  two  rotational  visco- 
meters,  it  is  believed  that  the  Alberta  data  from  40°F 
to  80°F  represents  the  best  determination  of  the  apparent 
viscosity.  Below  40°F,  the  rotational  viscometer  data  is 
high,  possibly  as  a  result  of  the  apparent  viscosity  being 
some  function  of  the  shear  rate. 

A  composite  curve  consisting  of  the  Alberta  data 
above  40°F,  and  the  pilot  pipeline  data  below  this  tempera¬ 
ture  is  presented  in  Fig.  (34),  as  a  compilation  of  the 
best  data. 


COMBINED  BLEND 

Pilot  Pipe  Line 

The  Combined  blend  consists  of  80%  by  volume  of 
the  weathered  Saskatchewan  blend  plus  20%  by  volume  Redwater 
crude.  The  apparent  viscosity  data  was  obtained  in  laminar 
flow  only,  with  the  Reynolds  numbers  varying  from  an  average 
at  70°F  of  1,200  to  less  than  1.0  at  0°F.  The  average 
Reynolds  numbers  and  shear  rates  at  the  wall  at  various 
temperatures  are:  70°F,  1200  and  11,000  sec,”^;  at  40°F 
96  and  2,300  sec.-1;  and  at  20°F,  75  and  1,400  sec/1, 

It  was  possible  to  obtain  more  consistent  data 
with  this  crude  because  of  better  temperature  control,  and 
the  lower  viscosity.  The  curve  was  drawn  through  the  lower 
body  of  the  data  in  Fig.  (16)  as  these  points  had  the  small¬ 
est  changes  in  temperature  from  inlet  to  outlet,  and  the 
best  thermal  equilibrium.  At  one  temperature,  the  data  lie 
between  plus  15%  and  minus  10%  of  the  smooth  curve  value. 
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TEMPERATURE 


If  the  apparent  viscosity  values  are  assumed  to  be  correct, 
the  temperature  correction  is  plus  3°F  or  minus  1°F. 

The  increase  in  the  apparent  viscosity  as  a 
result  of  the  non* isothermal  flow  at  70°F  to  80°F  is 
estimated  to  be  35%  to  40%.  Below  10°F,  it  is  expected 
that  the  oil  temperature  approaches  the  temperature  of  the 
pipe  wall  at  the  low  flow  rates  used. 

There  is  no  indication  from  the  pilot  pipeline 
data  that  the  rate  of  shear,  or  pipe  size,  has  any  effect 
on  the  apparent  viscosity  of  this  crude. 

MacMichael  Viscometer 

The  apparent  viscosity  data  obtained  on  the 
MacMichael  viscometer,  indicates  that  the  blend  is  pseudo* 

Q 

plastic  below  70  F,  as  shown  in  Fig.  (17).  With  the 
Saskatchewan  crude  oils  tested,  pseudoplastic  behavior 
was  not  evident  above  a  certain  temperature,  and  below 
this  temperature  the  pseudoplasticity  increases  with 
lower  temperatures.  The  constant  variation  exhibited  on 
the  MacMichael  viscometer  in  the  apparent  viscosity  may 
be  a  result  of  an  instrumental  error,  such  as  the  bob 
oscillating  about  the  center  of  the  cup,  or  variation  of 
the  room  temperature.  The  erratic  data  below  30°F  is 
believed  to  be  a  result  of  incorrect  oil  temperatures.  The 
blend  exhibited  some  weathering  at  the  higher  temperatures, 
as  the  apparent  viscosity  increased  with  time. 
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Alberta  Viscometer 

The  apparent  viscosity  data  obtained  are 
plotted  on  Fig.  (18),  with  the  indication  that  this 
blend  is  pseudoplastic  below  55°F.  For  an  increase  in  the 
shear  rate  from  3.3  to  90/sec.,  at  20°F,  there  is  a  100% 
increase  in  the  apparent  viscosity.  On  the  three  low 
temperature  runs,  some  thixotropic  behavior  was  indicated, 
as  the  same  apparent  viscosity  was  not  obtained  for  the 
same  shear  rate,  and  in  all  cases  the  apparent  viscosity 
was  lower  when  the  previous  shear  rate  history  was  higher. 

The  discrepancy  of  the  data  at  30°F  is  believed 
to  be  a  result  of  using  the  same  oil  temperature  correction 
for  three  and  four  inch  depths  of  immersion  at  the  colder 
temperatures.  Some  erratic  data  was  obtained  at  62°F, 
but  this  is  assumed  to  be  accidental  error. 

Weathering  of  the  crude  occurred  at  the  higher 
temperatures  as  the  apparent  viscosity  Increased  with  time, 

Correlat ion 

The  smoothed  curve  of  the  data  of  the  MacMichael 
and  Alberta  viscometers  resulted  in  a.  good  correlation  from 
50°F  to  70°F,  and  at  80°F  the  crude  in  the  MacMichael  appears 
to  have  been  weathered.  The  divergence  of  the  data  below 
45°F  in  Fig.  (19)  is  believed  to  be  a  result  of  pseudo¬ 
plasticity  with  the  15  sec,-"  curve  of  the  MacMichael 
viscometer  being  equivalent  to  a  rate  of  shear  of  5  to  8 
sec."^-  on  the  Alberta  viscometer.  The  variance  may  have 


been  caused  by  incorrect  oil  temperature  measurement  on  the 
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MacMichael.  As  the  Alberta  curve  represents  more  consistent 
data,  it  will  be  used  for  correlation. 

The  smoothed  curve  representing  the  pilot  pipeline 
data  is  38%  higher  than  the  Alberta  viscometer  curve  at  75°F 
and  is  50%  lower  at  20°F.  As  this  crude  has  exhibited 
pseudoplastic  properties  at  low  shear  rates,  the  smoothed 
data  taken  at  high  shear  rates  in  the  pilot  pipeline  should 
be  lower  than  the  Alberta  data.  If  the  estimated  correction 
of  35%**4-0%  for  non- isothermal  flow  were  applied  to  the 
pilot  pipeline  data  at  75°F  to  80°F,  these  two  curves  would 
correlate  within  5%  above  50°F.  The  effect  of  pressure  on 
the  apparent  viscosity  may  have  increased  the  data  of  the 
pilot  pipeline  7%  (See  Fig.  35). 

By  combining  the  Alberta  Viscometer  smoothed  data 
above  50°F,  and  using  the  pilot  pipeline  curve  below  50°F, 
the  resulting  curve  may  represent  the  apparent  viscosity 
of  the  Combined  blend  at  a  high  rate  of  shear.  This  comb  inac¬ 
tion  of  data  exhibits  a  straight  line  on  the  Saybolt  tempera¬ 
ture  paper  on  Fig.  (3^),  which  is  an  indication  of  a 
Newtonian  liquid.  The  temperature  at  which  the  composite 
curve  exhibits  pseudoplastic  properties  on  this  plot  is 
30°F,  for  the  two  straight  lines  intersect  at  this  tempera¬ 
ture  on  the  Saybolt  viscosity  temperature  paper.  The 
effect  of  pressure  in  the  pilot  pipeline  may  have  increased 
the  value  of  the  apparent  viscosity  as  much  as  9%  at  20°F 
and  6%  at  80°F,  compared  to  the  value  which  would  be  obtained 
at  atmospheric  pressure  (Fig,  35). 
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Combined  Blend 

Effect  of  Pressure  on  the  Apparent  Viscosity 

The  high  pressure  runs  have  been  used  to  compute 
the  pressure  effect  constant  HC"  by  determining  the  percent¬ 
age  increase  in  the  apparent  viscosity  per  100  psi.  increase 
in  pressure.  These  constants  have  been  averaged  for  temper¬ 
ature  ranges  of  10°F.  A  plot  of  this  data  is  presented  in 
Fig.  (20)  where  the  apparent  viscosity  curve  of  the  pilot 
pipeline  at  655  psi.  is  shown  as  compared  to  the  lower 
curve  representing  the  values  at  105  psi.  The  effect  of 
the  temperature  on  the  pressure  effect  as  determined  with 
this  crude,  is  presented  in  Fig.  (35),  with  the  highest 
values  of  C  occurring  at  the  lowest  temperatures. 

FOSTERTON 

Pilot  Pipeline 

The  apparent  viscosity  data  is  plotted  in  Fig,  (21) 
from  2.5°F  to  75°F.  All  the  tests  on  the  Fosterton  crude 
occuri^ed  in  the  laminar  region  of  flow,  and  at  various 
temperatures  the  average  values  of  the  Reynolds  number  and 
shear  rate  at  the  wall,  are  as  follows:  at  ?0°F,  523  and 
4,000  sec."1;  at  40°F,  73  and  1,730  sec.~x;  at  20°F,  33  and 
1,700  sec.”1;  at  3°F,  the  Reynolds  number  varies  from  0.5 
to  9.0  and  the  average  rate  of  shear  is  160  sec.  .  For 
the  range  of  shear  rates  obtained,  there  is  no  indication 
that  the  shear  rate  or  pipe  size  has  any  effect  on  the 
apparent  viscosity  of  the  Foster-ton  crude. 
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The  estimated  increase  in  the  apparent  viscosity 
as  a  result  of  non-isothermal  flow  is  estimated  to  be  35 % 
to  bQ%  at  70°F  to  75°F,  and  to  20%  at  40°F.  The  devia¬ 
tion  of  the  data  from  smoothed  curve  of  Fig.  v2l)  is 
approximately  plus  or  minus  5 %>. 

The  pilot  pipeline  data  for  this  crude  forms  a 
straight  line  down  to  20°F  on  the  ASTM  Saybolt  viscosity 
temperature  paper,  indicating  that  this  crude  oil  may  be 
Newtonian  above  20°F. 

Fosterton  Crude 
MacMichael 

The  data  obtained  is  plotted  in  Fig.  (22)  and 
the  best  fit  smooth  curve  has  been  constructed  with  the 
discrepancy  at  80°F  assumed  to  be  a  result  of  weathering. 
The  data  indicates  that  this  crude  is  slightly  pseudo¬ 
plastic  over  the  full  range  of  temperature.  The  constant 
variation  in  the  apparent  viscosity  may  be  a  result  of 
instrumental  errors,  such  as  the  bob  going  off  center  at 
the  low  rates  of  shear,  as  the  runs  as  60°F  show  no  effect 

on  the  apparent  viscosity  by  the  shear  rate.  A  small 

o 

amount  of  weathering  occurred  w ith  time  above  55  F.  The 
change  in  curvature  at  40°F  may  indicate  the  transition 
temperature  from  Newtonian  to  pseudoplastic  behavior,  as 
determined  at  these  low  rates  of  shear. 

Fosterton 

Alberta  Viscometer 

The  apparent  viscosity  data  obtained  with  the 
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Alberta  viscometer  results  in  very  little  scattering  from 
the  smoothed  curve,  except  at  80°F,  where  the  apparent 
viscosity  increased  rapidly  with  time,  as  a  result  of 
weathering. 

The  effect  of  weathering  was  apparent  down  to 
50°F,  with  the  result  that  any  pseudoplastic  effects  may 
have  been  hidden  by  the  increasing  values  of  the  apparent 
viscosity . 

The  variation  in  the  shear  rate  below  40°F 
results  in  a  small  change  in  the  apparent  viscosity.  On 
the  runs  at  32°F  and  24°F,  there  are  indications  that  the 
limiting  value  of  the  apparent  viscosity  may  be  obtained 
by  a  small  increase  in  the  shear  rate  above  78.0  sec.”\ 

The  data  obtained  indicate  that  the  Fosterton  crude  is 
slightly  pseudoplastic,  with  the  effect  being  most  apparent 
below  45°F. 

Fosterton  Crude 
Correlation 

The  data  of  the  Alberta  and  MacMichael  viscometers 
correlate  within  five  percent  above  40°F,  and  below  this 
temperature  the  correlation  is  poorer,  possibly  as  a  result 
of  the  larger  shear  rates  possible  in  the  Alberta  viscometer. 
At  80°F  the  data  obtained  by  both  methods  is  assumed  to  be 
too  high,  as  a  result  of  weathering.  The  data  obtained 
above  60°F  may  be  higher  as  a  result  of  weathering. 

The  correlation  of  the  Alberta  and  pilot  pipeline 
data  on  the  Fosterton  crude  may  be  better  because  the  oil 
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is  less  viscous  and  less  pseudoplastic  than  the  other 
blends.  The  good  correlation  between  the  Alberta  and 
pilot  pipeline  data  at  25°F  may  be  an  indication  that  the 
shear  rate  in  the  Alberta  viscometer  is  high  enough  to 
obtain  a  constant  value  of  the  apparent  viscosity.  It 
is  estimated  that  the  apparent  viscosity  data  of  the 
pilot  pipeline  is  35%  to  40%  high  at  70°F  as  a  result  of 
the  oil  to  pipe  temperature  gradient,  so  that  the  curve 
of  the  Alberta  data  would  be  5%  lower  than  the  correlated 
data.  This  small  discrepancy  of  5%  may  be  a  result  of  the 
pressure  effect  on  the  apparent  viscosity. 

By  using  the  smoothed  data  for  the  Alberta 
viscometer  above  25°F,  and  the  pilot  pipeline  data  below 
25°F,  the  resulting  composite  curve  is  a  linear  relation¬ 
ship  on  the  ASTM  Saybolt  viscosity  temperature  graph  paper, 
on  Fig.  (34),  which  may  indicate  that  a  shear  rate  above 
80  sec.”'1'  has  no  effect  on  the  apparent  viscosity  at  30°F 
or  higher,  for  this  crude. 

Effect  of  Pressure  on  the  Apparent  Viscosity 
of  the  Fosterton  Crude 

The  apparent  viscosity  data  on  the  pilot  pipeline 
at  pressures  varying  from  70  to  above  800  psi.  are 
presented  in  Fig.  (25).  As  a  result  of  the  scattering  in 
the  data,  it  was  necessary  to  average  the  pressure  effect 
coefficients  "C"  listed  in  Table  5-E,  for  10°F  ranges  in 
temperature.  These  values  were  used  to  determine  the 
position  of  the  high  pressure  curve  in  Fig.  (25),  and  the 
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plot  of  the  pressure  effect  coefficient  as  a  function  of 
the  temperature  in  Fig.  (35).  The  pressure  effect  coefficient 
at  70°F  is  3$  t  1%  and  at  22°F  it  is  10%  t  1%  per  100  psi. 

The  variation  in  the  temperature  for  the  data  may  he 
approximately  -  2°F,  so  that  any  correction  in  the  downstream 
temperature  reading  would  he  within  the  experimental  error 
for  most  of  the  high  pressure  runs. 

PEMBINA  CRUDE 

Pilot  Pipe  Line 

The  data  obtained  on  the  pilot  pipeline  indicates 
that  the  apparent  viscosity  of  the  Pemhina  crude  is  a  function 
of  the  shear  rate  helow  10,000  sec."1  from  10°F  to  80°F,  and 
the  change  in  the  apparent  viscosity  above  10,000  sec.”1 
is  very  small.  The  data  was  best  correlated  as  a  function 
of  the  shear  rate  at  the  wall  with  temperature  as  a 
parameter,  as  in  Fig.  (26).  There  were  275  runs  made  on  the 
Pembina  crude,  and  this  data  was  improved  by  omitting  the 
test  runs  with  large  temperature  drops,  variable  tempera¬ 
tures  during  a  run,  and  an  unsteady  rate  of  flow. 

To  determine  the  position  of  the  isothermals 

from  the  data,  it  is  necessary  to  average  the  data  for 

various  logarithmic  increments  in  the  shear  rate  for  each 
o 

10  F  range  in  temperature.  The  average  deviation  of  the 
data  is  "i"  20%  from  the  smoothed  curve,  with  some  of  the 
data  varying  as  much  as  100$. 
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The  data  indicates  that  the  apparent  viscosity  is 
a  function  of  the  shear  rate  history  and  the  test  pipeline 
diameter,  so  that  the  Pembina  crude  may  be  highly  thixotropi 
In  the  pilot  pipeline  tests,  the  shear  rate  history  was 
dependent  on  the  pump  used,  the  amount  of  liquid  by-massed 
at  the  pumps  or  pre-cool  section,  and  the  shear  rate  history 
of  the  oil  in  the  supply  tank. 

The  effect  of  the  shear  rate  on  the  apparent 
viscosity  increases  below  70°F  to  a  maximum  "at  45°F.  The 
slope  of  the  isothermal  curve  for  50°F  from  100  sec.“^  to 
10,000  sec/^  in  Fig.  (26)  is  approximately  -0.3,  so  that 
the  apparent  viscosity  is  approximately  inversely  proportion 
al  to  the  cube  root  of  the  shear  rate.  The  slope  of  the 
isothermal  curve  for  20°F  over  the  same  limits  is  -0.5,  so 
that  the  apparent  viscosity  is  approximately  inversely 
proportional  to  the  square  root  of  the  shear  rate. 

The  smoothed  data  of  Fig.  (26)  was  used  to  plot 

the  apparent  viscosity  as  a  function  of  the  temperature 

with  the  shear  rate  as  a  parameter  in  Fig.  (30).  On  this 

-1 

plot  there  is  only  data  for  the  80,000  sec.  shear  rate 
curve  down  to  40^F,  because  it  was  not  possible  to  obtain 
turbulent  flow  below  this  temperature.  The  apparent 
viscosity  data  for  the  100  to  200  sec.“^  rate  of  shear 
curves  is  erratic  because  it  "was  not  possible  to  maintain 
a  constant  rate  of  flow  during  a  run  at  the  low  rates. 

Also,  it  was  very  difficult  to  obtain  reliable  data  for 
the  100  and  200  sec.”1  rate  of  shear  curves  above  50°F  and 
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60°F  respectively,  because  of  the  very  low  pressure  drops 
necessary  for  the  low  rates  of  flow. 

The  apparent  viscosity  curve  for  a  shear  rate 
of  80,000  sec.”1  results  in  two  straight  lines  intersecting 
at  75°F  on  the  ASTM  Saybolt  viscosity  temperature  graph 
in  Fig.  (34).  This  indicates  that  the  oil  may  change  from 
a  Newtonian  to  a  non-Newtonian  liquid  at  75°F.  The  smoothed 
curve  for  a  rate  of  shear  of  3,000  sec.-1,  results  in 
three  straight  lines  intersecting  at  70°F  and  30°F.  The 
reason  for  the  two  points  of  intersection  is  not  known,  and 
it  can  only  be  propounded  that  there  may  be  two  changes  in 
the  structure  of  the  Pembina  crude,  or  the  effect  of  a 
temperature  gradient  in  laminar  flow  may  vary  over  the 
temperature  range  covered. 

Pembina  Crude., Alberta  Viscometer 

In  Fig.  (27),  the  apparent  viscosity  data  obtained 
at  four  different  rates  of  shear  at  each  temperature,  were 
used  to  construct  four  smooth  curves  for  the  range  of 
temperatures  covered.  From  this  figure,  it  may  be  seen 
that  the  pseudoplasticity  increases  from  a  minimum  at  80°F 
to  a  maximum  at  45°F,  and  then  decreases  to  a  minimum  value 
at  30°F  and  lower.  The  degree  of  thixotropy  exhibited  was 
found  to  be  a  maximum  at  45°F,  as  determined  by  the  Alberta 
viscometer  data. 

Because  the  apparent  viscosity  is  so  dependent 
on  the  shear  rate  history,  as  illustrated  by  a  comparison 
of  the  runs  at  50°F  and  60°F,  each  sample  was  subjected  to 
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a  rate  of  shear  of  approximately  64  sec.-*1,  for  a  sufficient 
period  of  time  to  reach  equilibrium  prior  to  the  determination 
of  the  apparent  viscosity  at  the  lower  rates  of  shear.  At 
the  higher  temperatures,  this  was  not  possible  because  of 
the  high  degree  of  weathering.  The  effect  of  time  at  one 
rate  of  shear  is  illustrated  in  Fig.  (28),  where  the 
upcurve ,  plateau  at  64.5  sec.  1,  and  the  downcurve  form  a 
type  of  hysteresis  curve.  This  is  also  illustrated  in 
Fig.  (26)  at  25.2°F.  Minimum  values  for  the  apparent 
viscosity  may  also  be  obtained  by  subjecting  the  crude 
oil  to  a  higher  shear  rate,  and  then  determining  the 
apparent  viscosity  at  decreased  values  of  the  shear  rate. 

An  example  of  this  phenomena  is  illustrated  in  Fig. (28) 
by  the  42.8°F  curve. 

The  reduced  degree  of  thixotropy  below  45°F 
may  be  the  result  of  two  instrumental  errors.  First, 
since  this  crude  has  a  high  apparent  viscosity  at  low 
rates  of  shear,  the  shearing  stress  across  the  annulus 
may  not  be  high  enough  to  shear  all  the  oil  in  the 
annulus.  This  would  result  in  a  high  value  of  the 
apparent  viscosity  with  very  little  change  in  the 
apparent  viscosity  with  shear  rate.  Second,  with  the 
low  rates  of  shear  on  the  end  of  the  bob,  the  apparent 
viscosity  as  a  result  would  be  very  high,  so  that  the 
end  effect  may  be  increased.  This  would  result  in  tne 
apparent  viscosity  values  at  low  rates  of  shear  being 
more  consistent,  as  the  shear  rates  in  the  annulus 
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would  be  approaching  the  shear  rate  at  the  end  of  the  bob. 
This  may  be  the  reason  for  the  60  sec.  1  shear  rate  curve 
on  Fig.  (27),  not  exhibiting  similar  curvatures  and 
inflection  points  as  the  3.7  sec.-1  curve. 

MacMichael  Viscometer 

The  apparent  viscosity  data  is  plotted  as  a 
function  of  the  temperature  in  Fig.  (29),  with  the  shear 
rate  as  a  parameter.  This  instrument  detected  a  wide 
spread  in  the  apparent  viscosity  from  80°F  to  20°F  with 
a  small  change  in  the  shear  rate.  The  scattering  of  the 
higher  shear  rate  points  at  60°F,  50°F,  and  44°F,  is 
indicated  to  be  a  result  of  the  variation  in  the  shear 
rate  history.  The  effect  of  the  shear  rate  on  the  apparent 
viscosity  increases  from  80°F  to  53°F,  and  then  decreases 
to  a  constant  variation  at  30°F  and  lower. 

The  limitations  of  the  MacMichael  viscometer 
for  a  crude  of  such  pseudoplasticity  and  thixotropy  are 
the  small  range  in  the  shear  rate,  and  the  large  annular 
clearance  between  the  bob  and  cup. 

Correlation  of  the  Pembina  Crude 

A  comparison  of  the  MacMichael  and  Alberta 
viscometer  data  is  presented  in  Fig.  (29).  It  may  be 
seen  that  the  apparent  viscosity  curve  for  15  see.-1 
in  the  MacMichael  viscometer  is  roughly  equivalent  to 
the  3.7  sec."1  curve  for  the  Alberta.  This  poor  correla¬ 
tion  may  be  a  result  of  the  broad  clearances  between  the 


' 


v  "i 

. 


;■  .  l-  ...  ■  :  -  - 

. 

i  v.' 


. ,  .  • 


.  .  ... 


.  I 


r 

■ 

; 

• 

!  .. 

....  -  .  - 

■ 

- 

• 

.  ,  i.  ~  i  V 

' 

<  ■ 

. 


,, 

. 


- 

, 

.  V  • 

• 

APPARENT  VISCOSITY  X  IOO  LB./FT.SEC. 


go 


81 


MacMichael  cup  and  bob,  with  the  result  that  the  shear 
stress  in  the  annulus  was  not  large  enough  to  shear  all 
the  oil  present.  Because  of  the  higher  rates  of  shear, 
and  more  consistent  shear  rate  history,  and  because  of 
the  better  temperature  measurement,  the  smoothed  curves 
of  the  Alberta  viscometer  data  will  be  used  for  correla¬ 
tion. 

The  data  obtained  by  the  Alberta  viscometer  and 
pilot  pipeline  are  plotted  on  log-log  coordinates  in  Fig. 
(26).  It  may  be  seen  that  the  Alberta  data  obtained  at 
70°F,  60°F,  50°F,  and  40°F,  correlates  approximately  with 
the  data  obtained  on  the  pilot  pipeline,  so  that  the  range 
of  shear  rates  for  these  isothermal  curves  varies  from  3.7 
to  10,000  sec.“^.  The  discrepancy  in  the  Alberta  data  at 
35°F  and  25°F  is  believed  to  be  a  result  of  the  instrumental 
errors  as  above. 

On  the  Alberta  viscometer,  the  shear  rate  can 
be  more  precisely  determined,  and  the  shear  rate  history 
better  controlled,  than  on  the  pilot  pipeline.  The  shear 
rate  history  in  the  pilot  pipeline  is  highly  variable  as  a 
result  of  passing  the  liquid  through  one  of  two  gear  pumps, 
one  of  two  throttling  valves,  and  by  the  passage  through 
sections  of  high  turbulence  prior  to  entering  the  test 
lines.  As  a  result  of  the  high  turbulence  in  the  equipment 
prior  to  entering  the  test  pipeline,  the  rate  of  shear 
history  is  high,  so  that  the  measured  apparent  viscosity 
of  the  oil  is  a  minimum.  In  the  Alberta  viscometer,  the 
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shear  rate  history  is  lower  than  the  highest-  rate  of  shear 
attained,  so  that  the  apparent  viscosity  at  the  highest 
possible  rate  of  shear  will  be  higher.  The  two  methods 
may  be  determining  the  values  of  the  apparent  viscosity 
above  and  below  the  true  value  for  the  100  sec.”1  rate  of 
shear. 

In  Fig.  (30),  the  apparent  viscosity  data  is 

plotted  as  a  function  of  the  temperature  with  the  shear 

rate  as  a  parameter.  Good  correlation  is  indicated  by  the 

smoothed  data  of  the  Alberta  viscometer  and  the  pilot 

-1 

pipeline,  except  for  the  60  sec.  curve  below  50°F. 

The  general  curvature  of  the  3.7  sec.”1’  curve  is  very 
similar  to  the  curves  for  the  pilot  pipeline.  The  Pembina 
crude  tested  in  the  rotational  viscometers  may  have  been 
more  weathered  as  a  result  of  the  high  pressure  tests, 
more  tests,  and  the  long  storage  period  at  high  tempera¬ 
tures,  with  the  periodic  loss  of  light  ends  when  obtaining 
sample  s . 

Effect  of  Pressure  on  the  Apparent  Viscosity 
of  the  Pembina  Crude 

To  determine  the  pressure  effect  coefficient  for 
the  apparent  viscosity  of  the  Pembina  crude,  it  was  necessary 
to  obtain  tests  at  the  same  shear  rate  at  the  one  temperature 
for  different  values  of  the  average  pressure  in  the  section. 
As  there  is  no  quick  method  of  determining  or  setting  the 
rate  of  shear  for  any  one  test  on  the  pilot  pipeline,  trial 
and  error  methods  were  used.  The  data  obtained  were  good 
in  only  a  few  of  the  many  tests.  The  average  value  of  the 
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pressure  effect  coefficient  was  14^/100  psi.,  and  the  data 
is  too  limited  to  indicate  the  effect  of  temperature  on  this 
coefficient.  This  value  is  shown  as  a  horizontal  line  in 
Fig.  (32) ,  over  the  temperature  range  covered. 

RED WATER  CRUDE 

Pilot  Pipe  Line 

The  apparent  viscosity  data  for  two  samples  of 
Redwater  crude  (Table  14-E) ,  were  obtained  by  A.  Masuda  (19) 
on  the  same  pilot  pipeline  as  used  in  this  investigation. 

The  apparent  viscosity  data  for  these  two  samples  has  been 
averaged,  and  the  resulting  values  have  been  used  to  construct 
the  curve  in  Fig.  (33).  There  is  no  indication  from  the 
original  data  of  pseudoplastic  behavior  in  the  Redwater 
crude,  so  that  one  value  of  the  apparent  viscosity  is 
sufficient  for  one  temperature  in  the  laminar  or  turbulent 
region.  Most  of  this  data  was  obtained  at  shear  rates 
above  1,000  see."-*-,  as  only  two  runs  were  in  the  region  of 
700  sec."1.  As  most  of  the  runs  above  60°F  were  in  the 
turbulent  region  of  flow,  the  effect  of  a  temperature 
gradient  on  the  apparent  viscosity  would  be  much  smaller 
than  in  laminar  flow.  The  data  obtained  results  in  two 
straight  lines  on  the  AST M  Saybolt  viscosity  temperature 
graph.  Fig.  (34),  indicating  that  the  transition  tempera¬ 
ture  from  a  Newtonian  to  a  non-Newtonian  liquid  may  be  at 
55°F. 
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MacMichael  Viscometer 

The  apoarent  viscosity  data  is  plotted  in 
Fig.  (31),  and  the  best  fit  curves  for  the  shear  rates  of 
15  sec.  ^  and  2.5  sec.  ^  are  illustrated.  The  data  above 
40°F  indicates  pseudoplastic  behavior,  but  the  variation 
in  the  apparent  viscosity  may  be  a  result  of  viscometer 
defects  and  also  weathering.  Below  35°F,  the  pseudoplastic 
behavior  is  more  evident,  but  poor  correlation  is  obtained, 
possibly  as  a  result  of  incorrect  oil  temperature. 

The  MacMichael  viscometer  is  seriously  under¬ 
damped  for  the  Redwater  crude  above  50°F.  There  is  no 
means  on  the  modified  MacMichael  to  dampen  the  bob  oscilla¬ 
tions.  Another  instrumental  error  detected  while  testing 
the  Redwater  crude,  was  the  uncontrollable  position  of  the 
top  of  the  bob.  With  a  B.S.  and  G-.  #  34  torsion  wire,  there 
is  a  clearance  of  0.028",  so  that  the  bob  is  free  to  move 
sideways  this  amount.  This  amount  of  movement,  plus  the 
slight  bending  caused  by  the  weight  of  the  stator,  leaning 
against  the  wire,  results  in  a  substantial  change  in  the 
deflection  as  illustrated  by  the  tests  at  78°F  and  63°F 
(Table  5-C). 

The  effect  of  weathering  was  very  serious  as  a 
result  of  the  time  necessary  to  attain  bob  equilibrium, 
and  as  a  result  of  the  quick  increase  in  the  apparent 
viscosity  with  time. 

Alberta  Viscometer 

The  apparent  viscosity  data  of  the  Alberta 
viscometer  is  plotted  in  Fig.  (32),  and  the  smoothed 
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curves  indicate  that  this  crude  is  a  Newtonian  liquid  from 
80°F  to  40°F,  and  at  the  lower  temperatures  small  pseudo¬ 
plastic  effects,  are  apparent*.  The  variation  in  the  value 
of  the  apparent  viscosity  above  4Q°F,  is  believed  to  be  a 
result  of  weathering,  for  the  effect  of  the  shear  rate  on 
the  apparent  viscosity  was  not  determinate,  as  its  value 
increased  with  time.  The  scattering  of  the  data  below 
40°F  is  believed  to  be  a  result  of  the  temperature  measuring 
methods.  The  value  of  the  apparent  viscosity  at  17.5°F 
for  a  shear  rate  of  3.7  sec.-1  is  reduced  by  19%  at  a  shear 
rate  of  70  sec,  \ 

Redwater  Crude 

Correlation  of  the  Apparent  Viscosity  Data 

The  smoothed  curves  of  the  three  methods  are 
illustrated  in  Fig.  (33).  The  data  of  the  MacMichael  at 
a  shear  rate  of  15  sec.-1  is  correlated  with  the  data  of 
the  Alberta  viscometer  at  60  sec.-1.  Because  the  data 
obtained  on  the  Alberta  is  more  consistent,  and  the  method 

of  measuring  the  oil  temperature  is  more  accurate,  the 

-1 

Alberta  data  at  the  shear  rate  of  60  sec.  is  believed 
to  be  more  reliable. 

The  Alberta  viscometer  data  is  approximately  20^ 
higher  than  the  pilot  pipeline  from  40°F  to  70  F.  Because 
of  the  high  gravity  of  this  crude,  the  discrepancy  may  be 
a  result  of  the  crude  being  weathered  either  prior  to,  or 
during  the  tests.  The  increasing  spread  between  the  pilot 
pipeline  and  Alberta  data  below  40°F  is  believed  to  be 
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a  result  of  the  pseudoplastic  behavior  and  the  large 
variation  in  the  rate  of  shear.  This  data  was  not 
obtained  on  the  same  sample  of  Redwater  crude,  and  the 
history  of  either  sample  is  unknown. 

CONCLUSION 

The  Saskatchewan  blend  was  found  to  be  a 
Newtonian  liquid  for  the  range  of  shear  rates  investigated 
on  the  pilot  pipeline.  At  60°F  to  80°F,  the  pipeline  data 
was  ~50%  to  45%  high,  because  of  the  high  thermal  gradients 
necessary  to  keep  the  oil  at  a  constant  temperature. 

The  MacMichael  data  indicated  that  the  oil  was 
pseudoplastic  over  the  full  temperature  range,  but  this 
may  have  been  a  result  of  an  instrumental  error.  The 
Alberta  viscometer  indicates  that  the  oil  is  pseudoplastic 
below  40°F,  and  Newtonian  above  this  temperature. 

A  composite  curve  of  the  Alberta  data  above 
40°F,  and  the  pipeline  data  below  40°F,  exhibits  a  smooth 
curve  on  the  ASTM  viscosity  temperature  graph  paper  (Fig. 34). 
This  linear  relationship  may  be  an  indication  that  the  oil 
exhibits  Newtonian  characteristics  above  4C°F.  The  value 
of  the  rate  of  shear  above  which  there  is  no  effect  of  the 
shear  rate  on  the  viscosity  may  lie  between  50  sec.“^  and 
400  sec.”1  for  this  blend. 

Combined  Blend 

On  the  pilot  pipeline,  this  crude  exhibited  no 

pseudoplastic  effects  for  the  range  of  shear  rates  investi¬ 
gated.  At  80°F  the  high  thermal  gradient  may  have  resulted 
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FIG.  3  4 

THE  SAYBOLT  VISCOSITY  OF  THE  OILS 
FROM  THE  DATA  OBTAINED  ON 
THE  ALBERTA  VISCOMETER  AND  PILOT 
PIPE  LINE 

- o-SASK.  BLEND  ALTA . 8  PILOT  PIP  E  LINEVISCS. 

- A” COMBINED  BLEND  ALTA. 8  PILOT  PIPE  LINE  DATA 

- 9-FOSTERTON  CRUDE  ALTA. 8  PILOT  PIPE  LINE  DATA 

--□-PEMBINA  CRUDE,  PI  LOT  PIPE  LINE  80000  SEC-1 
— V~ PEMBINA  CRUDE, PILOT  PIPE  LINE  3,000  SEC 
'/^“PEMBINA  CRUDE,  ALTA.  VISC.  6  4  SEC-' 

•  -  REDWATER  CRUDE  ,(MASUDA)PILOT  PIPELINE 
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in  the  apparent  viscosity  being  an  estimated  35%  to  40% 
higher  than  the  isothermal  value. 

The  data  obtained  on  the  MacMichael  indicates 
pseudoplastic  behavior  for  this  crude  below  70°F.  The 
scattering  of  data  about  the  smooth  curve  may  be  a  result 
of  the  method  used  to  determine  the  average  oil  temperature, 
or  an  instrumental  error. 

The  data  of  the  Alberta  viscometer  indicates 
that  this  blend  is  pseudoplastic  below  55°F.  The  data  of 
the  Alberta  viscometer  above  55°F  is  believed  to  be  more 
reliable  than  the  pilot  pipeline  data,  so  that  the  combina¬ 
tion  of  this  data  and  the  pilot  pipeline  data  below  45°F 
may  represent  the  apparent  viscosity  curve  for  an  infinite 
ra.te  of  shear  as  in  Fig.  (34). 

Fosterton 

The  apparent  viscosity  data,  obtained  from  the 
laminar  flow  tests  on  the  pilot  pipeline  indicates  no 
pseudoplastic  behavior.  As  it  was  not  possible  to  obtain 
turbulent  flow  at  the  high  temperatures,  the  high  thermal 
gradient  may  have  caused  the  apparent  viscosity  data  at 
70°F  to  75°F  to  be  35%  to  40%  higher  than  the  isothermal 
data. 

The  variation  of  the  apparent  viscosity  data 
determined  by  the  MacMichael  viscometer  from  80°F  to  20°F 
is  believed  to  be  the  result  of  an  instrument  error. 

The  data  of  the  Alberta  viscometer  indicates 
that  this  crude  is  pseudoplastic  below  45°F.  ihe  variation 
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in  the  apparent  viscosity  exhibited  could  be  caused  by  a 
3.0°F  change  in  temperature.  It  is  known  that  the  oil  temp¬ 
erature  cannot  be  changed  this  much  in  the  short  time 
required  to  cover  the  full  range  of  shear  rates  possible e 

A  composite  curve  of  the  Alberta  data  above  25°F 
and  the  pilot  pipeline  data  below,  plotted  on  the  ASTM 
graph  paper  Fig.  (34-),  indicates  that  this  curve  may 
represent  a  Newtonian  liquid  to  20°F,  and  non-Newtonian 
below  this  temperature. 

Pembina  Crude 

The  apparent  viscosity  data  of  the  MacMichael 
viscometer  was  obtained  over  a  small  range  in  the  rate 
of  shear.  Because  the  rate  of  shear  history  was  a  maximum 
of  15  sec."\  the  data  was  high.  The  MacMichael  is 
insufficiently  dampened  for  this  crude  at  the  high 
temperatures  and  at  the  low  temperatures  there  is  no  way 
of  accurately  determining  the  average  oil  temperature. 

Where  possible  the  shear  rate  history  of  the 
Alberta  data  was  64  sec,’\  and  this  resulted  in  a  better 
correlation.  Below  45°F,  the  data  indicates  that  the  shear 
stress  in  the  annulus  is  not  high  enough  to  completely  shear 
all  the  oil.  The  apparent  viscosity  data  correlates  rather 
well  with  the  data  obtained  on  the  pilot  pipeline.  The 
plot  of  the  smoothed  data  in  Fig.  (34)  for  the  60  sec.-1 
curve  results  in  two  points  of  inflection,  with  one  at  70°F 
and  the  other  at  30°F.  The  reason  for  two  transition  tempera¬ 
tures  for  a  pseudoplastic  crude  is  not  known. 
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A  large  number  of  tests  were  made  on  the  Pembina 
crude  over  a  wide  range  of  shear  rates.  The  poor  correla¬ 
tion  obtained  with  the  data  is  believed  to  be  a  result  of  the 
pseudoplastic  and  thixotropic  behavior  of  the  crude.  It 
may  have  been  possible  to  obtain  a  correlation  of  the  data 
of  the  three  test  lines,  but  the  variation  in  the  shear  rates 
of  flow  from  the  two  pumps,  results  in  a  large  variation  in  the 
shear  rate  history.  It  was  very  difficult  to  maintain  a 
constant  rate  of  flow  at  the  low1  pressure  drops  necessary 
for  low  rates  of  shear. 

In  the  turbulent  region  of  flow,  the  variation 

-1  «1 

of  the  apparent  viscosity  at  10,000  sec.  to  80,000  sec. 
is  small,  indicating  that  this  may  be  the  rate  of  shear 
above  which  there  is  no  effect  on  the  apparent  viscosity. 

Redwater  Crude 

No  pseudoplastic  behavior  was  indicated  by  the 
pilot  pipeline  data  obtained  by  Masuda  (19).  The  effect 
of  thermal  gradients  would  be  very  small  at  the  high  tempera¬ 
ture  because  turbulent  flow  was  obtained. 

The  data  obtained  with  the  MacMichael  visco¬ 
meter  was  erratic,  as  the  instrument  is  seriously  under¬ 
damped  for  such  a  low  viscosity  liquid  as  the  Redwater 
crude.  This  crude  exhibits  pseudoplastic  effects  below 
40°F  on  the  Alberta  viscometer,  but  the  variation  in  the 
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apparent  viscosity  is  small  over  the  range  of  shear  rates 
used.  It  is  believed  that  rapid  weathering  of  the  crude 
resulted  in  the  determination  of  high  values  for  the 
apparent  viscosity  above  60°F  in  the  Alberta  viscometer. 

The  Pressure  Effect  Coefficient  of  the 
Apparent  Viscosity 

The  data  obtained  for  the  Combined  blend,  the 
Fosterton  crude,  and  the  Pembina  crude,  are  plotted  in 
Fig.  (35).  In  general,  the  correlation  of  the  data  at 
high  pressures  was  rather  poor,  as  the  average  deviation 
varies  from  10%  to  30%.  This  may  be  a  result  of  low 
pressure  drops,  and  a  larger  relative  error  in  the 
pressure  drop. 

The  smoothed  data  for  the  Combined  blend  and 
Fosterton  crude  indicates  that  the  effect  of  pressure 
on  the  apparent  viscosity  decreases  with  increasing 
temperatures.  The  Pembina  data  was  so  scattered  that 
no  trend  could  be  determined. 

Pilot  Pipe  Line 

In  the  Pilot  Pipeline,  it  is  possible  to  obtain 
a  wide  variation  in  the  rate  of  shear  and  the  range  of 
temperatures  is  from  0°F  to  80°F  plus.  For  highly  viscous 
oils  and  at  high  rates  of  flow,  the  cooling  capacity  is 
found  to  be  inadequate.  As  a  result  it  was  necessary  to 
operate  with  a  large  temperature  gradient  from  the  oil  to 
the  bath,  which  has  resulted  in  increased  values  of  tne 
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friction  factor  and  the  apparent  viscosity.  To  better  investigate 
the  pseudoplastic  behavior  of  oils,  some  method  of  controlling  and  measur¬ 
ing  the  instantaneous  rate  of  flow  is  required.  The  use  of  thermocouples 
in  contact  with  the  flowing  oil  may  improve  the  accuracy  of  the  temperature 
measurement . 

MacMichael  Viscometer 

The  operation  of  the  MacMichael  viscometer  is  straight¬ 
forward  and  quick.  The  ease  with  which  the  torsion  wires  can  be  changed 
results  in  a  system  with  an  'uncontrollable  movement  of  rhe  bob  from 
side  to  side,  which  results  in  large  variations  for  low  viscosity 
liquids.  Because  there  is  no  means  of  dampening  the  oscillations  of 
the  bob,  the  instrument  is  limited  to  the  more  viscous  liquids.  The 
range  of  shear  rates  is  presently  limited  by  the  speed  of  the  motor, 
and  the  dimensions  of  the  cup  and  bob. 

Alberta  Viscometer 

It  is  possible  to  determine  the  angular  velocity 
quickly  and  accurately,  and  the  average  oil  Temperature  can  be 
measured  within  2.0°F  il°F.  The  range  of  shear  rates  is  theoretically 
as  high  as  152  sec.__L,  but  the  instability  of  the  bob  has  been  the 
limitation.  There  has  been  some  indication  that  the  annulus  is  too 
wide,  so  that  all  the  highly  pseudoplastic  liquid  in  the  annulus  is 
not  being  sheared.  The  dampening  of  the  bob  is  inadequate  for  light 
crudes,  such  as  the  ftedwater,  especially  at  temperatures  where  the 
weathering  is  very  serious. 
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Either  a  reduction  of  the  annular  clearance,  or 
the  installation  of  a  bearing  on  the  bob,  may  increase  the 
range  of  shear  rates  possible.  The  method  of  measuring  the 
temperature  may  be  improved  by  the  installation  of  more 
thermocouples,  and  the  thermal  gradients  may  be  lowered  by 
effectively  insulating  the  bob  and  top  of  the  cup. 
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NOMENCLATURE 

constant  in  Williamson  equation 


sec . 


-1 


coefficient  of  thixotropic  breakdown  with  time 
lb  ./ft .  sec2. 


A  &  B  constants  in  the  Von  Karman  friction  factor 
equation 


dimensional  conversion  constant 
cubic  feet  per  minute 
specific  heat  BTU/lb.°F 

coefficient  of  pressure  effect  on  viscosity 
%/l00  psi. 

pipe  diameter  ft. 

sum  of  the  friction  losses  To.f  ft ./lb. 

loss  in  the  shearing  force  per  unit  area  per 
unit  increase  in  the  shear  rate  It /-ft. sec. 

F  increase  in  yield  value  per  unit  decrease  in 

the  coefficient  of  rigidity  sec."1 

f  friction  factor  for  isothermal  flow 

#  * 

f^jf^jf^  non-isothermal  friction  factors 


c 

cf  m. 

CP 

C 

D 

dF 

E 


Sc 

ha 


J 

k 

K 

1 

In 

L 


acceleration  of  gravity  ( lb . lb„ f ) (ft ./sec .  ) 

equivalent  depth  of  immersion  in.,  and  cm. 

l/y  the  ratio  of  the  mixing  length  to  the 
distance  from  the  wall 

some  function  of 

thermal  conductivity  BTU . ft . /ft ,^hr . °F 

torsion  wire  constant  for  the  Alberta  in. cps ./sec . °A 

for  the  MacMichael  cm. cps. /see. 

mixing  length 

natural  logarithm 

length  of  conduit  section  ft. 


' 


. 


„ 


. 


■ 


.  .  . 


. 

■ 

. 


■ 


99 
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n 

N 

%e 
A  P 

Pr 

Q 

q 

r 

R 

S 

t 


v 

V 
w 

X 

y* 

y 

Y 


torsion  wire  torque  lb.^ft. 
exponential  constant 
revolutions  per  second 
Reynolds  number 
pressure  drop  lb.^/in.2 
Prandtl  number 
flow  rate  ft.^/mi n. 

constant  or  complex  number  for  exponential 
form  of  equation 

radius  as  measured  from  the  axis  ft. 
radius  of  the  cup,  bob  (in.),  and  pipe  (ft.) 
R2/Rl 

temperature  °F 

temperature  in  the  laminar  sublayer  °F 
temperature  of  the  oil  °C,  and  °F  respectively 

cumulative  time  of  shear  min, 
shear  stress  lb.^/ft.^ 

shear  stress  at  the  yield  value  lb . f/ft . 2 
apparent  yield  value  lb.^/ft." 

point  velocity  ft. /sec. 
bulk  velocity  ft. /sec. 
angular  velocity  2TfN  sec.  ^ 
distance  above  a  datum  ft. 

ypu*-/u  group  similar  to  the  Reynolds  number 

distance  as  measured  from  the  wall  of  a  pipe  ft. 

constant  or  complex  number  characteristic  of 
a  pseudoplastic  liquid 
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u * 

friction  velocity 

u+ 

u/u^  point  velocity  to  friction  velocity 

(dv/dr ) 

rate  of  shear,  shear  rate,  or  velocity 
gradient  sec.“l 

Greek 

Letters 

a 

coefficient  of  rigidity,  or  plastic  viscosity 
lb ./ft . sec . 

s* 

coefficient  of  viscosity,  or  viscosity  lb, /ft. sec. 

y^lp,y^hp 

viscosity  at  low  and  high  average  pressures 
respectively  lb. /ft. sec. 

a 

apparent  viscosity  as  defined  by  the  Hagen- 
Poiseuille  law  lb. /ft. sec. 

a 

apparent  viscosity  for  a  constant  rate  of 
shear  lb. /ft. sec. 

/a 

apparent  viscosity  defined  by  the  equation  for 
the  rotational  viscometer  lb. /ft. sec. 

/** 

apparent  viscosity  as  defined  by  the  Darcy-Weisbach 
equation  lb. /ft. sec. 

apparent  viscosity  at  zero  rate  of  shear  lb. /ft. sec 

apparent  viscosity  at  an  infinite  rate  of 
shear  lb. /ft. sec. 

”ir 

2 

kinematic  viscosity  ft. "/see. 

P 

density  of  flowing  liquids  lb. /ft. 

net  deflection  of  the  rotational  viscometers 
in  °A,  and  °M 

apparent  fluidity  ft.  sec. /lb 

Subscripts 


1,2 

on  a  conduit  upstream  and  downstream  sections;  on 
the  rotational  viscometer,  the  bob  and  cup  surfaces 
respectively 

r 

the  condition  at  any  radius  as  measured  from  the 
axis  ft. 
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w  the  condition  at  the  wall  of  a  conduit  ft. 

y  the  condition  at  any  distance  as  measured  from 

the  conduit  wall  ft* 
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APPENDIX  A 

GENERAL  GLASSIFICATION  AND  FLOW  BEHAVIOR  OF  LIQUIDS 

(for  nomenclature  see  page  ^8 ) 

Newtonian  Liquids 

Liquids  are  classified  according  to  the  general  type  of  relation- 
i  ship  determined  between  the  shear  stress  and  the  rate  of  shear  in  laminar 

■  IB 

flow  (32) .  The  Newtonian  liquids  are  characterized  by  a  constant  value 
of  the  viscosity  which  is  defined  as  the  ratio  of  the  shear  stress  to 
the  shear  rate.  A  typical  shear  stress-shear  rate  diagram  (often  referred 
to  as  a  consistency  curve)  for  such  a  material  is  illustrated  in  Fig.  36  (a) . 
The  value  of  the  viscosity  at  any  temperature  is  the  slope  as  follo\/s(l), 

>  S  -Tgc  /  (dv/dr) 

where 

v  s  point  velocity  ft. /sec. 

Some  examples  of  Newtonian  liquids  are ;  light  lubricating  oils, 
stove  oil,  and  some  crude  oils  over  limited  ranges  of  temperature. . 

Plastic  Liquids 

The  Bingham  plastics,  or  as  they  are  more  commonly  called  plastic 
liquids,  exhibit  rigidity  at  zero  rate  of  shear,  and  in  laminar  flow 
the  shearing  stress  in  excess  of  the  yield  value  Ty,  is  directly  proport¬ 
ional  to  the  shear  rate ( l) .  The  equation  of  the  typical  flow  curve  of  a 
plastic  liquid  illustrated  in  Fig. 36(b)  is 

ry  *  -(T-Ty)gc  /  (dv/dy)  (19) 

coefficient  of  rigidity  or  plastic 

viscosity  lb. /ft. sec. 

2. 

Tvs  yield  value  lb.f/ft. 

^  /  2 
T  —  shear  stress  at  the  shear  rate  (dv/o y)  Ib.^/ft. 


where 
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T 

tan 

VISCOSITY  T 

cotan  P 

COEFFICIENT  OF 

RIGIDITY 

tan  •©■ 

.  /  APPARENT 

/  VISCOSITY 

/ 

,  t  a  n  -e 

^  APPARENT  VISCOSITY 

( dv/dy) 

a)  N  E WTO N JAN 

(dv/dy) 

b)BINGHAM  PLASTIC, or  PLASTIC 

T 

// 
/ /o 

/tan©  APPARENT 

VISCOSITY  at 

// 

// 

indicated  (dv/dy) 

T 

(dv/dy) 

c)  PSEUDOPLASTIC 

(dv/dy) 
dl  DILITANT 
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ihe  flow  behavior  of  a  plastic  liquid  can  be  defined  by  an  apparent 
viscosity(l,3) .  The  apparent  viscosity  is  the  ratio  of  the  shear  stress 
to  the  shear  rate,  at  one  shear  rate,  but  the  value  varies  with  the 
change  of  the  shear  rate.  The  coefficient  of  rigidity  has  also  been 
defined  as  the  plastic  viscosity  (12),  and  this  value  is  a  constant 
for  a  plastic  liquid.  Examples  of  plastic  liquids  are  suspensions  such 
as  sand  in  water,  bentonite  in  water,  or  barytes  in  water(l,12). 
Pseudoplastic  Liquids 

Pseudcplastic  liquids  in  laminar  flow  exhibit  a  non-linear  relation¬ 
ship  between  the  shear  stress  and  shear  rate.  For  a  large  number  of 
pseudoplastic  liquids,  the  apparent  viscosity  asymptotically  approaches 
a  constant  value  with  increasing  shear  rates(l,36),  to  form  a  consistency 
curve  as  illustrated  in  Fig.  36(c).  Only  empirical  equations  have  been 
determined  to  express  the  shear  stress-shear  rate  relationship  for 
pseudoplastic  liquids  in  laminar  flow.  The  form  of  a  general  equation 
advocated  by  Reiner  (13),  to  describe  this  type  of  curve  is 

(dv/dy)  -  aTbgc  12.0} 

where 

a,b,»  dimensionless  constants,  the  value 
of  which  are  determined  by  the 
consistency  curve 

These  constants  cannot  be  associated  with  any  characteristic  oj. 
flow  behavior,  and  are  applicable  to  only  one  liquid  at  one  temperature. 

The  Willianson  equation  (36)  has  found  wide  application  for  a 
large  number  of  pseudoplastic  liquids,  for  which  the  apparent  viscosity 

practically  approaches  an  asymptotic  value  at  high  rates  of  shear. 
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Cne  form  of  this  equation  applicable  to  the  laminar  flow  region  is 
=  [Tygc  /(a1/  dv/dr)J  /  u^ 

The  limit  of  this  equation  as  (dv/dr)  approaches  zero  is(36) 

limits  =  >iD  =  T'  gc  /  (a)  /  U<x>  (U) 

dv/  dr  _>  O 
where 

:  apparent  viscosity  at  zero  rate  of  shear 
lb. /ft.  sec. 

^oo  =  apparent  viscosity  at  an  infinite  rate  of 
shear  lb. /ft.  sec. 

a1  z  constant,  as  determined  by  the  curvature  near 
the  origin  sec* 

Ty  =  intercept  on  the  shear  stress  axis  determined 
by  the  extrapolation  of  the  linear  part  of 
the  apparent  viscosity- shear  rate  curve. 

The  constants  aJTy,  and  are  functions  of  the  system,  and  the 

. 

temperature,  and  not  of  the  rate  of  shear,  or  the  dimensions  of  the 

tube (36) . 

The  decrease  in  apparent  viscosity  with  increasing  shear  rate 
as  exhibited  by  pseudoplastic  liquids  has  been  attributed  to  the  decrease 
of  the  internal  friction  caused  by  the  long  molecules  being  in  a  random 
distribution  at  the  low  rates  of  shear,  or  at  rest.  Under  increasing 
rates  of  shear,  these  long  molecules  are  aligned  to  the  streaming 
flow,  resulting  in  a  decreased  internal  resistance  between  the  flowing 
layers (12) .  Examples  of  this  are  the  long  molecules  synthesized  for 
artificial  rubber,  wax  chrystals  in  crude  oil,  bitumens,  and  resinous 
materials. 


Dilatant  Liquids 

A  dilatant  liquid  exhibits  an  "inverted*1  shear  stress  -  shear  rate 
relationship,  resulting  in  an  increased  apparent  viscosity  with  an 
increased  shear  rate,  ihis  flow  behavior  is  a  result  of  the  liouid  and 
solid  particles  occupying  a  minimum  volume  when  at  rest.  On  the  applic¬ 
ation  of  a  shearing  force,  it  is  necessary  for  the  mixture  to  increase 
in  volume  before  the  particles  can  slide  past  each  other.  This 
behavior  gives  a  high  internal  friction,  which  increases  with  increasing 
shear  rates.  A  typical  curve  for  such  materials  as  sand  and  water, 
clay  and  water  is  illustrated  in  Fig.  36(d). 

Plastic  and  Pseudoplastic  Thixotropic  Liquids 

The  internal  resistance  to  flow  of  a  thixotropic  plastic  or 
pseudoplastic  liquid  is  a  function  of  both  the  shear  rate,  and  the 
previous  shear  rate  history.  The  variation  in  the  consistency  is 
attributed  to  an  internal  structure  often  referred  to  as  the  thixotropic 
structure.  It  has  been  determined  experimentally  that  the  thixotropic 
structure  may  be  reduced  to  a  minimum  value  by  placing  the  liquid  under 
a  constant  rate  of  shear  for  a  sufficient  period  of  time.  The  minimum 
consistency  value  represents  the  equilibrium  value  at  which  the  break¬ 
down  by  shearing  forces,  and  the  buildup  of  the  thixotropic  structures 
are  equal.  At  this  equilibrium  value,  the  only  way  to  reduce  the  value 
of  the  consistency  is  by  an  increased  rate  of  shear.  This  is  illust¬ 
rated  in  Fig.  37(a)  and  (b) ,  for  thixotropic  plastics  and  pseudoplastics 
respectively. 

For  thixotrpic  plastics,  the  linear  downcurves  or  thixotropic  levels 
obtained  at  one  rate  of  shear  may  be  plotted  against  the  logarithm  of 
the  respective  cumulative  time  as  in  Fig.  37.  This  results  in  a  linear 
relationship  of  the  plastic  viscosity  and  logarithm  of  the  cumulative 
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time  of  shear,  the  slope  of  which  is  defined  as  the  coefficient  of 
thixotropic  breakdown  with  time (12), 

B  =  (t?i-t\2)/lnt2/t1  (2^0 

where 

B  -  coefficient  of  thixotropic  breakdown  with 
time  Ib/ft.  sec?" 

ti,  t^  ”  the  cumulative  time  interval  for  the 
plastic  viscosity  values  Tj|,  and 
respectively 

For  the  thixotropic  plastics  it  is  possible  to  determine  two  other 
coefficients,  for  better  defining  the  liquid  properties.  The 
coefficient  of  thixotropic  breakdown  with  increasing  shear  rate  may 
be  computed  as  follows (12), 

E  =  2(T]1  -  n2)  (In  wjj  /  w^)  (23) 


where 

E  s  loss  in  shearing  force  per  unit  area 

per  unit  increase  in  shear  rate  lb. /ft.  sec. 

Wp  w^  -  angular  velocity  for  tji  and  respectively  sec:1 
The  other  coefficient  which  may  be  determined  for  a  thixotropic  plastic 
is  the  increase  in  the  yield  value  for  a  decrease  in  the  plastic  viscosity. 

F  &  (T2  ”  T-^J  /  "  ^2/  (24) 


where 


F  s  increase  in  yield  value  per  unit  decrease 
in  the  coefficient  of  rigidity,  sec. 

T  ,  T_  s  intercept  in  the  shear  stress  axis  by 
the  thixotropic  level  for  the  plastic 


viscosity  values  7|][  and  respectively. 


(dv/dy) 


Thixotropic 
Level  / 


Downcurve 


(dv/dy) 


pcurve 


7  T  T 

FIG.  37 

a)TH,X'CTROPIC  PLASTIC  b)  THIXOTROPIC  PSEUDOPLASTIC 


( dv/dr) 


Plastic 

Viscosity 


a)  THIXOTROPIC  BREAKDOWN 
UNDER  A  CONSTANT  RATE  OF  SHEAR 


b)DECREASE  IN  THE  PLASTIC 
AT  A  CONSTANT  RATE  OF 
SHEAR 


FIG.  38 
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These  coefficients  which  can  only  be  determined  for  the  thixotropic 
plastic  liquids,  make  it  possible  to  better  determine  and  define  the 
thixotropic  properties,  for  comparison  and  evaluation (12) . 

No  analysis  of  thixotropic  pseudoplastics  has  been  performed  similar 
to  that  above,  as  the  dowcurve  obtained  is  not  a  linear  relationship. 
Because  of  the  complexity  of  the  variables  connected  with  these  liquids 
an  experimental  qualitative  approach  has  been  used  (l) .  By  flowing  a 
thixotropic  pseudoplastic  liquid  through  pipes  of  varying  diameters 
and  lengths  it  is  possible  to  determine  the  scale  up  factor  for 
applying  laboratory  data  to  practical  applications.  The  information 
obtained  may  be  used  to  estimate  the  size  of  lines  and  probable 
pressure  drop(l). 

There  have  been  two  theories  advocated  to  account  for  the 
thixotropic  behavior  of  these  liquids(12).  If  the  molecules  in  the 
liquid  possess  polar  charges,  long  chains  of  molecules  will  be  formed 
creating  a  gell-like  structure  with  a  yield  value.  Under  a  shearing 
action,  the  long  chains  will  be  broken  gradually,  and  the  internal 
friction  of  flow  will  be  reduced  with  time.  The  slow  buildup  of  the 
consistency  with  a  reduction  of  the  shear  r  ate  is  the  result  of  the 
time  necessary  for  the  molecules  to  realign  the  polar  charges  as  a 
result  of  Brownian  movement.  The  thixotropic  behavior  of  pseudoplastics 
is  believed  to  be  a  result  of  long  molecules  being  alligned  under  the 
action  of  shearing  forces,  and.  the  slow  buildup  is  a  resuxt  of  the 
time  necessary  for  the  long  chain  molecules  to  reach  the  initial  random 
order (1,12) . 
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Rheopectic  Liquids 

Rheopectic  liquids  are  those  which  will  set  up  or  exhibit  an 
increase  in  the  apparent  viscosity  very  rapidly  upon  being  rythmically 

shaken  or  tapped.  Some  examples  of  rheopectic  materials  are  gypsum 

m ' 

suspensions  in  water,  bentonite  sols,  and  vanadium  pentoxide  sols(l). 


EQUATIONS  OF  FLOW  OF  LIQUIDS  IN  PIPES 

Consider  a  Newtonian  liquid  flowing  isothermally  through  a 
differential  distance.  Then  the  work  effects  are  localized  and  the 
mechanical  energy  balance  for  each  element  is  (32,27), 


where 


dP/p  /  dx  /  VdV  =  -dF  ft.  Ib/lb. 
§c 


(25) 


dP/p  -  change  in  pressure  lb.fft.'/ft.  lb. 

dx  r  change  in  elevation  f^  lb.^/lb. 

2  2 

VdV/gQ  ~  change  in  the  kinetic  energy  (ft  /sec."-) 

(lb.£  sec.^/lb.ft.) 
dF  =  sum  of  friction  losses 

Assume  that  the  cross  section  is  constant,  that  the  pipe  is 
horizontal,  and  that  the  velocity  profiles  at  each  end  are  identical,  then 
dP/p  =  dF  (26) 

The  friction  loss  of  liquid  flow  lor  a  dif ferential  length 
dL  is  defined  by  the  following  relationship (32) 
dF  -  fdLV2/2Dgc 

Equating  equation  (26)  and  (27) ,  the  expression  for  the  head  loss 
is  obtained 

dP/p  -  fdL?2  /  2Dgc  ^23' 

Assume  that  D  is  constant,  and  that  the  velocity  profile  is  not  changing. 
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For  an  essentially  incompressible  liquid  flowing  through  the  finite 
section  (1)  to  (2)  with  no  exchange  of  work,  the  integration  of  this 
relationship  is  i. 

AP  =  V2p  VfdL 


(29) 


2°gc  % 

If  the  friction  factor  f  is  very  nearly  constant,  the  sum  of  the 
friction  losses  is  expressed  by  the  integrated  form  of  the  equation 
which  is  known  as  the  Darcy-Ueisbach  equation  (9,27) 

a?  ;  fL  V?  /  2Dg0  (30) 

The  preceding  analysis  is  based  on  the  thermodynamic  and 
empirical  approach  to  fluid  flow. 

The  analysis  of  the  shearing  stresses  evident  between  the 
cylindrical  layers  leads  to  an  improved  Understanding  of  the  internal 
flow  mechanism.  The  maximum  shear  stress  in  the  liquid  may  be  calcul¬ 
ated  by  taking  a  force  balance  on  the  ring  of  fluid  in  contact  with 
the  pipe  wall  (4,19,31), 


Tw  s  &PD  /  4L 


(31) 


where 

2 

Tw  -  shearing  stress  at  the  wall  lb.^/ft. 

By  solving  for  the  shearing  forces  on  any  other  cylinder  of 
flowing  liquid,  it  can  be  shown  that  Tr,  the  shear  stress  at  any  radius 
r,  is  linear  with  the  radius(lO),  and 

Tr/Tw  =  r/R  (32) 

The  shear  stress  varies  from  a  maximum  at  the  pipe  wall  to  zero 
at  the  axis.  This  is  true  no  matter  what  the  properties  of  the  liquid 


are,  or  the  type  of  flow. 


Laminar  Flow  of  Newtonian  Liquids 

For  a  Newtonian  liquid  in  laminar  flow,  the  relationship  of  the 
shear  stress  for  any  layer,  to  the  rate  of  shear  is  determined  by 
equation  (4)  page  (3) .  Substituting  the  general  expression  for  T  in 
terms  of  pressure  drop  and  pipe  length,  and  rearranging  gives  (27), 

dv  :  -g  A  P  rdr  /juL2.  (33) 

The  constant  of  integration  is  evaluated  at  r  =  R  where 
vr  -  o  and  the  integrated  form  is 

vr  r  "gcA?(lT-r2)  /  4juL 

The  average  or  bulk  velocity  over  the  entire  cross  section  may 
be  shown  to  be  one-half  the  maximum  velocity  which  occurs  at  the  center, 
as  the  plot  of  the  point  velocity  describes  a  parabola.  The  substit¬ 
ution  of  D  =  2R  into  the  expression  for  the  bulk  velocity,  yields  the 
well  known  relationship  called  the  Hagen-Poiseuille  Law(9), 

AP  =  32>LV  /  gcD?  (34) 

By  comparing  this  equation  with  the  Darcy-weisbach  equation,  the 
value  of  f  for  isothermal  laminar  flow  of  Newtonian  liquids  is  found 
to  be 

f  -  64^1  /  DVp 
=  64  /  %e 

This  has  been  conclusively  proven  by  experimental  data  (4,32), 
obtained  both  in  rough  and  in  smooth  pipes  for  steady  state  laminar 
flow. 


By  the  substitution  of  the  Hagen-Poiseuille  equation  for  pres¬ 
sure  drop  into  the  ex  pression  for  the  shear  stress  at  the  wall, 
eqN. (31) ,  the  following  expression  for  the  shear  rate  at  the  wall  is 


derived 
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(dv/dr)w  r  -8V/t) 

This  formula  shows  that  the  shear  rate  in  laminar  flow  is 
directly  proportional  to  the  volume  rate  of  flow  for  one  pipe. 

Turbulent  Flow  of  Newtonian  Liquids 

When  a  liquid  enters  a  pipe,  a  boundary  layer  developes  from  the 
leading  edge,  and  results  in  a  velocity  gradient  perpendicular  to 
the  stream  flow.  At  the  outer  edge  of  the  boundary  layer,  the  velocity 
gradient  is  zero  as  at  the  axis  of  a  pipe  in  laminar  or  turbulent 
flow.  The  range  of  flow  rates  defined  by  the  Reynolds  number  from 
0  to  2,000  represents  the  laminar  flow  region  for  most  systems,  and 
the  boundary  layer  flow  pattern  is  governed  by  the  viscous  forces 
only.  The  range  of  Reynolds  number  from  2,000  to  4,000  is  called  the 
transition  region  and  the  boundary  layer  changes  from  a  laminar  to 
a  turbulent  flow  pattern  in  which  the  velocity  profile  is  governed 
largely  by  the  momentum  exchange.  In  some  cases  the  transition  region 
1  starts  at  a  Reynolds  number  as  low  as  1,000(24,26),  but  in  most  cases 
it  occurs  at  a  value  of  2,000  or  higher  as  defined  by  the  first 
indication  of  sinousoidal  flow  at  the  axis  which  is  not  dampened  with 
l.ength(26) .  As  the  Reynolds  number  is  increased  above  2,000  the 
Turbulent  core,  characterized  by  the  sinousoidal  flow  pattern,  increases 
in  area  with  increasing  Reynolds  number.  In  fully  developed  turbulent 
flow  the  boundary  layer  consists  of  a  thick  turbulent  core,  a  thin 
buffer  layer,  and  a  very  thin  laminar  sublayer. 

In  the  turbulent  flow  region  it  is  assumed  that  the  shear  stress 


at  any  radius  In  a  pipe  can  be  represented  by  the  sum  of  a  viscous 
and  a  turbulent  agent  as  follows (4) 5 


/  drl  /  (Pi2  /  S  )  (dv/dr)2 

J  0  y 


(35) 


where 


Ty  I  Shearing  stress  at  any  distance  y  from  the 
pipe  wall,  lb.^/ft, 

1  -  may  be  interpretted  as  a  mixing  length 
The  first  term  represents  the  viscous  forces  in  effect  in  the 
laminar  sublayer  and  buffer  layer.  The  second  term  represents  the 
shearing  forces  resulting  from  the  momentum  interchange  in  the  buf¬ 
fer  layer  and  turbulent  core.  The  second  terms  analogous  to  the  kin¬ 
etic  energy  exchanged  from  layer  to  layer,  if  the  average  velocity 
components  in  two  dimensional  flow  can  be  represented  by  the  product 
of  the  velocity  gradients  and  the  mixing  lengths.  In  the  laminar 
sublayer  where  the  viscous  forces  control  fluid  flow,  the  first  term 
is  applicable,  and  the  second  term  is  zero.  From  the  buffer  layer 
to  the  turbulent  core,  the  effect  of  the  viscous  forces  become  so 
small  that  they  may  be  neglected.  For  this  central  core,  equation 
(35)  reduces  to 

V  TvS(/p  -  1(dUy  /dr)  (36^ 

z  u*  (37) 

The  quantity-^/ Tw  gc  /p  has  the  dimensions  of  velocity  and  has 
been  given  the  name  of  friction  v  elocity  and  the  symbol  u  .  This 

ij/t- 

equation  relates  the  friction  velocity  u  to  the  mixing  length,  and 
the  velocity  gradient  at  a  point  y. 

From  the  Darcy-Weisbach  equation  (4) >  and  the  expression  for  T^, 
it  can  be  shown  that 

Tw  =  fpV2  /  8gc  (38) 

upon  rearranging  and  introducing  the  definition  of  the 
friction  velocity  equation  (37) 

u*  z  vm 


(39) 
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The  friction  velocity  u  is  therefore  seen  to  be  proportional  to 
the  product  of  the  bulk  velocity  times  the  square  root  of  the  friction 
factor. 

Prandtl  found  from  the  investigation  of  experimental  results 
that  the  ratio  of  the  mixing  length  to  the  radius  was  some  function 
of  the  ratio  y/r.  To  express  this  more  concisely,  a  new  relation¬ 
ship  has  been  introduced  (7) , 

I  r  1/y  (40) 

Then  substituting  into  equation  (36) 
u*r  Iy(dvy/  dr) 

This  equation  does  not  apply  at  the  axis  where  Ty  is  zero  or  in  the 
laminar  sublayer  where  the  viscous  forces  control  the  flow  pattern. 
Rearranging  and  integrating  for  the  velocity  VT  r 

vrv  I  u*ti  lny  /  o)  (41) 

y  I 

where 

c  s  is  a  constant  of  integration.  It  may  be  shown 
(4,7)  that  the  constant  of  integration  is 

c  -  I*  -  1  In  (/y pu*)  (42) 

I 

where 

I'=  a  constant 

Because  the  viscous  forces  have  been  neglected,  equation (41) 
gives  a  velocity  at  the  walls  equal  to  minus  infinity,  instead  of 
zero,  and  similarly  the  shear  rate  at  the  axis  is  never  zero.  By 
determining  the  constants  of  this  equation  from  experimental  data 
(4  ),  the  resulting  equation  for  the  point  velocity  in  the  turbulent 


core  is 

=  u*(2.5  lnypu*/AT  /5.5) 


(43) 


This  equation  may  be  written  in  a  simple  form  by  the  introd 
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uction  of  a  dimensionless  velocity  term  u^  and  a  dimensionless  position 

.  * 
term  y 

where 

u/  z  which  is  a  dimensionless  parameter  of  the  (44) 
point  velocity  and  the  friction  velocity 
y ^  Z  ypu*  which  is  a  friction  distance  parameter  (45) 

with  a  dimensionless  structure  analogous  to 
the  Reynolds  number 

For  smooth  pipes,  the  point  velocity  in  the  turbulent  core  may  be 
determined  by 

vTs  2.5  lny*V  5.5  30<$*  (46) 


Because  the  viscous  effects  are  neglected,  this  relationship  is 
only  an  approximation  in  the  buffer  layer,  where  y*  varies  from  5  to 
30.  With  rough  pipes  the  same  equation  is  used  for  the  turbulent  core 
but  the  constant  is  less  than  5.5  because  the  maximum  velocity  in 
the  laminar  sublayer  is  lower. 

In  the  laminar  sublayer,  the  viscous  forces  control  the  flow 
pattern,  and  since  this  layer  is  very  thin,  the  parabolic  curve  may 
be  represented  by  the  linear  relationship. 


Tw  z  (ju  xr  )  /  (y  gc  ) 

\r  z  (Tw  y 5c)  /  ) 

Rearranging  and  introducing  the  dimensionless  quantities 


(47) 

(48) 

x/  and  y1'" 


* 


(49) 


Where  y^  is  less  than  5,  this  relationship  is  applicable,  and 
equation  (49)  indicates  that  a  higher  laminar  sublayer  velocity  will 
be  obtained  with  a  thicker  sublayer  j  i.  e.  with  smooth  pipes,  the 
laminar  sublayer  is  thicker  than  with  rough  pipes,  so  the  bulk 
velocity  is  higher. 

The  equation  for  the  buffer  layer  is  derived  by  interpolating 
between  the  laminar  sublayer  equation  and  the  turbulent  core  equation. 
Between  y^  equal  to  5  and  30,  the  equation  is 

u^  =  5.01n/-  3.05  (50) 

The  curves  depicting  the  equation  for  the  laminar  sublayer, 
the  buffer  layer,  and  the  turbulent  core  are  reproduced  in  Fig. (39) 
and  the  complete  curve  is  called  the  Universal  Velocity  Profile. 


UNIVERSAL  VELOCITY  PROFILE  FOR  A  SMOOTH  PIPE 


The  Velocity  Gradient  in  the  Turbulent  Flow 

Through!  the  laminar  sublayer,  the  velocity  gradient  is  a 
constant  defined  as  follows 

(d  v/dr)w  --v 

=  “Tw  g  Jjx  (51) 

By  substituting  in  equation  (38)  into  equation  (51),  the  practical 
form  of  the  velocity  gradient  equation  for  the  laminar  sublayer  is 
obtained, 

(dv/dr)  z  “APD  g  /4/xL 
w  c 


In  the  buffer  layer,  the  velocity  gradient  of  the  Universal 
velocity  profile  from  r  5  to  y^z  30,  is 


(du^/dy^) z  5.0  /  y 

(52) 

and  simplifying 

(dv/dr)  z  -(5.o/y)  ^/Twgc/p 

(53) 

By  substituting  the  relation  for  the  shear  stress  T  ,  the  velocity 
gradient  becomes 

(d\r/dr)  z  _(5.o/y)\J  APDgQ/ 4jo.pL  (54) 

For  y^  greater  than  30,  the  velocity  gradient  in  the  turbulent 
core  is 

du  /dj/  -  2.5  In  3/ 

z  (2.5/y)V  Vo/P  (55) 

Simplifying 

(dv/dr)  z  -(2.5/r)V^'Ogc / bjy 

The  velocity  gradient  in  the  turbulent  core  is  independent  of 
the  relative  roughness  of  the  pipe  wall,  as  this  effect  onlj  x educe., 
the  maximum  velocity  obtained  at  the  inside  limits  of  ode  laminar 


sublayer . 


In  turbulent  flow,  the  velocity,  as  determined  by  the  Universal 
velocity  profile  is  the  average  velocity  over  a  period  of  time  at 
one  point.  It  has  been  determined  experimentally  that  at  any  one 
position  in  the  buffer  layer  or  turbulent  core,  the  instantaneous 
velocity  will  oscillate  about  the  mean  value  in  a  simusoidul  manner. 
Therefore  the  shear  rate  as  computed  above  does  not  include  the 
shear  rate  resulting  from  the  oscillatory  action  of  the  eddies  (36). 
Friction  Factor  in  Turbulent  Flow 

In  the  transition  region  of  flow,  the  plot  of  the  Darcy- 
Weisbach  friction  factor  as  a  function  of  the  Reynolds  number, 
describes  a  loop,  the  position  of  which  must  be  determined  experiment¬ 
ally  (32) .  This  region  of  flow  is  also  characterized  by  oscillatory 
pressure  variations  (24,26)  and  rapid  changes  in  the  friction  factor 
with  a  small  change  in  the  Reynolds  number. 

In  the  turbulent  flow  region,  the  average  velocity  of  flow  V, 
may  be  determined  by  the  use  of  the  following  relationship  obtained 
by  integrating  the  velocity  profiles  and  neglecting  the  volume  of 
flow  in  the  laminar  sublayer,  as  it  is  insignificant  (4) • 

V/u*  =  1.43  /  5.75  log  Ru*/xr  (57) 

where 

V  -  average  or  bulk  velocity  ft. /sec. 

“vr-  kinematic  viscosity  -  u/p  ft*  /sec. 

The  integrated  form  of  the  velocity  profile,  as  above,  was  used 
by  von  Karman  to  find  the  friction  factor  equation  for  smooth  tubes. 
The  relationship  of  the  friction  velocity  and  the  Darcy- v'/eibbach 
friction  factor  are  substituted  into  equation  (57;  and  the  rearrange¬ 
ment  gives  the  general  form  of  the  friction  facuor  equation  as  follows 
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1/y[^~z  k/ JIT /  3/^8”  (log  Ru*  /-y-  )  (58) 

The  values  of  the  constants  were  determined  by  experimental  data, 
and  the  friction  factor  equation  for  turbulent  flow  in  smooth  pipes  is 
I/VT *  2.0  log  -y/T  -  0.8  (59) 

This  equation  describes  the  friction  factor  for  a  smooth  tube  such 
as  glass,  drawn  brass,  copper,  or  lead  and  its  range  of  application 
is  from  a  Reynolds  number  of  4,000  to  3,000,000.  The  curve  describes 
the  lowest  curve  plotted  on  many  friction  factor  charts. (4,26) 

There  are  many  plots  of  the  friction  factor  versus  Reynolds 
number  which  are  the  result  of  correlating  the  friction  factor  data 
of  a  large  number  of  investigators.  One  of  the  most  complete 
consists  of  1,000  points  correlated  with  a  maximum  deviation  of  five 
percent  from  the  mean  for  liquids  with  wide  variations  in  the  viscosity 
and  density (32). 

Empirical  equations  have  been  computed  to  fit  this  experimental 
data  to  facilitate  the  solution  of  problems.  One  of  the  earliest 
was  the  equation  of  Blasius  for  smooth  tubes  which  is  valid  from  a 
Reynolds  number  of  3,000  to  100,000  (10,9,27).  It  is 

f  =  0.316A’Re)  1//4 

Another  equation  computed  from  empirical  data  on  pipe  lines  of  six 
to  twelve  inchs  diameter  transporting  crude  oil,  is  the  Heltzel 
equation ( 14) . 

0.265 

f  =  0.364/ (NR  )  2,500  <  NRe<  57,600 

e  0.188 

f  r  0.157/(NRe)  57,600  <  NRq 

The  Heltzel  equation  agrees  with  the  Blasius  equation  from  a 
Reynolds  number  of  4>000  to  above  10,000.  The  values  obtained  by 
the  von  Karman  equation  are  3$  lower,  for  the  range  of  Reynolds 
numbers  from  4,000  to  above  20,000  than  the  Heltzel  equation. 
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The  data  used  in  the  determination  of  the  Heltzel  equation,  has  been 
correlated  with  the  Stanton  and  Pannell  data(2) . 

Flow  of  Plastic  Liquids 
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From  Bingham’s  basic  equation  for  the  coefficient  of  rigidity 
for  plastic  liquids  equation  (18),  Buckingham  derived  a  general 
expression  for  the  laminar  flow  of  plastic  liquids  in  a  pipe (21). 
The  equation  in  the  simplified  form  is (21) 


v=  g0n 


32Lrj 


f  i-  A 

3 


4L  T 
Y 

Enr 


(60) 


To  develope  and  simplify  this  equation,  the  following  assump¬ 
tions  were  made*  the  slip  velocity,  or  lubricting  layer  was  assumed 
to  be  negligible  or  non-existent (12, 21) ;  there  is  a  layer  of  liquid 
which  is  flowing  in  laminar  motion  between  the  pipe  wall  and  the  bound¬ 
ary  layer  of  the  central  core;  the  central  core  or  plug  is  a  semi¬ 
liquid  structure  because  the  shear  stress  at  the  edge  of  the  plug 
in  laminar  flow  is  less  than  the  yield  value  shear  stress. 

The  application  of  the  Buckingham  equation  has  been  simplified 
by  the  use  of  dimensionless  ratios  to  describe  the  flow  pattern  of 
plastic  liquids  and  by  the  use  of  monograms  to  more  readily  solve 
the  resulting  simultaneous  equations ( 21) .  Other  methods  have  been 
advocated  for  solving  the  problems  of  flow  for  plastic  liquids  and 
these  may  be  found  in  the  literature (31, 22, 28) . 
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Laminar  Flow  of  Pseudoplastic  Liquids 

m  9 

The  pseudoplastic  liquids  are  characterized  by  a  non-linear 
relationship  of  the  shear  stress  and  the  shear  rate  in  laminar  flow, 
which  can  be  expressed  by  the  following  equation ( 13) 

z  J  (T/dv/dr)  (61) 

where  -  apparent  viscosity  Ib/ft.sec, 

j  -  some  function  of 

T  s  shear  stress  at  a  point  which  is  a  unique 
function  of  the  shear  rate  (dv/dr) 

Equation  (6l)  shows  that  the  apparent  viscosity  is  some  function 
of  the  shear  rate,  with  the  result  that  it  varies  from  a  minimum 
value  at  the  wall  of  the  pipe  to  a  maximum  at  the  axis (36).  For 
one  flow  rate  in  laminar  flow,  the  determination  of  the  absolute 
value  of  the  apparent  viscosity  or  the  shear  stress-shear  rate  entails 
the  accurate  determination  of  the  velocity  profile  across  the  annulus 
or  radius  by  an  independent  method  (13,25).  In  a  few  cases  the 
velocity  profile  has  been  determined (24)  but  generally,  most  invest¬ 
igators  obtain  an  average  value  of  the  apparent  viscosity  for  one 
flow  rate.  The  variation  of  the  apparent  viscosity  with  shear  rate  is 
determined  by  obtaining  data  over  a  wide  range  of  flow  rates(l,3,l6,36) , 
in  a  capillary,  pipe  line,  or  rotational  concentric  cylinder  viscometer. 

Since  the  flow  analysis  of  a  pseudoplastic  liquid  is  so  comp¬ 
licated,  no  rigorous  theoretical  equation  has  been  derived.  However, 
several  empirical  methods  have  been  developed (1,3)  for  the  correlation 
of  shear  stress-shear  rate  relationships  and  apparent  viscosity- shear 
rate  relationships.  One  of  these  methods (13)  assumes  that  the  consistency 
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curve  for  a  pseudoplastic  liquid  can  be  described  by  the  general 
relationship 

Z/M'  (dv/dr)  (62) 

in  which  y  is  a  simple  or  complex  number 
characteristic  of  the  liquid 

Then  for  laminar  flow  in  a  capillary  tube,  the  following 
relationship  for  the  shear  stress  and  shear  rate  may  be  derived 

GP  R/2 L)^  ->’V(Y/3)/R 
”  8. 

where 

Y  -  a  constant  or  a  complex  number  for  the 
liquid 

T  rAPR/2L  for  the  shear  stress  for  a  straight 
tube  /  -ft* 

The  Williamson  equation  has  been  used  for  a  wide  variety  of 
pseudoplastic  liquids (36)  to  relate  either  the  shearing  stress  and 
rate  of  shear,  or  else  the  apparent  viscosity  and  rate  of  shear  in 
the  viscous  flow  region.  This  equation  was  derived  for  the  flow 
relationship  between  infinite  parallel  planes,  where  the  rate  of 
shear  is  constant.  However,  the  form  of  this  equation  is  useful  for 
laminar  flow  in  cylindrical  tubes,  and  the  constants  of  the  Williamson 
equation,  determined  from  a  small  number  of  selected  points,  are 
independent  of  the  instrument  used.  The  most  useful  form  of  this 
equation (36) 

where 

jx  -  apparent  viscosity  at  an  infinite  rate 
of  shear  lb. /ft. sec. 

a  z  a  constant;  the  value  is  defined  by  the  equation 


J  (a  /  dv/dr )1  / 
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Ty  r  intercept  on  the  shear  stress  axis  obtained 

by  extrapolating  the  tangent  to  the  consist- 

2 

ency  curve  lb.f  /  ft. 

Although  this  equation  can  only  be  used  to  determine  the 
apparent  viscosity  at  any  shear  rate  in  the  laminar  region,  the 
constants  can  be  determined  from  a  small  number  of  selected  experimen¬ 
tal  data,  whereas  the  other  methods  require  the  complete  determination 
of  the  consistency  curve. 

Alves  et  al  (1)  advocate  the  construction  of  a  plot  exhibiting 
the  shear  stress-shear  rate  relationship  for  one  liquid  from  data 
obtained  by  pilot  pipe  line  or  laboratory  viscometer.  This  graph¬ 
ical  relationship  can  then  be  used  directly  for  computations,  or  the 
curve  may  be  used  to  obtain  mathematical  relationship  for  the  shear 
stress  and  shear  rate.  By  substituting  this  mathematical  equation 
into  the  appropriate  relationship,  it  is  possible  to  obtain  a  general 
flow  equation  which  can  be  applied  directly  to  fluid  flow  problems 
within  the  limits  of  the  integration. 

Another  method  of  determing  the  equation  of  the  flow  curve 
for  non -Newtonians  has  been  outlined  by  Krieger  and  Maron(l6) . 

The  shear  stress  and  apparent  fluidity  values  are  determined  from 
capillary  flow  data  with  the  appropriate  Nexrtonian  equations.  The 
computed  values  are  then  used  in  the  following  expression. 

(dv/dr)w/  Tw~  gc  /  I A  d 

where 

^  -  apparent  fluidity  r  iAn  ft. sec. /lb. 
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The  apparent  fluidity  $  and  shear  stress  T  are  determined 

a  w 

independently  of  an  instrument  constant,  so  that  the  equation  is 
completely  general. 

By  assuming  that  Newton’s  hypotheisis  applies  to  non-Newtonians, 
it  is  possible  to  obtain  analogous  information  for  pseudoplastics 
and  a  basis  of  comparision  with  Newtonian  liquids.  In  the  laminar 
region  of  flow,  the  apparent  viscosity  of  a  pseudoplastic  liquid 
may  be  defined  by  the  rearrangement  of  the  Poiseuille  equation (3, 36) ; 

_M, a  Z  Pg0E2  /  32L? 

The  apparent  viscosity  data  obtained  on  pilot  pipe  lines  may 
be  correlated  by  means  of  the  shear  rate  at  the  wall  of  the  pipe 
determined  by  the  following  equation  derived  for  Newtonian  liquids. 

(dv/dr)w  =  -S  V/D 

where 

(dv/dr)w  z  shear  rate  at  the  wall  sec  ^ 

Since  this  equation  is  theoretically  applicable  to  Newtonian 
liquids  only,  the  shear  rate  determined  is  an  apparent  value.  The 
correlation  of  the  apparent  viscosity  with  the  shear  rate  is  consist¬ 
ent  for  one  liquid  if  the  same  consistency  curve  is  obtained  on  pipe 
lines  of  different  diameters  and  lengths.  If  different  curves  were 
obtained,  the  liquid  may  be  a  thixotropic  pseudoplastic,  for  which 
the  apparent  viscosity  is  a  function  of  the  shear  rate  history. 

A  method  of  determining  the  shear  rate  at  the  wall  of  a  pipe 
without  prior  assumption  of  a  consistency  curve  has  been  developed  by 
Rabinowitsch  (l), 
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-  (dv  /dr)  =  3 
±  w 


(dv/dr)  /  T 


(dv/d.r) 

dT 


This  equation  is  used  on  pseudoplastic  liquids (1) 

Turbulent  Flow  of  Pseudoplastic  Liquids 

There  has  been  no  equation  developed  for  the  flow  of  pseudo¬ 
plastics  in  the  turbulent  region  of  flow.  The  method  generally  used 
consists  of  determining  the  value  of  the  apparent  viscosity  by  means 
of  the  Darcy-Weisbach  equation  and  the  appropriate  friction  factor 
in  the  turbulent  region  for  several  shear  rates  (1,36).  This  value 
may  then  be  applied  to  any  rate  of  flow  in  the  turbulent  region  as 
it  has  been  determined  that  this  value  is  very  nearly  a  constant. 
iThe  value  of  the  apparent  viscosity  as  determined  at  very  high  rates 
of  shear  in  a  capillary  or  rotational  viscometer,  may  also  be  used 


(1,36). 


Effect  of  a  Temperature  Gradient  on  Newtonian  Liquids 
Laminar  Flow 

In  the  laminar  flow*a  Newtonian  viscous  liquid,  the  effect  of 
a  temperature  gradient  from  the  liquid  to  the  pipe  wall  results  in  a 


higher  friction  factor  (15) •  Because  the  heat  is  transferred  from 
the  oil  to  the  pipe  wall  by  conduction  only,  the  velocity  proiile  is 
distorted  as  a  result  of  the  change  in  viscosity  with  temperature 
as  in  Fig.  40 (10). 


FIG.  40  DISTORTION  OF  THE  PARABOLIC 

PROFI  LE 

Effect  of  a  temperature  gradient  on  a 
viscous  liquid  in  laminar  flow 
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It  may  be  shown  that  the  derivative  of  the  temperature  gradient 
at  the  pipe  wall  is  not  linear  within  a  cylindrical  tube,  but  is 
a  function  of  the  radius  as  follows  (10) 

(d2t/dy2)w  ;  1  (dt/dy)  (65) 

r 

where 

t  2  temperature  at  the  distance  y  from 
the  pipe  wall 

The  reason  the  temperature  profile  is  not  linear  is  a  result  of  the 
heat  flow  area  being  reduced  with  increasing  distance  from  the  pipe 
w  all.  For  liquids  with  small  temperature  coefficients  of  viscosity, 
the  parabolic  velocity  profile  is  only  slightly  distorted  under  a 
thermal  gradient,  and  after  a  sufficient  length  to  attain  thermal 
equilibrium,  the  temperature  profile  is  very  nearly  a  parabola. 

For  very  low  Reynolds  numbers  and  large  diameter  tubes,  it  has  been 
shown  that  free  convection  influences  the  heat  transfer  and  thus 
the  velocity  profile.  The  analysis  of  non-isothermal  fluid  has  resulted 
in  an  equation  which  can  be  used  to  determine  the  pressure  drop  in 
laminar  flow(l8) . 

This  equation  gives  the  pressure  drop  as  a  function  of  the 
length  from  the  entrance.  Martinelli  (18)  found  that  this  equation 
gave  a  favourable  correlation  with  pressure  drop  data  obtained  on 
liquids  flowing  under  large  temperature  gradients,  and  at  Reynolds 
numbers  high  in  the  laminar  region  of  flow.  Graphs  obtained  from  the 
data  of  this  study  indicated  that  the  pressure  drop  approached 
a  linear  relationship  with  distance,  after  a  calming  length  of 
approximately  200  diameters. 
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Other  pressure  drop  data,  obtained  while  determining  heat 
transfer  coefficients,  have  been  correlated  by  empirical  equations 
(15,29).  An  empirical  relationship  for  viscous  oils  being  cooled 
was  derived  from  data  obtained  on  a  section  consisting  of  100 
diameters  calming  length,  and  105  diameters  test  section,  by  Keevil 
and  McAdams (15), 

f  t/f  :  1  /  0.0133 (t0  -  tp)°G  (66) 

where 

f  -  non-isothermal  friction  factor 

o 

t  -  temperature  of  the  pipe  C 

o_ 

t  Q  -  temperature  of  the  oil  in  the  stream  G 
Another  empirical  equation  for  non-isothermal  laminar  flow 
is  presented  as  a  modified  form  of  the  Reynolds  number (32) 


where 

jjJ  z  viscosity  taken  at  temperature  s  t  - 

(t-tJ/2 

w 

t  2  temperature  of  the  wall 
t  -  average  or  bulk  temperature 
f*  -  non-isothermal  friction  factor 
From  the  experimental  data  obtained  on  viscous  oils  with  heat 
exchangers,  Seider  and  Tate (29)  computed  the  following  empirical 
relationship  for  the  laminar  flow  region, 

q/f  r  1.1  y%/jiv)  for  NRe<  2,000  (68) 

f'  -  non-isothermal  friction  factor 


where 


■ 


I 


I 


t 


I 


J- 


JXg  =  viscosity  at  the  axis  temperature 
=  viscosity  at  the  wall  temperature 

The  data  used  for  this  correlation  was  obtained  from  liquids 
with  a  wide  range  of  viscosities  which  resulted  in  a  large  spread 
of  the  data.  The  last  two  equations  were  derived  for  heating  condit' 
ions,  and  they  are  based  on  the  accurate  determination  of  the  film 
temperature.  For  these  reasons  both  equations  are  not  as  applicable 
to  this  work  as  that  of  Keevil  and  McAdams (15). 

A  theoretical  analysis  of  the  temperature  profile  of  a  liquid 

entering  a  pipe  indicates  that  the  length  required  to  obtain  thermal 

equilibrium  is  100  to  1000  times  longer  than  the  dynamic  calming 

length  in  laminar  flow(10) . 

Turbulent  Flow 

When  a  liquid  is  flowing  under  turbulent  flow  in  a  pipe  which 
is  colder  than  the  oil,  the  main  resistance  to  the  transfer  of  heat 
is  the  very  thin  laminar  sublayer.  The  heat  is  transferred  through 
this  thin  layer  by  conduction  only  with  the  result  that  the  temp¬ 
erature  gradient  is  linear  from  the  inner  edge  of  the  buffer  layer 
to  the  pipe  wall.  Because  of  the  eddy  action  and  momentum  exchange, 
the  temperature  gradient  from  the  buffer  layer  to  the  axis  is  very 
nearly  linear,  and  of  a  constant  value.  These  properties  result  in 
a  relatively  small  increase  of  the  friction  factor  for  a  temperature 
gradient  from  the  flowing  liquid  to  the  pipe  wall(l0;. 

An  analysis  of  heat  flow  in  the  turbulent  region  has  resulted 
in  an  equation  for  the  temperature  profile  as  a  function  of  the 
respective  velocities  in  each  of  the  liquid  layers,  and  the  Prandtl 
number  as  follows (10), 
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\  ~  *w  =  -^  c  °p  (  \ 

tS  "  ^  k  TVg“ 


(69) 


where 


tfc  r  temperature  of  laminar  sublayer  at  y  r  5.0 
tw  =  temperature  of  the  pipe  wall 
tg  -  temperature  of  the  flowing  liquid 

r  maximum  velocity  in  turbulent  flow  ft. /sec. 
r  velocity  of  laminar  sublayer  @  y7^-  5.0  ft. /sec. 
Pr  -  Prandtl  number.  (Cp>*/k) 


With  this  equation  the  temperature  profile  may  be  determined 


if  the  velocity  profile  is  known.  This  temperature  profile  computed 
by  this  equation,  is  illustrated  in  Fig.  (41)  with  the  effect  of 
increasing  values  of  the  Prandtl  number  shown  in  Fig.  42  (10) . 

From  this  illustration  it  may  be  seen  that  the  temperature  drop 


in  the  laminar  sublayer  is  the  essential  one  for  the  viscous  liquids 
with  large  Prandtl  numbers. 

The  decreased  effect  of  a  temperature  gradient  on  the  friction 
factor  is  attributed  to  the  greater  mixing,  decreased  effect  of  the 
viscosity  on  the  pressure  drop,  and  the  nearly  linear  temperature 
profile  in  the  turbulent  core.  Friction  factor  data  obtained  on 
viscous  oil  flowing  non-isothermally  in  pipes  have  been  correlated 
by  the  following  equation  (29), 


(70) 


where 


z  viscosity  at  the  mean  temperature 
-  viscosity  at  the  wall  temperature 


EQUATIONS  OF  LIQUID  FLOW  IN  CONCENTRIC  CYLINDER  VISCOMETERS 


Flow  of  Newtonian  Liquids 

The  geometrical  flow  pattern, obtained  when  a  bob  is 
within  a  cup  rotating  concentrically,  has  made  it  possible 
rigorous  mathematical  relationships  for  Newtonian  liquids. 

i 

system  it  may  be  shown  that  the  torque  on  the  bob  is (31), 

2  2  o 

M  -  8Trhe^N  (R  1  R2)(c) 

g=  (Rj  -  4) 


suspended 
to  derive 
For  this 


(71) 


where 


M  -  torque  lb.f  (ft.) 
c  -  dimensional  constant 
he  z  equivalent  depth  of  immersion  in. or  cm. 

Rl  r  radius  of  the  bob  ft. 

R0  r  radius  of  the  cup  ft. 

jlf  equation  (71)  is  rearranged  to  express  the  viscosity  explicitly, 
then  for  one  system  and  one  wire,  the  equation  may  be  reduced  to 

ja  ;  ^  K  /  he  N 

This  simplified  form  of  the  equation  is  used  for  the  determin¬ 
ation  of  the  wire  constants  from  the  calibration  runs  and  for  the 
determination  of  the  viscosity  or  apparent  viscosity  from  experimental 
data.  For  the  pseudoplastic  liquids,  the  apparent  viscosity  is 
correlated  with  the  shear  rate  at  the  bob  wall(l,17),  or  at  the 
midpoint  of  the  annulus  (35,31).  For  Newtonian  liquids  the  value  of 
the  shear  rate  at  the  midpoint  of  the  annulus  is  very  nearly  the  same 
as  the  average  weighted  value (31) •  This  value  is  determined  with  the 
following  equation,  by  letting  r  -  (R]_  /  Rg)  (13,35/'. 
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(dv/dr)r  =  (R^  ) 

~<4  *  R1  > 


where 


r  r  any  radius  between  and  R0 


(72) 


A  plot  of  equation  (71)  is  presented  in  Fig.  (43)  for  a  cup 
speed  of  one  revolution  per  second  on  the  Alberta  Viscometer. 

Flow  of  Bingham  Plastics 

The  flo characteristics  of  the  rotational  viscometer  make  it 
possible  to  determine  the  coefficient  of  rigidity  quickly  and  accur¬ 
ately.  When  the  rotational  speed  is  high  enough  and  the  shear  stress 
across  the  annulus  is  sufficient  to  cause  ail  the  liquid  to  flow, 
the;  linear  relationship  of  the  shear  rate  and  plastic  viscosity  may 
be  obtained (12) .  The  plastic  viscosity  is 
i/r\  -  M. 


(  12  "  h  )  ”  T,r  ln(Rc/Rb) 
~2 

2TP 


StrHThe  Rj  Ro  Y 


(73) 


The  rate  of  shear  is  determined  with  equation  (71) ,  for  the 
coefficient  of  rigidity  has  a  constant  value  if  the  shear  stress  is 

larger  than  T  .  Since  it  is  possible  to  determine  a  linear  consist- 

*y 

ency  curve  for  plastic  materials,  this  instrument  is  widely  used  for 
their  investigation  (12,34) • 

Flo\-j  of  Pseudoplastic  Liquids 

There  has  been  no  equation  derived  for  tne  flow  of  pseudoplasmc 
liquids  in  the  rotational  viscometer.  In  place  of  these,  numerous 
semi-theoretical  and  empirical  equations  have  derived  to  oetter 
interpret  experimental  data,  and  to  obtain  analogous  information 

(1,17). 
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FIG.  41  TURBULENT  FLOW 
TEMPERATURE  AND  VELOCITY 
PROFILE  WITH  LAMINAR  SUBLAYER 


FIG.  42  TEMPERATURE  PROFILES 
FOR  VARIOUS  PRANDTL  NUMBERS 
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From  the  basic  relationship  for  the  dear  rate  at  the  bob  wall, 
an  equation  has  been  derived  for  the  shear  rate  at  the  wall  with 
the  assumption  that  the  viscosity  is  an  inverse  function  of  the  shear 
rate  and  incorporating  M  =  -g/T  (dvr/dr) ,  the  equation  for  the  rate 
of  shear  at  the  bob  wall  is  (l) 

-(dv/dr)  ;  2TX  (dw/d^  )/(  1  -  (R^)  )  (74) 


where 

z  =  d(logw)/d(logT1)  (75) 

The  exponent  z  is  the  slope  of  the  log-log  plot  of  w  vs.  T]_,  at  the 
point  in  question.  This  value  of  the  shear  rate  is  used  to  determine 
the  value  of  the  shear  stress  at  the  bob  surface  from  the  general 
equation  jx>  z  -Tg  (dv/dr) . 

An  equation  has  been  derived  analytically  which  gives  a  direct 
determination  of  the  flow  curve  from  rotational  viscometer  data(l7) . 
It  is  based  on  the  assumption  that  there  is  some  functional  relation¬ 
ship  between  the  shear  rate  and  shear  stress  as  follows, 

(dv/dr)  =  ,i(T) 

=  fdw/dln  r) 


(76; 


where 


j  =  is  some  function  of 

For  a  constant  deflection  for  one  speed,  the  following  equation 
for  the  rate  of  shear  is  developed,  which  is  the  first  three  terms 
in  a  series  solution  of  eauation  (7  6)//'7;)/, 

2  2 

(dv/dr)  =  w  (  1  /  In  s *d  (lnw)_  /  (In  s)  d _ (w^  ^  (77) 

lns(  d(lnT1)  '  3v  d  (lnT^T  ) 

where 

s  “  ^2/r1 
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w  z  angular  velocity  -  2'Tr'N  Se*.-' 

T1  2  shear  stress  at  the  bob  wall  \\>.x  I- Ptl 
For  concentric  cylinder  systems  with  a  low  R2/Rl  ra'tio,  only  the 
first  two  of  the  terms  in  the  bracket  will  have  to  be  obtained  for 
an  error  of  less  than  one  percent.  This  is  a  direct  determination 
of  the  flow  curve  without  prior  assumption  of  a  flow  equation,  and 
the  equation  determined  is  independent  of  an  instrument  constant, 
and  therefore  completely  general(l7). 

The  velocity  profile  of  a  highly  pseudoplastic  liquid  flowing 
betw  een  rotating  concentric  cylinders  was  determined  by  a  hot  wire 
anemometer ( 25 ) .  This  investigation  indicated  that  the  rate  of  shear 
varied  linearly  from  the  cup  wall  to  close  to  the  bob  w  all.  Near 
the  bob  wall  the  rate  increased  very  quickly  so  that  the  value  of  the 
apparent  viscosity  as  determined  will  be  higher  than  the  actual  value. 
The  conclusion  of  this  investigation  was  that  the  determination  of 
the  shear  rate  by  any  other  method  than  the  experimental  method  as 
mentioned  above,  is  only  an  approximation. 

To  obtain  the  apparent  viscosity  directly  from  laboratory  data, 
the  assumption  has  been  made  that  the  equations  for  Newtonian  liquids 
may  be  applied  to  pseudoplastic  liquids  (3,13,22).  The  plots  of  the 
apparent  viscosity  and  the  shear  rate  are  then  used  for  the  solution 
of  flow  problems  in  pipes,  or  to  obtain  graphical  correlations  of 
rotational  data  and  capillary  data(l,22). 

The  Correlation  of  the  Concentric  Cylinder  and  Pipe  Line  Data 

There  have  been  several  methods  advocated  for  the  correlation 
of  the  consistency  data  obtained  on  capillary  and  rotational  viscometers, 
with  the  observed  data  from  the  model  and  prototype  pipe  lines. 
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Two  methods  have  been  outlined  by  Alves  et  al  (l)  which  were  used  to 
correlate  data  from  three  concentric  cylinder  viscometers  and  a 
pilot  pipe  line.  For  the  concentric  cylinder  viscometers,  the  shear 
rate  was  determined  at  the  bob  wall  with  equations  (73)  and  (74) , 
and  a  plot  was  constructed  exhibiting  the  shear  stress  as  a  function 
of  the  shear  rate.  These  plots  were  used  to  determine  the  relation¬ 
ship  of  these  variables  in  equation  form  such  that 

(dv/dr)  =  j  (Tx)  (78) 

and  this  equation  was  substituted  into  the  following  general  relation¬ 
ship  (1) 

i  (T)=  3(dv/dr)  /  T  (  d'  (dv/dr)  ) 

J  (  d  (T)  )  (79) 

In  the  above  investigation,  equation  (79)  was  also  used  to  determine 

the  shear  rate  at  the  wall  of  a  pipe  for  pseudoplastic  liquids. 

Substituting  equation  {18  into  equation  (79)  and  formally  integrating 

over  the  same  range  of  shear  rates  as  those  of  equation  {18 )  gives 

the  relationship  for  the  flow  of  this  liquid  in  a  pipe. 

The  other  method  used  by  Alves  et  al  (l)  involves  the  graphical 

integration  of  the  following  equation 

Q  =1T2  C  C  g(T  )dr.rdr  (80) 

Jo  )o  1 

These  relations  are  applicable  to  laminar  flow  only  and  in 
some  cases  tjie  range  of  shear  rate  is  limited.  For  the  pseudoplastic 
liquids  investigated  by  Alves  et  al  (l) ,  this  method  resulted  in  a 
good  correlation  of  the  shear  stress  and  rate  of  shear  data  as  obtained 
by  the  rotational,  capillary  and  pipe  line  viscometers. 
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It  is  known  that  the  velocity  profile $  for  laminar  flow  in  a 
pipe  and  in  a  rotational  viscometer,  are  not  the  same  so  that  the 
respective  rates  of  shear  at  the  pipe  and  bob  wall  are  not  a  true 
criteria  .  This  may  indicate  that  the  apparent  viscosity  as  deter¬ 
mined  by  a  rotational  viscometer  is  not  directly  applicable  to  the 
flow  in  a  conduit  (22).  Because  of  this  condition,  a  method  has  been 
developed  which  makes  it  possible  to  correlate  the  apparent  viscosity 
data  for  plastic  and  pseudoplastic  liquids.  The  technique  is  based 
on  three  assumptions (22) , 

1.  The  shear  rate  in  a  pipe  is  a  unique  function  of  V/t),  which 
is  a  simplified  form  of  the  shear  rate  at  the  wall  of  the  pipe  (equa¬ 
tion  5.) 

2.  The  shear  rate  in  a  rotational  viscometer  is  a  unique 
function  of  the  rotational  speed  w, .  This  is  a  simplified  form  of 
equation  (17) 

3.  The  functional  forms  of  both  the  relationship  stated  are 
identical,  or 


(81) 


shear  rate  in  the  pipe  r  j  (V/wD) 
shear  rate  in  the  vise. 


This  method  has  been  used  to  correlate  the  pressure  drop 
calculated  from  rotational  viscometer  data  with  the  observed  pressure 
drop  data  in  a  conduit.  To  correlate  the  apparent  viscosity  by  this 
method,  it  is  necessary  to  obtain  data  on  both  a  rotational  viscometer 
and  on  a  conduit  for  the  same  liquid.  This  data  is  then  correlated 
as  follows 


(82) 


(Mp)  "  (Dw) 


where 


Z  apparent  viscosity  in  the  rotational  viscometer 
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/*p  r  apparent  viscosity  as  obtained  on  a  conduit 
It  is  onlj  necessary  to  obtain  a  small  number  of  points  with 
the  rotational  viscometer,  but  it  is  necessary  to  obtain  a  large 
number  of  points  on  the  conduit  to  fully  investigate  the  laminar 
region  of  flow.  The  correlation  of  the  data  by  equation  (81)  is 
best  determined  graphically  as  in  Fig.  (44) •  This  plot  can  then  be 
used  to  determine  the  pressure  drop  in  any  size  of  conduit. 

FIG.  44 

Graphical  Correlation  of  the 
Apparent  Viscosity  Obtained 
on  a  Rotational  Viscometer 
and  a  conduit 

(  V  /  D  w  ) 

EFFECT  OF  PRESSURE  ON  THE  VALUE  OF  VISCOSITY 

In  the  derivation  of  the  Darcy- Wei sbach  equation  of  liquid  flow 
it  is  assumed  that  the  viscosity  throughout  the  full  length  of  the 

|| 

section  is  a  constant  value.  Where  there  is  a  high  average  pressure 
in  the  section,  or  a  large  pressure  drop,  the  effect  of  the  pressure 
on  the  viscosity  of  crude  oils  is  of  sufficient  order  to  be  measured(2). 

The  effect  of  pressure  on  the  viscosity  of  organic  liquids  has 
been  investigated  by  Warber  and  Sacks,  Pause,  and  Bridgman (13) .  Their 
work  showed  that  the  viscosity  increased  linearly  with  pressure  up  to 
1,500  atmospheres  pressure,  for  various  organic  compounds.  Similarly 
for  crude  oils,  the  effect  of  pressure  on  the  viscosity  has  been  found 
to  be  linear  up  to  1,000  psi(2).  These  investigations  also  show  that 
the  effect  of  pressure  on  viscosity  is  larger  when  the  liquids  are 
colder,  and  also  for  the  lower  specific  gravity  crude  oils  than  for 
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the  more  dense  oils.  Therefore  the  effect  of  pressure  on  viscosity- 
may  be  determined  by  comparing  the  values  of  the  viscosity  at  the 
average  pressure  in  a  conduit  at  the  low  and  high  pressures.  The 
following  equation  is  being  proposed  for  the  determination  of  the 
pressure  effect  on  viscosity 


C 


(8?) 


where 


(  /*lp  )  (?hp-Plp)PSi« 


0 


coefficient  of  pressure  effect  %/lOOpsi 


/*lp 


viscosity  at  the  low  average  pressure  lb./ft.sec. 


viscosity  at  the  high  average  pressure  lb./ft.sec. 


P 


lp 


average  low  pressure  in  psi. 


average  high  pressure  in  psi. 


'  J  '  '  !  '  '  '•  r  -j  or  r-x^oo-e,  ,  ... 


'  i,  .  ,t 
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APPENDIX  B 

EQUATIONS  USED  TO  COMPUTE  THE  PILOT  PIPE  LINE  DATA 
(For  nomenclature  see  Pgfc«) 

This  section  is  concerned  with  the  derivation  of  the 
formulae  used  on  the  5/8"  test  line,  and  the  use  of  these 
equations  on  one  test.  The  same  procedure  is  used  for  the 
l/2"  and  3/8"  test  lines  data,  and  the  formulae  for  all  the 
test  lines  are  summarized  in  Table  1-B, 

Development  of  the  Equations 
Average  Oil  Temperature 


FIG.  45  LOCATION  OF  PRESSURE  TAPS  AND  THERMOMETERS 
ON  THE  PILOT  PIPE  LINE 

The  inlet  and  outlet  oil  temperatures  are  obtained  at 
Tiand  T2  respectively,  as  illustrated  in  Fig.  (45).  Since 
the  oil  flow  is  perpendicular  to,  and  also  along  the  thermo** 
meters,  it  is  assumed  that  the  bulk  oil  temperature  is  obtained. 

At  the  present  time  there  is  no  method  of  determining  the 
temperature  gradient  curve  between  T^  and  T2»  so  it  is  assumed 
that  it  is  linear.  Referring  to  Fig.  (45),  the  average  oil 


temperature  is 
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^ave  -  -+’  Tc)/2  (l«*b) 

where 

Tt)  s  bulk  temperature  at  P^,  °F 
Tc  -  bulk  temperature  at  P2,  °F 
Because  the  oil  temperature  is  not  measured  at  points 
b  and  c,  the  average  temperature  is  determined  by  the  ratio 
of  lengths  between  the  two  points  Tj  and  T2. 

For  .>  T2 


Tb  +  To 

s  Ti  - 

(0.5  +  8.38  +  9.8 )  (T]_  -  T2) 

2 

29.63 

Tave 

» 

rH 

EH 

Si 

0*629  (Ti  -  T2) 

(2~b) 

For  T2  >  Tq 

Tave 

S  T1  + 

0.629  (T2  *  Ti) 

(3-b) 

To  obtain  the 

corrected  downstream  temperature, 

a  call 

brated  chart  was  used  fcr  T2  on  the  Saskatchewan  blend  runs. 
Recalibration  of  this  thermometer  indicated  a  general 
correction  of  0.4  °F,  with  a  deviation  of  approximately  0.2  °F. 
For  the  remaining  crudes  the  correction  of  0.4  °F  was 
introduced  into  the  equations  (2®b)  and  (3-b),  and  the 
following  relationship  was  used. 


For  TX  >  T2  Tave 

-  Tx  -  0*63  AT  + 

o 

• 

ro 

(4-b) 

For  TX <  T2  Tave 

=  T1  t  0.63  AT  + 

0.2 

(5-b) 

where 

1ST 

«  T]_  -  T2  or  T'2 

"  Tx 

Bulk  Velocity 

The  bulk  velocity  is  determined  as  follows  from  the 
measured  rate  of  flow  Q, 


V 
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=  Q/A 

where 

V  m  bulk  velocity  ft. /sec. 

Q  s  flow  rate  ft.Vmin.,  or  cfm. 

A  a  area  of  internal  cross  section  ft.2 
substituting  in  the  area  of  the  5/8"  test  line 

V  *  £  ft»Vmlh.  mln./60  secx _ 

'rf  (0, 386 )2ft .2/  4  x  144  in. 2 

V/  m  Q,  x  102/4.87  ft, /sec.  ( 6«*b ) 

Friction  Factor 

Rearranging,  the  Darcy«*Weisbach  equation,  and  expressing 
the  friction  factor  explicitly 

f  s  2  Dgc  &p/LV2p 

For  one  run,  the  variables  are  the  pressure  drop  and  the 
bulk  velocity ,  Then  for  the  5/8's  test  line,  the  friction 
factor  is 

2  x  0,386  in.  x  32.2(lb)(ft  _)  AP  Ibf  12  in 
f  _  _  ( llqg)  sec2)  Tn2  ft 

P  “7 3  x  19.5  ft.  v  — * 

ft-*  sec 

f  ■  15.32  ap/v2p  (7-b.) 

Laminar  Flow 

Reynolds  Number 

Nrq  g  64  /  f  (8»b<) 

Apparent  Viscosity 

Rearranging  the  Hagen«Poiseuille  equation,  the  apparent 
viscosity  is  obtained  directly, 
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nA  -  &c .  g2 

32.2  (lb  )(ft  )  0.3862  in2  AP  lbf 

(lbf)(sec2)  in2 

32  x  19.5  ft.  x~V  ft. 

sec . 

-^a  =  7.71  x  10-3  £p/v  lb/ft.  sec.  (9-b) 

Rate  of  Shear 

From  equation  (5)  the  shear  rate  at  the  wall  in  laminar 
flow  is 


(dv/dr)w  «  ~8V  /  D 

We  are  not  interested  in  the  algebraic  value  of  the  shear 
rate  at  the  wall.  It  may  be  computed  more  directly  from  the 
volume  rate  of  flow  as  follows, 


V  s 
(dv/dr )w  * 


Q4  ArD2 

4x84  ^ 

_ _ min  60  sec 

(0.386  in)3  ft3 

«(  12  Tn") 


(dv/dr )w  e 


5.10  x  103  Q,/sec . 


( 10-b ) 


Turbulent  Flow 

Reynolds  Number 

From  the  Heltzel  relationship,  the  Reynolds  number  is 
computed  from  the  friction  factor  as  follows, 

NRe  =  (0.364/f)  3-77  (H-b) 

Apparent  Viscosity 

Knowing  the  Reynolds  number,  the  viscosity  is  obtained 


with  the  following  equation, 
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jua  z  VDp/  NRe 

V  ft  0,386  In  ft  p  16 

_  sec _ 12  in.  ft 3 

NRe 

~  0,0322  (V  y  /  NRe)  lb/ft.  sec,  (12-b) 

Rate  of  Shear 

The  relationship  of  the  shear  stress  at  the  wall  of  the 
pipe  to  the  pressure  drop  is  (Equation  31,  Appendix  A), 

Tw  -  AP  D  /  4  L  (13«b) 

and  the  relationship  of  the  shear  stress  and  the  rate  of 
shear  is  (Equation  4), 

Tw  s  -  jm*(dv/dr)^c  (14-b) 

Substituting  Equation  (l4«b)  into  Equation  (13~b),  the 
relationship  for  the  rate  of  shear  is, 

(dv/dr)w  =  D  gc  /  4  L 1JLA 

A  P  lbf  x  Q*?86  in  32,2  lb  ft 

in^  Ibf  sec^ 

8  - - - — - - - - - - - - - - - 

4  x  19-5  ft-  ftrfe 

(dv/dr)w  .  1.912  (&P/^  )  sec?1  (15-b) 

Equations  Developed  for  the  Effect  of  Pressure 
on  the  Apparent  Viscosity 

To  determine  the  effect  of  a  high  average  pressure  on 
the  value  of  the  apparent  viscosity,  it  was  necessary  to 
determine  a  smoothed  curve  for  low  pressures,  and  then 
determine  the  increase  in  the  consistency  at  higher  pressures. 
The  smoothed  curve  for  the  low  pressure  data  was  checked 
between  the  high  pressure  tests  to  detect  any  change  as  the 
result  of  weathering.  The  pressure  effect  coefficient  "C" 
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was  computed  as  follows  for  each  high  pressure  point, 

0  z  tyh  “  100  x  100'J  /  Fx(ph  “  ;:V  ( 16— b) 

where 


V-h  -  apparent  viscosity  at  high  average  pressure 
Jdl  -  apparent  viscosity  value  at  the  low  pressure 
Ph  —  average  high  pressure  psi. 
p]_  z  average  low  pressure  psi. 

C  -  pressure  effect  coefficient  in  %  j  100  psi. 

The  pressure  effect  coefficient  was  averaged  for  each  high 
pressure  point,  and  the  values  were  averaged  for  each  10°F  range. 

The  downstream  pressure  was  controlled  by  a  throttle  valve 
approximately  one  foot  from  the  thermometer.  The  large  pressure  drop 
at  the  downstream  throttling  valve  may  result  in  an  increase  in  the 
oil  temperature.  If  it  is  assumed  that  the  oil  is  being  throttled 
adiabatically ,  the  general  energy  equation,  neglecting  potential  and 
kinetic  nergy  changes ,  is 

(pl  -  P2)  /  P  :  (17-b) 

The  general  equation  for  the  change  in  enthalpy  of  a  liquid  is 


dH  -  G  dT  /  (  c>  h)  dP 

"  p  (>P)T 

and  for  an  isothermal  process  the  expression 
(*H)  -  v  -  T  yv\ 

UP)  (Tf)p 


(18-b) 

for  the  differential  is 
(19-b) 


Combining  equations  (18-b)  and  (19-b)  for  an  adiabatic  process,  and 
assuming  that  the  specific  heat  is  a  constant  for  a  small  temperature 
rise,  the  integrated  expression  is 

c  (T  -  T  )  -  -  v(p  -  P  )  /  T(^v)dP  /  (2>T)  (20-b) 

jp  _i_  J— 

The  last  term  may  be  determined  by  taking  ohe  reciprocal  slope  for 
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the  respective  crude  in  Fig.  50.  For  the  Combined  blend, 

.  -> 

i  ^  v)  -  6.12  x  10  °  f»r  /  lb.°F 

rs  bp 

Evaluating  this  term  at  60  F  (520  R) ,  and  assuming  a  pressure 

drop  of  500  psi.,  a  specific  heat  of  0.45  BTU  f  lb.°F,  equation  (20-b)  is 

0.5(T  -  T 2)-  -  0.0175  ft.3(0  -  500)  x  144  lb  r  - 

ft.2 

°  -6  3 . 

520  R  x  b.12  x  10 _ ft.  (0  -  500)  x  144  lb.f 

lb.°R  ft.- 


(q  -  T2)  -  2.9°F 


The  actual  temperature  rise  may  be  smaller  than  the  adiabatic 
value  of  2.9°F  as  a  result  of  heat  conduction  to  the  cold  bath.  Because 
of  the  wiae  variation  in  the  flow  rates,  and  because  the  temperature  is 
only  determined  within  £_  2°F,  it  is  assumed  that  a  temperature  correction 
to  the  downstream  side  is  not  justified  for  the  high  pressure  runs. 


Sample  Computations 

Consider  a  run  on  the  Pembina  crude  in  which  the  flow 
condition  may  be  laminar  or  turbulent.  The  observed  data  as  obtained 


from  Test  No.  E-7  in  Table  6-E  is 

Pipe  size  5/Bh 

Pressure  drop  (152.7  -  11. B)  140.9 

Average  upstream  temperature  40. 2° F 

Average  downstream  temperature  38*5  F 

Oil  temperature  in  receiver  42.0°F 

3 

Volume  collected  5.15  it. 

Time  of  collection  z>*3~  mins. 


From  equation  (1-b)  the  average  temperature  of  the  teat  is 


Ti  -  (0.629  AT) +0.2 

40.2  -  (0.629) (40.2  -  38.5) 

40.2-1.1  +0.2 

39.3  °F 

The  density  of  the  crude  at  the  measuring  tank  temperature 

is  53.1  lb/ft3. 

The  density  of  the  crude  in  the  receiver  is  53.0  lb/ft3. 

The  rate  of  flow  Q,  corrected  to  pipe  conditions  is 

.  52*°  x  mis  ft3 

Q  3  53.1  3.32  min 

=5  1.55  ft3/min 

Sample  Computations 

The  bulk  velocity  from  equation  (6«b)  is 
V  »  Qx|C2/4.87  ft/sec 

«  31.8  ft/sec 

The  term  &P/V  appears  in  both  Equation  (7-b)  and 
Equation  (9-b),  so  this  was  computed  and  recorded, 

«  140.9  ^5  7  31.3  — 

=  4.43  lb. ft.  /  in2  sec. 

The  laminar  apparent  viscosity  is  then  computed  directly  as 
in  Equation  (9-b), 

A>i-a  =  7.71  x  10- 3  AP/V 

s  0.0342  lb/ft.  sec. 

The  turbulent  friction  factor  is  then  computed  as  in 
Equation  (7-b), 


^ave 


Tave  » 


V 


f 


f  a 

The  Reynolds  number  in  the  laminar  region  is 
NRe  ss  64  /  f 

s  64  /  0*0402 

%e  «  1»  550 

It  is  not  possible  to  determine  whether  this  is  laminar 
or  turbulent.  The  laminar  shear  rate  for  this  flow  rate  is 
obtained  by  Equation  (lo-b). 

(dv/dr  )w  s  5,10  x  105  q  /  Sec. 

s  5.1 O  x  103  x  1*55  /  sec. 

(dv/dr)w  a  7,900  secT^ 

The  Reynolds  number  is  now  computed  for  turbulent  flow  by 
Equation  (  P  ~b)  . 

NRe  s  (0.364  /  f)3s7F 

s  (0,364  /  0.0402  )V 77 
%e  *  4,010 

From  Equation  (lZ~  b)  the  apparent  viscosity  is 

>La  s  0.0322  (  Vp  /  NRe) 

S  0.0322  (31.-8  x  53.1  /  4,03.0) 

5  0.0136  Ib/ft , sec. 

The  shear  rate  at  the  wall  for  turbulent  flow  by 
Equation  (/5~-’b), 


(AP)  15.32 
(V  )  Vf 

4.43  x  15.32  /  3.18  x  53.1 
0.0402 
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(dv/dr )w  3  1.912  (  AP/^u  )  sect1 

2  1.912  (140.9  /  0.0136) 

s  19,700  sec*:1 

The  Reynolds  number,  and  the  value  of  the  apparent 
viscosity  indicate  that  the  flow  may  be  turbulent. 
Effect  of  Pressure  on  the  Apparent  Viscosity. 

Consider  Test  No.  C«63  for  the  Combined  Blend. 

Inlet  pressure  739.9  psi. 

Exit  pressure  58  2  psi. 


Average  oil  temperature 


86.2 


Apparent  Viscosity  O.O565  Ib/ft.sec. 
at  average  pressure  of  661  psi. 

Apparent  viscosity  at  10  5. .0  psi.,  at  86.2  °F  as  obtained 

from  Fig.  20,  is  0.043  lb/ft.  sec.  Then  by  Equation  (Mo~b) 

the  pressure  effect  coefficient  is 

( 0.0565  »  0.043 )  (  10 ^  %  )  lb, /ft  .sec. 

G  S  (  0.043  )  (661  «  105.0)psi.lb./ft .sec 

c  s  5.65  %  /  100  psi. 
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TABLE  1-3 

PILOT  PIPE  LINE  EQUATIONS 


Nominal  Diameter, in. 

3/8" 

1/2” 

5/8" 

Internal  Pipe  Diameter, in 

.  0.249 

0.311 

0.386 

Gross  Section  Area,  ft2 

338xl0“4 

5.28xl0~4 

8.12xl0"4 

Velocity,  ft/sec 

Qx  IQ2 

QxlO2 

QxlO2 

5787 

2.03 

3,17 

Friction  Factor 

9.88AP 

V2P 

12.35AP 

V2P 

15.32AP 

V2p 

Reyno Id  s  Number , laminar 

turbulent 

64/f 

(0.364/f )3. 77 

64/f 

(0.364/f ) 3. 77 

64/f 

(0.364/f )3. 77 

Apparent  Viscosity , lb/ft . 

sec. 

laminar 

turbulent 

3. 21x10”  ^AP) 

(  v) 

0.0208 (Vd) 

4.99xl0~3(AP) 

Cv) 

0.0259(Vp) 

7.71xlO-3(AP) 
(  V) 

0.0322  (Vo) 
<%e> 

Shear  Rate  at  Wall,  sec7^ 

laminar 

18.92xl05Q 

9.74x10-^0, 

5.10x105Q 

turbulent 

1.231(AP/m) 

1.541  (AP//J 

1  •  912  (4P^u  ) 

where 

Q,  a  flow  rate  ft?/min* 

A P  e  pressure  drop  lbf/in? 

p  »  density  lb/ft? 

a  bulk  velocity  ft. /sec. 


V 


:r 


■i- 
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TABLE  2  •»  B  NOMENCLATURE 

G 

pressure  effect  on  viscosity  jg/100  psl. 

Cp 

specific  heat  BTU/lb.°F 

f 

friction  factor 

D 

pipe  diameter  ft. 

6c 

acceleration  of  gravity  (lb. /lb. f ) (ft ./sec.2) 

L 

length  of  test  section  ft. 

%e 

Reynolds  number 

P 

density  lb. /ft. 3 

P 

pressure  at  a  cross  section  psi. 

Subscripts  L  and  P  indicate  high  and  low  average 
pressures . 

AP 

pressure  drop  psi. 

Q 

volume  rate  of  flow  ft.^/niin. 

T 

J’W 

shear  stress  at  wall  lb.f/in.2 

T 

xave 

average  oil  temperature  in  the  section 

Ti,  To 

temperature  at  the  upstream  and  downstream  sections 
respectively,  °F 

Tb »  Tc 

temperature  at  the  inlet  and  outlet  of  the  test 
section  respectively 

V 

bulk  velocity  ft. /sec. 

A* 

apparent  viscosity  lb. /ft.  sec. 

At 

apparent  viscosity  at  high  average,  and  low  average 
pressures  respectively 

(dv/dr)w  shear  rate  at  the  wall  sec. 
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APPENDIX  C 


CALIBRATION  OF  THE  TORSION  WIRES  AND  TORSION  WIRE  CONSTANTS 

TEMPERATURE  SURVEY  OF  THE  ALBERTA  VISCOMETER 


SAMPLE  COMPUTATIONS  FOR  THE  ALBERTA  AND  MACMICHAEL  VISCOMETERS 
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APPENDIX  C 

(For  nomclature  see  page  C -/$) 

Calibration  of  the  Torsion  Wires 

All  the  torsion  wires  used  on  the  MacMichael  and  Alberta 
instruments  were  calibrated  by  the  procedure  outlined  by 
Villiers  Fisher  (31) •  Each  wire  used  was  calibrated  with 
at  least  two  standardized  Newtonian  oils,  at  several  depths 
of  immersion.  These  oils  were  recalibrated  recently,  and  the 
changes  in  the  viscosity  were  less  than  one  percent. 

By  plotting  the  deflection  for  various  speeds  and  depths 
of  immersion,  the  torsion  wire  constants,  and  also  the  end 
effect,  were  computed  (12),  The  individual  constants  and 
the  average  weighted  constants  are  presented  in  Table  (1**C) 
and  (2~C),  The  end  effect  of  the  MacMichael  was  found  to  be 
0,24  T  0.04  cm,,  as  compared  with  0,25  cm,  determined 
previously  (30).  The  end  effect  on  the  Alberta  viscometer 
was  determined  to  be  0.16 "  —  0.03,  as  compared  with  0.18 
determined  by  Villier3  Fisher  (31) « 

The  correlation  of  the  calibration  constants  may  be 
determined  by  plotting  the  wire  constants  as  a  function  of 
the  wire  diameter  on  log»log  paper.  For  the  Alberta  and 
MacMichael  viscometers  the  slope  of  this  plot  is  0,248  and 
0,252  respectively,  as  compared  to  the  theoretical  value 

of  0.250. 

The  maximum  deviation  on  the  Alberta  torsion  constants 
was  5,9%  for  wire  no.  24»1,  as  a  result  of  its  being  bent, 
and  the  average  deviation  was  2,7%, 
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For  the  MacMichael  torsion  wires,  the  maximum  deviation 
is  with  an  average  deviation  of  1.4$. 

The  temperature  of  the  oil  was  maintained  at  70°FfOa2°F, 
with  a  room  variance  of  70°F  to  80°F.  There  was  no  method 
of  determining  the  average  oil  temperature  when  the  wire 
constants  were  being  determined,  so  the  size  of  this  error 
may  be  as  large  as  0C3°F,  as  a  result  of  the  temperature 
gradients . 

Because  the  torsion  wire  constants  are  dependent  on 
the  diameter,  length,  and  geometry,  of  the  bob  and  stator 
system,  it  is  possible  to  determine  the  theoretical 
correlation  between  the  MacMichael  and  Alberta  viscometers. 

The  ratio  of  the  torsion  wire  constants  may  be  determined 
with  the  following  expression, 

%  a  °M  (av/dr)RlA 
kM  °A  ( d  v/dr )  Rljj  R^m 

Substituting  in  the  expressions  for  the  shear  rates  at 
the  bob  surfaces  of  the  Alberta  and  MacMichael  viscometer,  the 
ratio  and  the  respective  bob  radii,  the  ratio  becomes 

Ka  /  K„  ,  16.6 

This  is  based  on  the  assumption  that  the  Alberta  wire  has 
a  freely  suspended  length  of  10,5  inches,  whereas,  it  is  less 

than  this,  so  that  the  actual  ratio  may  be  as  low  as  loa3.  Fox 

the  wires  which  may  be  correlated  as  above,  the  ratios  are 

16.2,  16.3,  and  16.7,  indicating  a  maximum  deviation  of  2$ 

between  the  instrument  constants  of  the  two  instruments. 
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WIRE  CONSTANT  ALBERTA  ( in.cps  /  °A  sec.  )  MACMICHAEL  (cm.cps./  °M  sec.) 
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TABLE  1  *  C 


TORSION  WIRE  CONSTANTS 

MacMIchael  Viscometer 


Gauge 

Wire 

No. 

&  Wire 
Diameter 
in. 

Standard 
Oil 
cps . 

Depth 

1.6cm 

of  Immersion 
2.6cm  3.6cm 

Weighted  Averaged 
Constant  Value 

34-1 

0.0062 

310.3 

1.19 

1.11 

1.17 

1.15 

205,5 

1.13 

1.15 

1.16 

1.15 

1.15 

34-2 

0.00615 

73,0 

1.085 

1.085 

1.124 

205.5 

1,144 

1.144 

32-4 

0.00787 

310.3 

2.95 

3*06 

3.10 

3.05 

205.5 

3.00 

3.04 

3.04 

3.03 

3.04 

32-2 

0,0079 

475.5 

3.08 

3.06 

3.13 

3.10 

310.3 

3.32 

3.12 

3.12 

3.11 

30-1 

0,0101 

475.5 

8,06 

7.8 

7,95 

7.92 

7.92 

28-1 

0.0119 

475.5 

15.6 

15,4 

15.6 

15.5 

310.3 

16.0 

15.3 

15.24 

15.4 

15.4 

End  Effect  m  *0.24  cm. 

Effective  Depth  hg  c  (h  4-  0.24)cm. 

*Units  of  the  wire  constants  are  (cm.)(N  rev  ./sec..JJ_c2s_. ) 

°M 


' 
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TABLE  2  -  C 
TORSION  WIRE  CONSTANTS 


Alberta  Viscometer 
Gauge  &  Wire  Standard 


Wire 

Diameter 

Oil 

Depth  of  Immersion 

Weighted 

^Averaged 

No. 

in. 

cps. 

1” 

2" 

3" 

4" 

Constant 

Value 

32-3 

0.00798 

10.92 

0.191 

0.187 

0.189 

0.190 

0.185 

5.477 

0.177 

0.184 

0.184 

0.189 

0.189 

0.187 

30-3 

0.010 

20.70 

0.490 

0.493 

0.496 

0.494 

10  ..92 

0.434  0.482 

0.489 

0.485 

0.485 

0.490 

26-1 

0.016 

205.5 

3.29 

3.28 

3.23 

3.26 

3.26 

310.3 

3.35 

3.26 

3.22 

3.18 

3.23 

3.25 

26—2 

0*0161 

20.70 

3.04 

3.12 

3.18 

3.19 

3.17 

73.51 

3.28 

3.23 

3.21 

3.24 

3.23 

205.-5 

3.18 

3.15 

3.13 

3.23 

3.18 

3»  19 

24-1 

0.0198 

475.5 

6.98 

7.15 

7.08 

310.3 

7.27' 

7.25 

7.26 

7.17 

22-1 

0..025 

475.5 

19.0 

19.5  : 

19.5  19.8  : 

19.5 

19.5 

20-1 

0.0321 

475.5 

54.4  53.8  ! 

54.1 

54.1 

End  Effect  -  -0. 

16 

Effective  Depth  - 

•  h  +  0 

...16 

*  Units 

of  the  wire  constant  ; 

are  (in.)(N 

rev  ,/i 

3ec .  )  (cps 

.) 

°A 
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Temperature  Survey  of  the  Alberta  Viscometer 

Preliminary  tests  were  conducted  with  the  Alberta 
rotational  viscometer  on  a  crude  which  was  previously  tested 
on  the  pilot  pipe  line,  over  a  wide  range  of  temperatures. 

These  tests  indicated  that  the  oil  temperature  was  different 
than  the  circulating  temperature,  and  that  the  difference 
increased  with  lower  circulating  temperatures. 

To  determine  the  average  temperature  of  the  oil,  a 
temperature  survey,  of  the  annulus  under  varying  temperature 
and  rotational  speeds  was  undertaken.  Thermocouples  were 
attached  to  the  bob  surface  in  such  a  way  as  to  measure  the 
surface  temperatures  of  the  bob  and  cup  at  the  critical 
points  for  various  depths  of  immersion,  as  illustrated  in 
Fig,  (48).  This  data  is  presented  in  Table  (3**C)»  and  one 
run  is  plotted  in  Fig.  (48),  Experimental  runs  were  made  at 
three  depths  of  immersion  for  the  Pembina  crude  and  the 
Saskatchewan  blend.  From  Fig.  (48)  it  may  be  seen  that  the 
temperature  gradient  within  two  inches  of  the  oil  surface,  is 
very  steep,  so  that  the  runs  at  this  depth  of  immersion  would 
be  very  susceptible  to  error. 

Since  the  flow  of  the  oil  in  the  annulus  was  laminar, 
it  was  assumed  that  the  heat  flow  was  a  result  of  conduction 
only,  and  that  any  heat  transferred  by  convection  currents 
was  negligible.  The  temperature  gradient  from  the  bob  wall 
to  the  cup  wall  was  assumed  to  be  linear  for  the  small 
distance  involved.  Then  by  plotting  the  temperature  profile 
of  the  bob  and  cup  surface,  Fig. (48),  it  was  possible  to  deter* 
mine  the  average  temperature  of  the  respective  surface  by 
integration. 
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FIG.  47 

TEMPERATURE  SURVEY  OF  THE  ALBERTA 
VISCOMETER 

C.F.  CIRCULATING  FLUID 
H.F.  DIRECTION  OF  HEAT  FLOW 

POSITION  OF  THERMOCOUPLES 

2  END  OF  BOB  1/8"  INTO  OIL 

3  2" FROM  END  l/8"lNT0  OIL 

4  3"FROM  END  1/8"  INTO  OIL 

6  3  3/4" FROM  END  , ON  BOB 

7  2  1/4"  FROM  END,  ON  BOB 
9  4  1/4"  FROM  END  ,  ON  BOB 
B  INSIDE  BOB,  ON  BOTTOM 


The  average  temperature  of  the  oil  from  the  survey 
plots  is 

Tave  =  R1  Tc  +  R£  Tb 
Rl  +  R2 

where 

Tc  «  average  weighted  temperature  of  cup  wall  °F 

s  average  weighted  temperature  of  bob  wall  °F 

Tave  s  average  weighted  oil  temperature  in  the  annulus 

The  average  weighted  oil  temperature  as  determined  above 
was  found  to  be  slightly  higher  than  the  temperature  obtained 
with  the  thermocouple  located  in  the  bottom  of  the  bob0 

The  correction  to  the  thermocouple  temperature  in  the 
bob,  or  the  bob  temperature,  was  correlated  with  the 
difference  between  the  average  circulating  temperature  and 
the  ambient  temperature.  This  plot  (Fig. 49),  indicates  that 
the  correction  to  the  bob  temperature  for  three  inch  and 
four  inch  immersion,  was  linear,  whereas  a  curve  was 
obtained  for  the  two  inch  immersion.  In  all  the  tests 
conducted,  only  three  and  four  inch  depths  of  immersion  were 
used,,  because  of  the  high  thermal  gradients  obtained.  During 
the  test  runs,  the  bob  temperature  was  used  to  indicate 
thermal  equilibrium,  and  also  as  a  method  of  empirically 
determining  the  average  oil  temperature. 
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TABLE  3  -  C 

Alberta  Viscometer  Temperature  Survey 
Observations 


4M  Immersion 
Thermocouple 

Pembina 

Crude 

Circul. 

Room 

2  3  4 

9  7 

6 

B 

*1 

Op 

t2 

Op 

Fluid 

Op 

Temp. 

Op 

34.5 

54.5 

58.0 

62.5 

56.0 

60.5 

55.0 

52.4 

53.6 

53.0 

70 

42.0 

43.5 

46.0 

54.5 

45.0 

52.5 

43.5 

43.8 

44.2 

44.0 

75 

34.5 

35.5 

39.5 

52.5 

38.5 

48.0 

36.5 

34.5 

34.9 

34.7 

75 

No  difference  : 

In  temperature  readings  from  0, 

,5  to 

1.0  R 

21.0 

24.5 

27 

41 

25.5 

37.5 

22.5 

22.3 

23.3 

22.8 

72 

3"  Immersion 

Pembina  Crude 

42.0 

46.0 

50.5 

56.5 

47.0 

55.5 

44 

40.2 

40.6 

40.4 

75 

48.0 

49.5 

55.0 

59.5 

52.0 

61.5 

50.5 

46.2 

46.6 

46.4 

75 

55.0 

56.5 

60 

58.5 

57.0 

54.4 

54.4 

54.4 

67 

61.5 

61.5 

62.5 

62 

63 

62 

60.2 

60.6 

60.4 

65 

27.0 

31.5 

36.0 

34.5 

29.0 

26.3 

27.3 

26.8 

72 

22.5 

28 

35.5 

30.5 

26.0 

22.3 

23.3 

22.8 

72 

TABLE  4  «  G 


Alberta  Viscometer  Temperature  Survey 
Observations 

3"  Immersion  Saskatchewan  Blend 
Thermocouple 


2 

3 

4 

7 

6 

24.5 

31.0 

40.5 

35 

35 

38.0 

51.0 

41.0 

42.5 

46.5 

52.5 

48.0 

48.5 

51.5 

55.5 

52.5 

61.0 

50.5 

54.0 

57.5 

54.5 

62.5 

2"  Immersion 

22 

35.5 

43 

41 

35.5 

47.0 

57.5 

50 

43.5 

52.0 

61.0 

55 

49.0 

56.5 

63.0 

58.0 

62.0 

64.0 

63.2 

64 

B 

Ti 

T2 

T 

ave 

Room 

Temp 

op 

op 

°p 

Op 

28 . 

23.3 

25.1 

24.2 

71.0 

38.0 

31.7 

32.5 

31.8 

70 

46.5 

38.5 

39.4 

39.0 

70 

50.0 

44.8 

46.6 

45.7 

70 

52.0 

47.4 

47.-8 

47.6 

70 

30 

23.0 

23.0 

71.5 

39 

30.8 

31.8 

31.3 

74 

47.0 

39.9 

41.0 

40.4 

70 

52.0 

47.8 

47.8 

47.8 

72 

63.0 

62.6 

62.6 

62.6 

72 

TABLE  5  -  G 
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Alberta  Viscometer  Temperature  Survey 
Computations 


Bob 

Integrated  Correction 

Circulation 

Room 

Circ  e^Room 

Temp. 

Oil  Temp*  to  Bob  Temp* 

Temp. 

Temp. 

Temp, 

Op 

Op 

Op 

Op 

Op 

Op 

4,!  Immersion 

Pembina 

55 

56.5 

1*5 

53 

70 

17 

43.5 

45.3 

1.8 

44 

75 

31 

36.5 

38.7 

2.2 

30*5 

75 

34.5 

22.5 

26.3 

3.8 

22.8 

72 

49 

3"  Immersion 

Pembina 

62.0 

62*0 

0 

60*6 

65 

5.4 

57  *0 

57*4 

0.4 

54*4 

67 

12*6 

50*0 

50*7 

0*7 

46.4 

75 

28*6 

44.0 

46,0 

2.0 

40*4 

75 

34*6 

29*0 

31*7 

2*7 

26*8 

72 

45.2 

26*0 

27*9 

2.9 

22*8 

72 

49.2 

Alberta  Viscometer  Temperature 

Survey 

Computations 

Ave. 

Oil 

Bob 

Correction 

Cabinet 

Circuit 

Circ.wRoom 

Temp* 

Temp. 

to  Bob  Temp. 

Temp. 

Temp. 

Temp* 

Op 

Op 

Op 

Op 

Op 

Op 

3"  Immersion 

Saskatchewan  Blend 

53.8 

52*0 

1*8 

70 

47*6 

22.4 

51.2 

50.0 

1*2 

70 

45.7 

24.3 

47.1 

46*5 

0*6 

70 

39.0 

31.0 

39.4 

38*0 

1.4 

70 

31*8 

38*2 

31.0 

28*0 

3.0 

71 

24.2 

46.8 

2"  Immersion 

63.3 

63 

0*3 

72 

62*6 

9.4 

53*2 

52 

1.2 

72 

47.8 

24*2 

48*4 

47 

1.4 

70 

40*4 

32*6 

40*2 

39 

1.2 

74 

31.3 

42*2 

30*3 

30 

0*3 

71.5 

23*1 

48.5 

SAMPLE  COMPUTATIONS  FOR  THE  ALBERTA  AND 
MACMICHAEL  VISCOMETERS 

(For  nomenclature  see  page  C-/g) 

Alberta  Viscometer. 

The  torsion  wire  constants  determined  by  calibration 
with  standardized  Newtonian  oils,  are  used  to  compute  the 
viscosity,  or  apparent  viscosity,  from  the  laboratory  data 
as  follows, 

K  in,.  N  rev.  cps.„  0*672  x  10  ^  lb. 

'  Q  A  '  — — — »  " 

*  A  (  sec . ) _  _ cps.  ft, sec,  (l 


(Ae  in.) 


N  rev., 
sec , 


simplifying 


_  K  (°A)  0.672  X  10-g  lb. _ 

■s»  (be)(N)  ft.  sec. 

Consider  run  no.  1  in  test  6  on  the  Saskatchewan  blend. 

The  following  is  data  obtained, 

Depth  of  immersion  4” 

Gauge  and  Wire  No.  24«»1 

Zero  Reading  298*4  °A 

Final  Deflection  Reading  316.3  °A 

Revolutions  timed  30  revs,  in  28.4  secs. 

Bob  Temperature  28  °F 

Circulating  Temperature  27  °F 

Ambient  Temperature  73  °F 

Computations 

Net  Deflection 

Because  the  direction  of  rotation  is  opposite  to  the 
degree  scale,  it  is  necessary  to  subtract 
between  the  two  readings  from  360°, 


the  difference 


Net  Defln 


s 


360  ~  (316.3  «  298.5) 
342.2  °A 


Rotary  Speed 


N  s  30  /  28.4 

N  m  1.057  rev, /sec 


Wire  Constant 


From  Table  2-C. 

K  s  7.17  (N  rev./sec. ) (he  In.) (ops.) 

°A 

Effective  Depth  of  Immersion 


he  =  4.0M  +  0.16” 


4.16” 


Apparent  Viscosity 
From  Equation  (2*c). 

7.17 (in. ) (cps. )  x  (342.2°A)  x.  0.6/2  x  10~3  lb. 

(°A)(sec.)  (cpsTT  ft. sec. 

(4.16  in.)  (1.057  rev. /sec.) 

a  0.0376  lb. /ft. sec. 

Rate  of  Shear 

From  Appendix  A,  Equation  (  IX)  the  rate  of  shear  at 
the  midpoint  of  the  annulus  is 


For  the  Alberta  Viscometer 


R1  «  0.7440  in 


r2  =  0.8815  in.  Rx  +■  R2 

2 


=  0.8128 


r 


\ 


I 

ii 
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Substituting  in  these  values,  the  rate  of  shear  at  the 
midpoint  is 

(dv/dr)m  s  (4)(2'ITN  rev . /sec ,) (0 ,7440 )2  (0.8815)2 

(0.8128)2  (0.88152  -  0.74402) 

“  36,, 6  N  /  sec »  (3«c) 

For  the  above  run,  the  rate  of  shear  at  the  midpoint  is 
(dv/dr)m  *  (36,6) (1.057)  /  sec, 

=  38.8  sec,'"1 

Oil  Temperature 

The  average  oil  temperature  is  determined  by  correcting 
the  bob  temperature  for  the  temperature  difference  between 
the  circulating  fluid  and  the  ambient  temperature. 

For  the  above  run,  the  temperature  difference  -^T,  is 

^T  ss  73  «  27 

m  46  °F 

Referring  to  Fig,  (39) 

Corr’n,  3  4-  2.8  °F 
And  the  average  oil  temperature  is 

Tave  «  28  +■  2.8 

s  30.8  °P 

MacMichael  Viscometer 

The  torsion  wire  constants,  as  determined  by  the  call** 
brat  ion,  are  listed  in  Table  1-0. 


< 


Apparent  Viscosity 

(K  (cm. ) (N  rev./sec. ) (cpa. )  ) ( °M ) ( 0 .672x10"*^ )lb. _ 

”  (  °M  ^  cpa.  ft„aec 

(he  cm.)  (  N  rev. /see.) 


In  Test  No.  5,  in  Table  (1«F),  Run  no.  1  of  the 
Saskatchewan  Blend,  the  observed  data  was, 

Depth  of  Immersion  3.6  cm. 

Gauge  and  Wire  No,  32«*4 

Zero  Reading  263.0  °M 

Deflection  261.0  °M 

Revolutions  timed,  20  revs,  in  29.2  secs. 

Bob  Temperature  47.2  °F 

Circulating  fluid  temperature  40.2°F 
Net  Deflection 

Net  Defln,  s  300  «  (263  -  261) 

9  298  °M 


Rotary  Speed 


N 


-3  20  rev.  /  29*2  sec* 

a  0,685  rev.  /  sec. 

Wire  Constant 

From  Table  1«*C,  the  average  weighted  value  is 
K 


9  3.04  (cm.)  N (rev./ sec . ) (cps.i 

°M 


Effective  depth  of  immersion 

For  all  the  experimental  tests,  the  depth  of  immersion 


was  kept  constant  at  3.6  cm 


4 


) 


l 


C-1T 


Then 


hQ  a  3.6  +  0.24 
s  3.84  cm. 


Apparent  Viscosity 

Substituting  into  Equation  (4«c) 

(3. .04  cm.  (rev. /sec.  Heps.) )  (298°M)  (0.672x10*^  )_lb _ _ 

(  om  )  ft. sec. 

(3.84  cm.)  (0.685  rev. /sec.) 

^3  r  0.232  lh./ft.sec. 


Rate  of  Shear 

For  the  MacMichael,  the  bob  and  cup  radii  are 


(dv/dr)m 


R1  a  0.394. in. 

P-2  2  0.536  in.  (Ri+-  R2)/2  *s  0.465  in. 

2(2  tt  N  rev. )  (0.536)2(0.394)2 

sec . 

(0.456)2  (0,536s  •  0..3942) 


=  19*6  N  sec."11 

For  the  above  test 

(dv/dr )m  -  13.4  see.*1 


(5-o) 


Average  Oil  Temperature 

As  a  result  of  the  temperature  gradient  across  the  annulus, 
it  was  necessary  to  determine  the  average  temperature  by  tne 
arithmetic  average  of  the  bob  and  circulating  fluid  tempera** 
ture .  For  the  above  run 

lave  -  (47.2  ^40.2)/2 


43.7  °F 


( 


C—  18 


°A 

h 

ka 

% 

°M 

N 

R1A 

rim 

T 

Tave 

(dv/dr) 

a 


NOMENCLATURE.  APPENDIX  C 

degrees  deflection  on  the  Alberta  viscometer  2<if-3600A 

depth  of  immersion  in.  or  cm. 

effective  depth  of  immersion  in.  or  cm. 

torsion  wire  constant  for  the  Alberta 

(in.)(N  rev. /sec . ) (cps . )/°A 

torsion  wire  constant  for  the  MacMichael 

(cm.)(N  rev/sec. ) (cps )/°M 

degrees  deflection  on  the  MacMichael  viscometer 
2 TT*  300°M 

J 

rotary  speed  rev./sec. 

bob  radius  on  the  Alberta  in. 

bob  radius  on  the  MacMichael  in. 

temperature  °F;  subscripts  1  and  2  indicate  bob 

and  cup  surfaces  respectively 

average  oil  temperature  °F 

rate  of  shear  sec."^ 

apparent  viscosity  lb  ../ft,  sec. 


APPENDIX  D 


THE  SOURCES  AND  GRAVITIES  OF  THE 


OILS  TESTED 
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APPENDIX  D 
TABLE  1  -  D 

Source  and  Gravities  of  Crude  Oils  Tested 


Type  of 

Crude 

Date  Rec'd 

Gravity 

°A.P.I. 

Test  Remarks 

Notation 

Oils  Tested  on  the  Pilot  Pipe  Line 

Saskatchewan  Blend 
consisting  of 

22.2-22.1 

B 

Fractions  are 
volume  percent 

39%  Fosterton 

Dec. 15,  1953 

24.5-24.0 

- 

Rec'd  from  Socony 
Vacuum  Fosterton  #1 

26.5%  Cantaur 

Dec. 15,  1953 
Feb.  9,  1954 

22.0-21.5 

- 

Rec'd  from  Socony 
Vacuum  Cantaur  #2-8B 

34*5$  Success 

Dec. 15,  1953 

22.2-21.7 

- 

Rec'd  from  Socony 
Vacuum  Success  #9-llB 

Combined  Blend 
consisting  of 

— 

24.0 

C 

Fractions  are 
volume  percent 

80%  Saskatchewan 

Blend  weathered 

sample  from 

above  makeup 

20%  Redwater 

Feb. 19,  1954 

34.6-34.1 

- 

Rec'd  from  I.P.I. 

Tank  No. 22  @  Redwater 

Fosterton 

May  26,  1954 

23.7 

D 

Rec'd  from  Socony 

Vacuum  Fosterton  No.l 

Pembina 

May  and  June 

35.9-36.0 

36.6 

E 

Rec'd  from  Imperial 
Tests  @  Excelsior  and 
also  from  Pembina  Field 

Pembina 

Tested  June  21 
Tested  June  28 

36.0 

35.7 

E 

Crude  weathered  as 

Tests  progressed 

Rotational  Viscometer 


At  the  completion  of  the  tests  on  the  pilot  pipe  line,  samples 
were  obtained  at  40°F.  to  60°F.,  and  sealed  in  five  gallon  cans. 
However  they  werg  stored  in  a  room  in  which  the  temperature  varied 
from  70°F.  to  85  F.,  and  they  were  opened  often  to  remove  samples  for 
testing.  Another  sample  of  Redwater  crude  was  obtained  from  the 
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Redwater  feed  line  to  the  Interprovincial  Pipe  Line  Co.  tanks  in 
Edmonton,  for  use  on  the  rotational  viscometers.  The  gravity  as 
determined  by  a  hydrometer  was  35.0°  A.P.I. 

Determination  of  the  Specific  Gravity  of  the  Oils  Tested 

The  relationship  of  the  specific  gravity  and  the  temperature 
for  the  Fosterton,  Success,  Cantaur  crudes  and  the  Saskatchewan  Blend 
and  the  Combined  Blend,  were  determined  on  the  Westphal  chain  gravi¬ 
metric  balance.  The  best  fit  curves  of  the  data  listed  in  Table 
2-D  are  presented  in  Fig.  (50).  The  procedure  consisted  of  obtaining 
a  sample  after  the  blends  and  crudes  had  been  circulated  thoroughly 
and  several  test  runs  had  been  made.  The  sample  was  split  and  one 
sample  was  cooled  immediately,  while  the  other  one  was  tested  and 
cooled  in  successive  steps.  The  cooled  sample  was  then  tested  at 
successive  increments  in  the  temperature.  The  instrument  was  checked 
with  distilled  water  before  and  after  each  test.  During  the  testing 
program  the  instrument  was  also  checked  with  Toluene  of  S.G.  0.8660 
at  20°C.  The  gravimetric  balance  can  also  be  checked  very  closely  by 
suspending  the  bob  in  air. 

The  data  obtained  on  the  Fosterton  crude  indicates  that  the 
oil  was  weathering  at  the  higher  temperatures.  The  Pembina  crude 
weathered  so  quickly  at  the  higher  temperatures  that  no  consistent  data 
was  obtained.  At  the  lower  temperatures  the  Pembina  crude  was  so 
thixotropic  that,  the  instrument  could  not  move  freely.  The  gravity 
was  determined  with  a  hydrometer  with  several  observations  at  various 
temperatures.  The  standard  API  gravity  temperature  tables  were  usea 
to  plot  the  specific  gravity  as  a  function  of  temperature. 
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FIG.  50 

v  THE 

VARIATION 

OF  THE  DENSITY  FOR 

THE  OILS 

TESTED 

o 

SA  SK.  BLEN  D  22  °  API 

• 

COMBINED 

BLEND  2  4  °  API 

o 

FOSTERTON 

CRUDE  23.  7°  API 

• 

PEMBINA 

CRUDE  3  5.  7°  A  P  1 

o 

R  EDWATER 

CRUDE  3  5.1  °  A  PI 
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APPENDIX  D 

TABLE  2  -  D 


Observed  Gravities  of  the  Oils  Tested 
Specific  Gravity  of  the  Saskatchewan  Blend 


Temp. 

Specific 

Temp. 

Specific 

°F 

Gravity 

V 

Gravity 

54.5 

0.9225 

Feb. 19,  1954 

45 

0.9215 

57.9 

0.9214 

50.5 

0.9200 

58.8 

0.9212 

53.0 

0.9191 

63.3 

0.9199 

55 

0.9188 

64.3 

0.9194 

63.5 

0.9153 

66.2 

0.9187 

65.3 

0.9152 

67.5 

0.9182 

68.0 

0.9145 

71.0 

0.9170 

71.0 

0.9135 

71.9 

0.9165 

72.3 

0.9125 

77.0 

0.9151 

75.8 

0.9118 

80 

0.9143 

79.5 

0.9105 

81.5 

0.9137 

84.0 

0.9131 

Gravity 

60/60  22 

.2  °API 

22, 

.1  °AP1 

Specific  Gravity 

of  the  Foster  ton  Crude 

Temp. 

Specific 

Temp. 

Specific 

Gravity 

Gravity 

76 

0.9032 

May  26 

40 

0.9180 

60 

0.9093 

47 

0.9146 

40 

0.9185 

53.5 

0.9130 

42.5 

0.9167 

58 

0.9113 

54 

0.9120 

60 

0.9100 

57 

0.9112 

70 

0.9069 

60 

0.9105 

74 

0.9053 

64 

0.9090 

83 

0.9028 

68 

0.9078 

O 

72 

0.9070 

Gravity  6O/6O 

24.0 

API 

77 

0.9059 

March  1 


After  Tests 
Temp.  Specific 
Gravity 

63.2  0.9070  June  12 

65.7  0.9060 

68.6  0.9052 

70.4  0.9049 

72.7  0.9041 

75.9  0.9032 

79.9  0.9022 


24.2  °API 


J 


■ 


. 
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TABLE  2  -  D 


Specific  Gravity  of  the  Combined  Blend 


Temp. 

Specific 

Temp. 

Specific 

F 

Gravity 

Op 

Gravity 

70.9 

0.9042 

34.8 

0.9170 

59.5 

0.9080 

37.0 

0.9158 

56.0 

0.9088 

41.0 

0.9149 

52.3 

0.9110 

44.0 

0.9134 

52.0 

0.9105 

46.9 

0.9128 

49.9 

0.9119 

53.0 

0.9115 

44.9 

0.9138 

57.2 

0.9100 

44.0 

0.9145 

60.5 

0.9085 

40.2 

0.9160 

66.3 

0.9073 

69.5 

0.9066 

73.0 

0.9049 

78.0 

0.9033 

O 

80.6 

0.9020 

Gravity 

60/60  24.0  API 

87.8 

0.8997 

92.8 

0.8950 

Gravity 

of  the  Pembina  Crude 

Temp. 

Observed  Corrected 

Temp. 

Observed 

Corrected 

oF 

Gravity  Gravity 

Op 

Gravity 

Gravity 

°API  °API 

°API 

°API 

June  18 3 

,  1954 

June  28 

,  1954 

54 

35.15  35.55 

59.0 

35.6 

35.7 

54 

35.2  35.6 

59.2 

35.65 

35.73 

53.5 

35.25  35.70 

59.5 

35.7 

35.75 

60.0 

35.7 

35.7 

Average 

weighted  gravity  35.6  API 

Average  value  35.7°API 

Gravity  Determination  of  the  Redvater  Crude 

Gravity  by  Hydrometer  Feb.  19 

34.6  °API 

Temp. 

Specific 

Temp . 

Specific 

op 

Gravity 

Op 

Gravity 

41 

0.8618 

74.2 

0.8468 

March  5 

44.5 
49.8 
56.0 
57.0 

61.5 
62.0 
69.0 
74.0 
79.0 


0.8600 

0.8585 

0.8562 

0.8560 

0.8530 

0.8529 

0.8518 

0.8507 

0.8502 


66.7 

63.0 

57.8 

53.5 
44.  8 

41.5 


0.8504 

0.8530 

9.8560 

0.8586 

0.8623 

0.8641 


Average  Gravity  60/60  34.1  API 


1954 


' 


V 


■ 


This  sample  was  used  to  make  up  the  Combined  Blend 
Sample  obtained  for  rotational  viscometer  tests  from  Inter' 
Provincial  Pipe  Line  Go.  feed  line  from  Red water  to  duplicate 
sample  No.  2  tested  by  A.  Masuda. 

Gravity  by  hydrometer  60/60  is  35.1°API 
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THE  OBSERVED  DATA  AND  COMPUTATIONS  ON  THE  PILOT 


PIPE  LINE  FOE 
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APPENDIX  1  -  E 


Data  Obtained  on  The  Pilot  Pipe  Line 


The  procedure  used  in  obtaining  the  data  was  changed  with 
different  crudes  and  as  a  result  of  the  previous  experience  obtained* 
For  the  Saskatchewan  Blend,  the  corrected  downstream,  temperature  has 
been  listed*  For  the  other  crudes,  this  correction  was  incorporated 
in  the  formula  for  the  average  temperature  of  the  oil  (Appendix  B)» 

The  rate  of  flow  is  in  cubic  feet  per  minute  (cfm),  the  volume  cor¬ 
rected  for  the  temperature  difference  from  the  measuring  tank  to  the 
test  line  temperature*  For  the  Saskatchewan  Blend  the  temperature  of 
the  oil  in  the  measuring  tank  was  determined  at  several  points  on  the 
tank*  This  temperature  was  later  obtained  on  the  pipe  carrying  the 
test  oil  into  the  measuring  tank,  where  the  transfer  of  heat  was  better, 
and  a  more  representative  temperature  could  be  obtained* 

All  tests  on  the  Saskatchewan  Blend,  Combined  Blend  and  the 
Fosterton  Crude  were  in  the  laminar  region  of  flow.  For  the  Pembina 
Crude,  it  was  necessary  to  compute  the  apparent  viscosity  for  both 
the  laminar  and  turbulent  region.  For  the  runs  in  which  the  flow  could 
be  either,  the  laminar  calculations  are  first  and  the  turbulent  calcula¬ 
tions  are  on  the  line  below* 
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Observed  Data  and  Computations 
for  the  Saskatchewan  Blend  on  the  Pilot  Pipe  Line 
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Observed  Data  and  Computations 
for  the  Pembina  Crude  on  the  Pilot  Pipe  Line 
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TABLE  6  -  e 


Summary  of  Apparent  Viscosity  Data 
for  the  Saskatchewan  Blend 

Composite.  Data  of  the  Pilot  Pipe  Line  and  Alberta  Viscometer: 


Reciprocal 

of 


Apparent 

Kinematic 

Saybolt 

!emp. 

Density 

Viscosity 

Viscosity 

Universal 

°F 

Ib./i'U3 

lb. /ft. sec. 

sec. /ft. ^ 

sec. 

10 

50.5 

1.10 

53.2 

8,000 

20 

58.3 

0.60 

97.1 

4,450 

30 

58.0 

0.35 

166.0 

2,600 

AO 

57 » 8 

0.205 

282 

1,540 

50 

57.6 

0.125 

460 

940 

60 

57*4 

0.08 

717 

600 

70 

57.2 

0.055 

1040 

420 

80 

57.0 

0.0375 

1520 

287 

TABLE 

9  -  E 

Apparent  Viscosity  Data 

for  the  Combined 

Blend 

Combined  Data  of 

the  Alberta 

and  Pilot  Pipe  Line 

Viscometers: 

10 

57.7 

0.83 

69.5 

6,100 

20 

57.5 

0,36 

160 

2,700 

30 

57.3 

0.18 

318 

1,360 

40 

57.1 

0.108 

529 

820 

50 

56.9 

0.072 

790 

545 

60 

56.7 

0.050 

1,133 

380 

70 

56.5 

0.034 

1,660 

260 

80 

56.3 

0.0245 

02,300 

190 

TABLE 

10  -  E 

Apoarent  Viscosity  Data  for  the  Fosterton 

Crude 

Combined  Dat?  of 

the  Alberta 

and  Pilot  Pipe  Line 

Viscometers: 

10 

58.0 

0.55 

105.3 

4,050 

20 

57.8 

0.28 

206 

2,100 

30 

57.6 

0.163 

354 

1,230 

40 

57.3 

0.10 

573 

760 

50 

57.1 

0.065 

879 

495 

60 

56.8 

0.045 

1,263 

345 

70 

56.6 

0.031 

1,828 

240 

80 

56.3 

0.0227 

2,480 

177 

^  r 
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TABLE  ll  -  E 


Data 

Obtained 

Apparent  /iscositv  Data  for  the  Pembina  flmHi 

e 

at  a  Shear  Rate  of  3.000  see.-l 

on  the  Pilot  Pine 

T.i  rie  ; 

Reciprocal 

of 

Apparent 

Kinematic 

Saybolt 

Temp. 

Densitv 

Viscosity 

Viscosity 

Universal 

op 

lb/ft0^ 

lb. /ft. sec. 

sec. /ft. ^ 

sec. 

20 

53.6 

47.0 

1,140 

380 

30 

53*4 

32.0 

1,670 

260 

40 

53.2 

19.0 

2,800 

155 

50 

53.0 

11.6 

4,560 

100 

60 

52.7 

7.0 

7,530 

68 

70 

52.5 

5.1 

10,300 

55 

80 

•  52.2 

4.2 

12,420 

50 

TABLE  12 

-  E 

Data 

Obtained 

at  a  Shear 

Rate  of  80.000  sec.“^  on  the  Pilot  Pino  Trine: 

40 

53.2 

12.5 

4,250 

106 

50 

53.0 

8.6 

6,160 

79 

60 

52,7 

6.0 

8,770 

61 

70 

52.5 

4.5 

11,670 

51 

80 

52.2 

3.9 

13,400 

48 

TABLE  13 

-  E 

Data 

Obtained 

at  a  Shear 

Rate  of  64  sec.”^  on 

the  Alberta  Viscometer: 

20 

53.6 

0,78 

68.6 

6,200 

30 

53.4 

0.46 

116 

3,700 

40 

53.2 

0.14 

380 

1,300 

50 

53.0 

0.037 

1,430 

305 

60 

52.7 

0,0134 

3,930 

112 

70 

52.5 

0.0076 

6,900 

72 

80 

52.2 

0.0057 

9,150 

59 

17  -  E 


TABLE  14  ~  E 


Apparent  Viscosity  Data  for  the  Redwatcr  Crude 


This 

apparent  Viscosity  Data 

was  Obtained 

on  the  Pilot  Pipe 

Line  by  A.  Masuda 

Sample  No .  1 

Apparent 

Sample  No.  2 

Apparent 

Average 

Apparent 

Saybolt 

Temp  ■ 

.  Viscosity 

Viscosity 

Viscosity 

Universal 

op 

lb. /ft. sec. 

lb. /ft. sec. 

lb./ft.sec. 

sec. 

20 

19.72 

25.1 

22.4 

159 

25 

14.57 

17.60 

16.09 

119 

30 

12.49 

32.6 

12.55 

104 

35 

10.87 

10.41 

10.64 

93 

40 

9.50 

8.86 

9.18 

84 

45 

8.25 

7.67 

7.96 

75 

50 

•  7.09 

6.70 

6.90 

67 

55 

6.01 

5.93 

6.97 

60 

60 

5.05 

5.20 

5.13 

55 

65 

4.49 

4.56 

4.53 

52 

70 

4.21 

4.15 

4.18 

50 

75 

3.87 

3.77 

3.82 

48 

80 

3.56 

3.38 

3.47 

46 

APPENDIX  F 


THE  OBSERVED  DATA  AND  COMPUTED  RESULTS  ON  THE 
MACMICHAEL  VISCOMETER  FOR 


TABLE 

1-F. . , 

TABLE 

2-F. . , 

...COMBINED  BLEND 

TABLE 

3-F... 

...FOSTERTON  CRUDE 

TABLE 

4-F... 

. . .PEMBINA  CRUDE 

TABLE 

5-F. . , 

. . .RED WATER  CRUDE 

. 


....  • 


•• 


wi  I 


IF 


TABLE  1  -  F 


Test 


Wire 

Gauge 


1  34-1 


34-1 


34-1 


4  34-1 


32-4 


30-1 


S  A  S  K 

A  T  C 

H  E  W  A  N 

B 

LEND 

All  runs 

are  at 

,  an  immersion  of  3*6  cm, 

Avg. 

Net 

Rotary 

Shear 

App. 

Bob 

Circ. 

Oil 

Defln* 

Sneed 

Rate 

Vise, 

Temp. 

Temp. 

Temp. 

% 

im 

1/sec. 

Ib./ft.sec. 

Op 

°F 

op 

165 

0.741 

14*5 

0.0448 

80*6 

77.0 

78.8 

29 

0.118 

2.3 

0.0495 

73*5 

0.34 

6.6 

0.0465 

118*5 

0.54 

10.6 

0.0442 

142*5 

0*667 

13*1 

0.0430 

160*5 

0.735 

14.4 

0o  044 

30 

0.125 

2*4 

0.0483 

78.8 

181 

O.667 

13.1 

0*0546 

70,8 

69.4 

70.1 

199 

0*665 

13*0 

0*0603 

71*2 

69.0 

70.1 

159 

0.513 

10*0 

0. 0624 

103.5 

0*317 

6.2 

0.0658 

41 

0*116 

2*3 

0.0712 

104 

0*327 

6.4 

0.064 

165*5 

0.531 

10.4 

0.0628 

69*0 

230.5 

0,764 

14.9 

0*0609 

71*2 

70.1 

207 

0*525 

10.3 

0.0794 

63.6 

57.8 

60.6 

275.5 

0,672 

13*1 

0.0825 

281 

0,669 

13*1 

0.0845 

37+0 

0.319 

6*2 

0.0844 

60,6 

281.5 

0.672 

13*1 

0*0844 

372 

0*602 

11.8 

0.124 

54 

46 

50 

331 

0.521 

10.2 

0.128 

285 

0.433 

8,5 

0*133 

145 

0.221 

4*3 

0.132 

80 

0*126 

2*5 

0.128 

50 

418 

0*671 

13*1 

0.125 

298 

0,685 

13*4 

0*232 

47*2 

40.2 

43*7 

191 

0.431 

8*4 

0.236 

99 

0.222 

4.3 

0.238 

56 

0.123 

2*4 

0*244 

100 

0*222 

4*3 

0*24 

309 

0.701 

13*7 

0*235 

222 

0.79 

15*5 

0.39 

40 

27 

33*5 

3.50 

0*521 

10*2 

0.399 

66 

0.216 

4*2 

0*425 

42.5 

0.132 

2.6 

0.446 

160 

0*534 

10.4 

0.415 

234.5 

0.779 

15*2 

0*418 

33,5 

44.5 

0*131 

25*6 

0.470 

15  mins, 


later 


a 


2F 


Wire 

Avg. 

Test 

Gauge 

Net 

Rotary 

Shear 

App. 

Bob 

Circ. 

Oil 

No, 

&  No. 

Defln. 

Speed 

Rate 

Vise. 

Temp. 

Temp. 

Temp. 

Remarks 

on 

1/sec . 

3b. /ft.  sec. 

op 

£F 

7 

30-1 

299 

0.585 

11.4 

0.709 

39 

17.6 

28.4 

237 

0.433 

8.5 

0.759 

322 

0.516 

10.1 

0.866 

36 

17.6 

26.8 

115 

0.151 

3.0 

1.06 

177 

0.247 

4.8 

0.994 

297 

0.447 

8.8 

0.92 

Oil  in  Annulus 

343 

0.532 

10.4 

0.895 

26.8 

22°F 

8 

28 

298 

0.785 

15.4 

1.02 

31 

8.6 

19.8 

215 

0.518 

10.1 

1.12 

156 

0.339 

6.6 

1.24 

110 

0.213 

4.2 

1.39 

74 

0.127 

2,5 

1.56 

286 

0.664 

13.0 

1.16 

Oil  in  Annulus 

331 

0.796 

15.6 

1.12 

19.8 

18.5°F 

TABLE  2  - 

F 

C 

OMBINED 

B  L  E 

N  D 

9 

34-1 

16,5 

0.117 

2.3 

0.0284 

81 

79.4 

80.2 

63.5 

0,442 

8.7 

0.0289 

85.6 

0.602 

11.8 

0.0286 

114 

0.813 

15.9 

0.0282 

80,2 

10 

34-1 

113.5 

0.669 

13.1 

0.0342 

71.5 

70 

70,7 

116.5 

0.669 

13.1 

0.0351 

60 

0.32 

6.3 

0.0378 

23 

0.121 

2,4 

0.0382 

95.5 

0.535 

10.5 

0.036 

136.5 

0.765 

15.0 

0.0359 

22 

0.116 

2.3 

0.0383 

70.7 

11 

34-1 

155.9 

0.667 

12,9 

0.0471 

64 

60 

62 

177.9 

0.773 

15.1 

0.0464 

130.9 

0.514 

10.1 

0.0512 

87.9 

0.326 

6.3 

0.0543 

62 

47.9 

0.171 

3.3 

0.0566 

12 

34-1 

262.5 

O.669 

13.1 

0.079 

54 

46.2 

50.1 

311.5 

0.802 

15.7 

0.078 

212.7 

0.518 

10.1 

0.0826 

181,5 

0.425 

8.3 

0.086 

144.5 

0.324 

6.3 

0.0896 

50.7 

0.118 

2.3 

0.0862 

103.5 

0.225 

4.4 

0,0925 

50,1 

3F 


Wire 

Avg. 

Test 

Gauge 

Net 

Rotary 

Shear 

App. 

Bob 

Circ. 

Oil 

No, 

&  No. 

Defln. 

Speed 

Rate 

Vise. 

Temp. 

Temp. 

Temp, 

rps 

1/sec. 

lb, /ft. sec. 

op 

2F 

op 

13 

32-4 

162.2 

0.76 

14.9 

0.114 

49 

40 

44.5 

115 

0.515 

10.1 

0.119 

49 

0.215 

4.2 

0.122 

26.5 

0.118 

2.3 

0.120 

44.5 

14 

30-1 

72.5 

0.1195 

2.3 

1.16 

34 

14.9 

24.5 

340 

0.801 

15.7 

0.588 

290 

0.665 

13.0 

0.604 

24.5 

74 

0.109 

2.1 

0.769 

32.5 

34.9 

23.7 

132 

0.2195 

4.3 

0.834 

91 

0.119 

2,3 

1.06 

207 

0.344 

6.7 

0.835 

255 

0,45 

8.8 

0,786 

332 

0.665 

13.0 

0.693 

380 

0.789 

15.4 

0.667 

23.7 

15 

30-1 

132.5 

0.699 

23.7 

0,262 

40 

26 

33 

68 

0.328 

6.4 

0.288 

32 

0.149 

3.9 

0.298 

91 

0.451 

8.8 

0.28 

131 

O.676 

13.2 

0.254 

33 

Remarks 


TABLS  3  -  F 


16  34-1 


17  34*1 


F  0 

S  T  E 

RT0N 

C  R  U 

D  E 

93.5 

0.702 

13.7 

0.0268 

79 

81 

80 

101 

0.722 

14.1 

0.0282 

95.5 

0.652 

12.8 

0.0294 

48.8 

0.32 

6.3 

0.0307 

17.7 

0.117 

2.3 

0,0304 

49 

0.341 

6.7 

0.0289 

80 

94.2 

0.67 

13.1 

0.0283 

120.5 

0.771 

15.1 

0.0315 

70 

69.3 

69.7 

106.5 

0,655 

12.8 

0.0328 

53 

0.314 

6.1 

0.0340 

29.3 

0.169 

3.3 

0.0350 

63 

0.376 

7.4 

0.0338 

69.7 

113.8 

0.667 

13.0 

0.0342 

106 

0.449 

8,8 

0.0474 

57.4 

63.0 

60.2 

127 

0.517 

10.1 

0.0494 

153 

0,63 

12.3 

0.0489 

173.5 

0.718 

14.0 

0.0486 

IS  34-1 


4F 


Wire  Avg. 


Test 

Gauge 

Net 

Rotary 

Shear 

App. 

Bob 

Circ. 

Oil 

No. 

&  No. 

Defln. 

Speed 

Rate 

Vise. 

Temp. 

Temp. 

Temp. 

om 

£Es 

1/sec. 

lb. /ft. sec. 

°F 

op 

°F 

19 

34-1 

53.5 

0,593 

11.6 

0.0482 

63.5 

56 

59.8 

61.5 

0.679 

13.3 

0.0484 

70.5 

0.723 

14.1 

0.052 

20 

34-1 

163 

0.661 

12.9 

0,0496 

64 

58 

61 

188.5 

0.78 

15,3 

0.0486 

155 

O.665 

13.0 

0.0469 

125 

0.515 

10.1 

0.0488 

21 

34-1 

176 

0.662 

12.9 

0,0535 

62 

56 

59 

83.5 

0.324 

6.3 

0.0519 

32.5 

0.1227 

2.4 

0.0534 

146.5 

0.532 

10.4 

0.0554 

181,5 

0.704 

13,8 

O.O519 

22 

34-1 

203 

0,61 

11.9 

0.067 

53.6 

47.7 

50.6 

178 

0.534 

10.4 

0.0671 

70 

0.222 

4.3 

0.0635 

39.5 

0,1208 

2.4 

0.0659 

244 

0,441 

8.6 

0.0658 

50.6 

23 

32-4 

140.5 

0.78 

15  #3 

0.096 

41.2 

36 

41.6 

123.5 

0.658 

12. 9 

0.10 

99 

0.516 

10.1 

0.102 

83.5 

0.431 

8.4 

0.104 

63.5 

0.322 

6,3 

0.105 

43.3 

0,216 

4,2 

0.107 

157 

0.781 

15.3 

0.107 

41.6 

24 

32-4 

273.5 

0.786 

15.4 

0.186 

37 

28.4 

32.7 

234.5 

0,658 

12,9 

0.190 

184 

0.513 

10.0 

0.192 

112 

0,324 

6.3 

0.185 

76.5 

0,22 

4,3 

0.186 

45.5 

0.129 

2.5 

0.3-88 

25 

30-1 

146 

0,781 

15.3 

0,259 

37 

19 

28 

100 

0.519 

10.1 

0,267 

82 

0.415 

8.1 

0.274 

64 

0.328 

6,4 

0.27 

43 

0,217 

4.2 

0.275 

27 

0.132 

2.6 

0.285 

68.5 

0,342 

6,7 

0.278 

104.5 

0.545 

10.7 

0,266 

127.5 

0,662 

12.9 

0.267 

28 

149.5 

0.778 

15.2 

0,267 

37 

19 

Remarks 


5F 


TABLE  U  -  F 


Test 

No. 


26 


27 


28 


29 


30 


31 


P 

EMBI 

N  A  C 

RUD 

E 

Wire 

Avg. 

Gauge 

Net 

Rotary 

Shear 

App. 

Bob 

Circ. 

Oil 

&  No. 

Defln, 

Soeed 

Rate 

Vise. 

Temp. 

Temp. 

Temp. 

1/sec. 

Ib./ft^ec. 

op 

of 

of 

34-1 

76 

0.526 

10,3 

0.0291 

61 

59 

60 

40 

0.1276 

2,5 

0.0632 

57 

0.228 

4,5 

0.0504 

88.5 

0.535 

10.5 

0.0333 

102 

0.76 

14,9 

0.027 

35 

0.128 

2.5 

0.055 

60 

34-1 

25 

0. 663 

13,0 

0.0076 

79 

80.8 

79,9 

4 

0.102 

2.0 

0,0079 

12,5 

0.335 

6.6 

0.0075 

21 

0.529 

10.4 

0.0080 

31 

0.765 

15.0 

0.00816 

79,9 

34-1 

15 

0.322 

6.3 

0,00939 

70.4 

69.8 

70.1 

24.5 

0.526 

10.3 

0.00939 

26 

0,529 

10.4 

0.0099 

.  37 

0.679 

13,2 

0.011 

20 

0.321 

6.3 

0.0126 

8.5 

0.123 

2,4 

0.0139 

70.1 

34-1 

251.3 

0.755 

14,8 

0.067 

54,5 

46 

50.2 

254.3 

0.758 

14,9 

0,0676 

242.3 

0.666 

13.0 

0.0734 

171.3 

0.322 

6.3 

0.107 

34-2.3 

0.224 

4,4 

0.128 

124.3 

0.129 

2.5 

0.194 

137.3 

0,131 

2.6 

0.211 

304.3 

0.793 

15,5 

0.0773 

50,2 

223.7 

0.798 

15.6 

0.0564 

32-4 

186 

0.795 

15.5 

0.125 

49 

39-2 

44,1 

168 

0.654 

12,8 

0,137 

149 

0.51 

10.0 

0.156 

122 

0.318 

6.2 

0.205 

101 

0.218 

4.3 

0.248 

44,1 

81 

0,12 

2.3 

0.360 

28-1 

116 

0.775 

15.2 

0,404 

46 

27 

36,5 

78.5 

0.51 

10.0 

0,416 

49 

0.314 

6.1 

0.421 

34 

0.234- 

4.2 

0,429 

36.5 

18 

0.12 

2.3 

0.556 

Remarks 


Light  ends 
being  lost 


Stopped  for  2 
mins. 
3  minso  later 


c 


1 


6f 


Wire 


Test 

Gaupe 

Net 

Rotary  Shear 

App. 

Bob 

Circ. 

Avff. 

No. 

&  NOo 

Defln. 

Speed 

Rate 

Vise. 

Temp.  Oil  Temp.  Remarks 

om 

1/sec. 

Ib«/ft.sec 

t  °F 

2f 

Op 

32 

30-1 

168.0 

0.776 

15.0 

0.585 

42 

19.4 

30.7 

116.3 

0.518 

10.2 

0.607 

76.5 

0.317 

6.2 

0.652 

54.0 

0.218 

4.3 

0.669 

31.0 

0.122 

2.4 

0.686 

\ 

124.0 

0.534 

10.5 

0.627 

171*0 

0.785 

15.4 

0.588 

30.7 

33 

28-1 

291.0 

0.79 

15*5 

0.995 

9.5 

29.5 

19.5 

218.0 

0.518 

10.2 

1.14 

110.0 

0.217 

4.2 

1.37 

74*0 

0.129 

2.5 

1.55 

311.0 

0.772 

15.1 

1.09 

19.5 

TABLE  5  • 

-  F 

R  I 

SDWATER 

C  R  U  D 

E 

34 

34-2 

4*5 

0.1713 

3.4 

5.17 

78 

77*5 

77.8 

7.5 

0.392 

7.7 

3.77 

14*0 

0.572 

11.2 

4*  82 

8.0 

0..325 

60  A- 

4.85 

11.5 

0.443 

8.7 

5.12 

78 

77.8 

35 

34-1 

17*7 

O.664 

13.0 

5.36 

70 

70.5 

70.2 

19.0 

O.664 

13.0 

5*76 

15.0 

0.517 

10.1 

5*84 

36 

34-2 

9.5 

0.296 

5*8 

6.32 

64*5 

60.8 

62.7 

16.5 

0.53 

10.4 

6.12 

16.5 

0.519 

10.2 

6.26 

16.0 

0.522 

10.2 

6.03 

5.0 

0.187 

3.7 

5.27 

63.2 

58.6 

60.9 

34.0 

0.445 

8.7 

6.17 

18.0 

0.55 

10.8 

6.46 

37 

34-2 

9.5 

0.208 

4.1 

8.53 

53*2 

50.0 

51.6 

15.5 

0.401 

7.8 

7.6 

22.0 

0.532 

10.4 

8.15 

Bob  off 

21.5 

0.51 

10.0 

8.3 

10.7 

0.246 

4*8 

8.57 

38 

34-1 

44.0 

0.787 

15.4 

11.2 

46.0 

39.2 

42.6 

39*0 

0.66 

12.9 

11.9 

31.5 

0.516 

10.1 

12.3 

6.5 

0.116 

2,3 

11.3 

42,6 

31.5 

0.535 

10.5 

11.9 

c 


)■ 


7F 


Wire  Avg. 


Test 

Gauge 

Net 

Rotary 

Shear  App. 

Bob 

Circ. 

Oil 

Mo. 

&  No. 

Defln. 

Speed 

Rate 

Vise. 

Temp. 

Temp. 

Temp. 

°M 

ros 

1/sec. 

Ih/ft.sec. 

°F 

£F 

Op 

39 

34-2 

18 

0.206 

4.1 

17.2 

37 

27.5 

32.2 

36.5 

0.45 

8 .8 

16.0 

56 

O.664 

13.0 

16.6 

66.5 

0.785 

15.4 

16.7 

30 

0.344 

6.7 

17.2 

12 

0.1317 

2.6 

18.0 

32.2 

40 

34-2 

86 

0.61 

11.9 

27.8 

31.5 

17.6 

24.6 

95 

0.689 

13.5 

27.2 

62 

0.425 

8.3 

28.8 

49 

0.327 

6.4 

29*5 

36 

0.222 

4.3 

31.8 

23 

0.131 

2.6 

34*6 

24.6 

107 

0.675 

13*2 

31.3 

41 

34-2 

94*5 

0.662 

12.9 

28.2 

35.6 

9.5 

22.6 

88 

0.657 

12.9 

26.4 

58 

0.41 

8.0 

27.9 

50.5 

0.306 

6.0 

32.5 

19 

0.116 

2.3 

32.3 

22.6 

98 

0.693 

13.5 

27.9 

42 

32-4 

67.7 

0.675 

13*2 

53.6 

31*5 

8.6 

20 

53.5 

0.528 

10.3 

54.1 

24.8 

0.222 

4.3 

59*6 

15 

0.1213 

2.4 

66.0 

20 

81 

0.797 

15.6 

54.2 

Remarks 


c 


APPENDIX  G 


THE  OBSERVED  DATA  AND  COMPUTED  RESULTS  ON  THE 
ALBERTA  VISCOMETER  FOR 


TABLE 

1-G. . , 

TABLE 

2-G... 

...COMBINED  BLEND 

TABLE 

3-G... 

...FOSTERTON  CRUDE 

TABLE 

4-G. . , 

. . .PEMBINA  CRUDE 

TABLE 

5-G... 

...REDWATER  CRUDE 

1G 


TABLE  1  G 


SASKATCHEWAN  BLEND 


Wire 


Test 

Gauge  3mm.  Net 

Rotary 

Shear 

No. 

&  No. 

in  Defln. 

Speed 

Rate 

2A 

rps 

1/sec. 

1 

26-2 

4"  7,4 

0.0997 

3.6 

72.0 

0.945 

34.8 

109.0 

1.42 

52.0 

2 

26-2 

4"  96.4 

0.939 

34.3 

10.6 

0.102 

4.1 

98.1 

0.946 

34.6 

3 

26-2 

4"  15.0 

0.0993 

3.6 

32.5 

0.215 

7.9 

158.8 

1.058 

38.7 

326.0 

2.16 

79.0 

4 

26-2 

3"  19.0 

0.1014 

3.7 

109.2 

0.585 

21.4 

263.5 

1.41 

51.6 

399.5 

2.14 

78.3 

5 

24-1 

4"  250.3 

1.76 

64.4 

201.5 

1.41 

51.6 

16.0 

0.103 

3.8 

91.2 

0.583 

21.4 

272.0 

1.76 

64.4 

327.0 

2.12 

77.6 

36.3 

0.222 

8.1 

6 

24-1 

4"  342.2 

1.06 

38.8 

79.0 

0.216 

7.9 

36.0 

0.0964 

3.5 

384.5 

1.17 

4  2.8 

n 

( 

22-1 

4"  294.7 

1.41 

51.6 

139.5 

0.583 

21.4 

32.5 

0.0983 

3.6 

66.3 

0.222 

8.1 

368.0 

1.65 

6O.4 

8 

22-1 

3"  335.3 

1.4 

51.2 

158.5 

0.58 

21.2 

33.0 

0.0886 

3,2 

35 

0.0886 

3.2 

22-1 

3"  167.0 

0.575 

21.0 

344.0 

1.39 

50.9 

Avg. 

App.  Bob  Circo  Oil 

Vise®  Temp®  Temp®  Temp.  Remarks 

Iba/ft®  sec®  of  °F  °F 


0.0382 

0.0392 

0.0395 

78.5 

81.5 

78.5 

0.0529 

0.0535 

0.0534 

70.8 

70.2 

71.0 

0.0778 

0.0779 

0.0774 

0.0777 

60.3 

58.6 

61.5 

0.126 

0.127 

0.127 

0.127 

49.5 

47.7 

51.3 

0.165 

0.165 

44.5 

39.2 

46.5 

0.18 

0.181 

0.179 

0.178 

41.6 

39.0 

43.6 

0.189 

41.2 

39.0 

43.2 

0.376 

0.424 

0.433 

0.380 

23.0 

27.0 

30.8 

0.682 

0.754 

22.5 

20.3 

25,9 

1.04 

0.942 

0.701 

21.5 

20.0 

24.9 

0.992 

1.13 

1.54 

1.63 

1.20 

1.03 

16.5 

12.0 

19.1 

19.1 

19.1 

2G 


TABLE  2  -  G 


COMBINED  BLEND 


Wire  Avg. 


Test 

Gauge 

,  Net 

Rotary 

Shear 

App. 

Bob 

Circ. 

Oil 

No, 

&  No- 

in 

Defln, 

Speed 

Rate 

Vise. 

Temp.,. 

Temp. 

Temp. 

2A 

rps 

1/ sec. 

lb, /ft.  sec. 

Op 

op 

OF 

9 

30-3 

3" 

335*6 

1.42 

52.0 

0.0246 

80.5 

81.0 

80.5 

254.5 

1.07 

39*2 

0,0249 

54*2 

0,226 

8.3 

0.025 

10 

26-2 

4" 

45*3 

0.705 

25.8 

0.0332 

70,5 

69.4 

69*9 

7.0 

0.1027 

3.3 

0.0351 

69*5 

69.4 

69.5 

130.1 

1.905 

69.7 

0.0351 

11 

26-2 

4" 

105.  a 

1.061 

38.8 

0.0513 

61.5 

60,0 

62.1 

30.5 

0,216 

7.9 

0.0728 

20.2 

0.0995 

3.6 

0.1047 

212.0 

2.14 

78.3 

0.051 

12 

26-2 

4" 

314.0 

2.15 

73.8 

0.0752 

49.5 

43.0 

51*3 

162.5 

1.06 

33.3 

0.0786 

91.0 

0.533 

21.3 

0.0305 

17.0 

0.101 

3.7 

0.0865 

49.5 

47*8 

51*3 

13 

26-2 

3" 

35.3 

0.220 

3.1 

0.109 

41*6 

40.1 

44*1 

163.3 

1.06 

38.8 

0.105 

325,5 

2,14 

73,3 

0,1035 

215.3 

1,42 

52.0 

0.103 

41.6 

40.0 

43.7 

14 

24-1 

4" 

23*7 

0.098 

3.6 

0.326 

28.8 

27*5 

31.9 

52.2 

0.22 

3.1 

0.276 

131*2 

0.534 

21.4 

0.26 

359*2 

1.76 

60.4 

0.237 

49*7 

0,213 

8.0 

0,264 

51*7 

0,213 

8.0 

0.275 

31*9 

24,7 

0.0959 

3*5 

0,299 

27*5 

15 

24-1 

3" 

23.5 

0.0938 

3.6 

0.44 

26.0 

22.0 

29*4 

60.3 

0.213 

8.0 

0.421 

370,0 

1.76 

64.4 

0.32 

237*3 

1,06 

33.8 

0.342 

22.0 

29.4 

29*0 

0.0986 

3,6 

0.449 

16 

22-1 

4" 

30.7 

0.0973 

3*6 

1.02 

17*5 

14*0 

21,3 

59.0 

0.22 

8,1 

0.344 

126.0 

0,580 

21.2 

0.683 

292,0 

1.75 

64.0 

0,525 

Remarks 


3G 


Wire 

Avg. 

Test  Gauge  Rum.  Net 

Rotary 

Shear 

App. 

Bob 

Circ.  Oil 

No.  &  No.  in  Defln. 

Sneed 

Rate 

Vise, 

Temp. 

Temp.  Temp. 

2A 

rns 

1/ sec. 

lb.  /ft.  sec. 

Op 

op  op 

Con^d 

16  22-1  4"  362.5 

2.48 

90.8 

O.46O 

176.8 

1.06 

38.8 

0.524 

26.5 

0.0907 

3.3 

0.92 

17.5 

14.0  21.3 

Remarks 


TABLE  3  -  G 


FOSTER  TON  CRUDE 


17 

26-2 

4" 

87.0 

1.77 

82.5 

1.65 

18 

26-2 

4" 

56.7 

0.703 

41.0 

0.468 

104.8 

1.17 

150.5 

1.66 

19 

26-2 

4" 

86.2 

1.06 

49*2 

0.573 

121.2 

1.42 

153.7 

1.78 

20 

30-3 

3" 

313.7 

1.06 

32.0 

0 . 1044 

21 

26-2 

4" 

132.0 

1,07 

28.5 

0.221 

13.5 

0,1025 

75.0 

0,583 

136.0 

1.07 

277.0 

2.15 

327.5 

2.51 

78.0 

0.58 

141.0 

1,064 

22 

26-2 

3" 

163.5 

1.42 

128.0 

1.07 

14.0 

0.105 

139.0 

1.06 

23 

24-1 

4" 

14.5 

0.104 

78.3 

0.584 

268.0 

2.13 

29*3 

0.218 

225.7 

1.75 

78.7 

0,582 

14.5 

0.1009 

64.9 

0.0252 

79.8 

60,2 

0.0258 

25.7 

0.0414 

61,8 

17.1 

0.0451 

42,8 

0.0459 

60.7 

0.0466 

60.8 

39.0 

0,0416 

61.0 

57.6 

61.8 

21.0 

0.0441 

59.0 

57.6 

59.9 

52.0 

0,0439 

65.0 

0.0443 

59.0 

59.9 

38.8 

0,0307 

69*2 

69.7 

3.8 

0.0319 

39.2 

O.6636 

49.0 

47,7 

50.8 

8.1 

0,0664 

3*7 

0.0679 

21.3 

0.0662 

39.5 

0.0655 

78.7 

O.O664 

91.9 

0.0672 

21.2 

0.0693 

39.0 

0.068 

49.0 

47.7 

50.8 

52.0 

0.0783 

44.0 

42.0 

45.7 

39.2 

0.0812 

3.8 

0.0905 

41.6 

41.4 

43.3 

38.8 

0.0888 

3.8 

0.162 

29.2 

28.4 

32.2 

21.4 

0.155 

78.0 

0.146 

31.8 

8.0 

0.156 

64.0 

0.149 

21.3 

0.157 

3.7 

0.167 

28.8 

28,4 

31.8 

1 


4G 


Wire 


Test  Gauge  Inrnu  Net 

No.  &  No.  in  Defln. 

Rotary- 

Speed 

Shear 

Rate 

2A 

r£s 

1/sec. 

24 

24-1 

4"  343.0 

1.76 

64.4 

214.2 

1.07 

39.2 

49.5 

0.222 

8.1 

23.0 

0.098 

3.6 

50.5 

0.21 

7.7 

219.5 

1.07 

3^.2 

349.0 

1.76 

64.5 

25 

24-1 

3"  27.5 

0.1045 

3.8 

249.0 

1.06 

38.8 

395.0 

1.75 

64.I 

246.5 

1.05 

38.4 

27.0 

0.1015 

3.7 

Avgo 

App.  Bob  Circ.  Oil 

Vise.  Temp.  Temp.  Temp.  Remarks 

lb.zit.sec.  °F  °F  °F 

0.225  21.0  24.4  23,9 

0.232 
0.253 
0.272 
0.266 
0.238 
0.224 

0.402  14.8  14.5  13.5 

0.358 
0.344 
0.359 
0.406 


TABLE  4  -  G 


PEMBINA  CRUDE 


26  30-3 


27  30-3 


28  30-3 


29  30-3 


72,2 

0.945 

34.6 

6.05 

8.5 

0,1027 

3.8 

6.55 

17.7 

0.22 

8.1 

6.36 

45.3 

0.578 

21.2 

6.2 

74.0 

0.947 

34.7 

6.18 

75.0 

0.947 

34.7 

6.25 

46.5 

0.578 

21.2 

6.36 

11.5 

0.10 

3.7 

9.1 

24.8 

0.219 

8,0 

8.94 

65.9 

0,580 

21.2 

9.0 

190.0 

1.67 

61.1 

9.0 

169.0 

0.932 

34.1 

14.3 

111,0 

O.566 

20.7 

15.5 

50.0 

0.214 

7.8 

18.5 

31.0 

0.0987 

3.6 

24.8 

62.0 

0.214 

7.8 

23.0 

124.0 

0.563 

20,8 

17.4 

215.0 

1.52 

55.6 

11,2 

34.5 

0.1018 

3.7 

26.8 

67.5 

0.216 

7.9 

24.7 

140.0 

0.575 

21.0 

19.2 

218.0 

1.06 

38.8 

16.3 

269.0 

1.41 

51.5 

15.1 

256,5 

1.41 

51.5 

14.4 

287.0 

1.64 

60.0 

13.8 

102.0 

0.575 

21.0 

14.1 

43.5 

0fl2l6 

7.9 

16.0 

24.5 

0.10 

3.7 

19.4 

78.8  30.6  73.3 

Loss  of  light 

80.6  78.8  ends 

69.0  69.6  69.6 

Loss  of  light 

69.6  69.6  ends 

58.5  59.6  59.4 

58.5  59.6  59.4 

59.5  60.5 

60.5 

60.0  61.0  After  7 k  mins 

61.0 


....  *. 


o 


A? 


A. 


Test 

Wire 

Gauge 

3nm 

..  Net 

Rotary 

Shear 

App. 

Bob 

Circ. 

Avg. 

Oil 

&  Ho. 

in 

Defln. 

Sneed 

Rate 

Vise. 

Tfimn. 

2A 

1/sec. 

Hwit.sec. 

of 

,L  11  , 

of 

30 

26-2 

4" 

118.0 

1.66 

60.7 

36.6 

49.0 

47.7 

50.7 

102.0 

1.42 

52.0 

37.0 

81.0 

1.066 

39.0 

39.2 

51.0 

0.585 

21.4 

44.8 

26.0 

0.224 

8,2 

59.6 

28.0 

0.105 

3.8 

83.2 

48.5 

47.7 

50.2 

31 

26-2 

4" 

252,0 

2.35 

86.0 

55.2 

40.5 

39.0 

42. 8 

138.0 

1.068 

39.1 

66.5 

92.0 

0.584 

21.4 

31.2 

48 

0.22 

8.1 

0.112 

32 

0.101 

3.7 

0.163 

39.0 

42.8 

32 

26-2 

4" 

231.5 

1.77 

64.8 

67.3 

42,0 

39.0 

44.2 

153.5 

1.07 

39.2 

73.7 

105.5 

0,586 

21.4 

92.8 

61.5 

0.226 

8.3 

141.0 

44.0 

0.107 

3.9 

212,0 

232.5 

1,77 

64.8 

67.7 

211.5 

1.78 

65.1 

61,3 

19+5.5 

1.06 

38.8 

70.9 

97,5 

0.583 

21.4 

36,2 

56.0 

0,224 

8.2 

129.0 

38.5 

0.1037 

3.3 

192.0 

42.0 

39.0 

44.2 

33 

22-1 

4" 

276.0 

2,33 

37.0 

370.0 

31,5 

30.2 

34.4 

212.0 

1.77 

64.8 

377.0 

77.0 

0.579 

21.2 

399.0 

30.5 

0.218 

8,0 

440. 0 

14,0 

0. 0982 

3,6 

443 

34.4 

Equal  intervals  of 

Thixotropic  run  -  time  to  take  equal  number  of 

readings  both  in  an  increasing 
speed  direction  and  then  in  a 
decreasing  speed  direction,  as 
follows: 


5G 


Remarks 


Breakdown  after 
8  mins. 


34  22-1  4”  14.0 

31.5 

73.5 

136.5 

267.5 

285.5 

215.0 

134.5 
77.0 
30,0 

13.5 


0.0932 

3.6 

448.0 

0,223 

8.2 

444,0 

0.532 

21,3 

424. 0 

1.06 

33.8 

405.0 

1.77 

64.8 

388.0 

2,38 

87.0 

377.0 

1.77 

64.8 

383.0 

1.05 

33.4 

403.0 

0,579 

21.2 

418.0 

0.218 

8.0 

433.0 

0,0982 

3.6 

433.0 

31.5  30.2  34.4 


31.5  30.0  34.4 


o 


6g 


Wire 

Test  Gauge  Imm, 
No.  &  No.  in 


35  22-1  31 


36  20-2  4" 


37  22-1  4" 


38  22-1  3' 


39  32-3  4" 


Net 

Defln. 

Rotaiy 

Sneed 

Shear 

Rate 

App. 

Vise. 

Bob  Circ. 
Temp.  Temp. 

Avg. 

Oil 

Temp. 

r£s 

1/sec. 

Ib^fLsec. 

Op 

op 

Op 

234.5 

2.43 

89.0 

401.0 

28.5 

20.3 

32.2 

178.0 

1.77 

64.8 

417.0 

112.5 

1.06 

38.8 

439.0 

64.7 

0.581 

21,2 

461.0 

25.5 

0.216 

7.9 

490.0 

10.8 

0.0893 

3,3 

502.0 

244.5 

2.41 

88,2 

421.0 

28.5 

20.0 

32.2 

9.0 

0.1005 

3.7 

782.0 

21.5 

14,7 

25.2 

22.0 

0.21 

7.7 

810.0 

52.0 

0.585 

21.4 

777,0 

136.7 

1.75 

64.0 

682.0 

178.5 

2.40 

87.8 

650.0 

167.7 

2.38 

87.0 

615.0 

128.5 

1.75 

64.1 

643.0 

20.0 

0.22 

8.1 

795.0 

9.5 

0.0971 

3.6 

855.0 

14.7 

25.2 

415.0 

1.4 

51.2 

935.0 

14.0 

11.7 

16.2 

328.0 

1.05 

38.4 

987.0 

198.0 

0.572 

20.9 

1090.0 

88.0 

0.212 

7.8 

1310,0 

41.0 

0.0893 

3.3 

1450.0 

14.0 

11.7 

16.2 

358,5 

In  73 

63.3 

859.0 

14,0 

17.7 

294.5 

1.41 

51.6 

868.0 

230.5 

1.06 

38.8 

916.0 

136.5 

0.568 

20.8 

997.0 

57.5 

0.208 

7.6 

1150.0 

26.5 

0.0852 

3.1 

1290.0 

17,7 

TABLE  5  - 

G 

R 

E  D  W  A  T  E  R 

CRH 

D  E 

68.2 

0.469 

17.2 

4.61 

76.0 

80.4 

76.0 

176.7 

1.18 

43.1 

4.72 

15.0 

0.1022 

3.7 

4.72 

76.0 

26 5.2 

1,64 

60.0 

5.09 

80.4 

Remarks 


4s  Min So  later 
8  Mins,  later 


Loss  of  light 
ends 


40  32-3  4"  90.3 

15.9 
35.2 
96.5 
343.2 


0.56  20.5 
0.0926  3.4 
0.208  7.6 
0.565  20.7 
1.89  69.2 


5.08  67.6  70.4  68.0 

5.41 

5.33 

5.38 

5.75 


7G 


wire  Avg. 

Test  Gauge  Iran,  Net  Rotary  Shear  App,  Bob  Giro,  Oil 

No*  4  Ho*  in  Define  Speed  Rate  Vise,  Team.  Ternr.  Temp. 

2A  £B§,  1/ sec.  Ib./frt.sec,  2L  2Z  £F 


41 

32-3 

4" 

212,3 

1.17 

38,5 

0.21 

17  »  8 

0.097 

199.8 

1.046 

320.3 

1.63 

42 

32-3 

4" 

300.0 

1.29 

25.7 

0.107 

146*0 

0,592 

366,5 

1.42 

43 

30-3 

4!' 

65.0 

0.591 

11,8 

1.04 

26,3 

0.224 

226.0 

1.773 

29,8 

0.224 

14.3 

0,105 

44 

30-3 

4" 

342,5 

1.76 

117*3 

0,586 

21,5 

0.1018 

378,5 

1.77 

128,0 

0.586 

23,0 

0,1018 

45 

30-3 

4" 

85,0 

0.22 

433*0 

1.17 

90,5 

0.222 

42.0 

0,104 

46 

30-3 

qn 

325.0 

1.41 

25.0 

0.15 

347.0 

1.41 

47 

30-3 

3" 

316.5 

1.42 

135,0 

0,586 

54.5 

0.222 

48 

26-2 

4" 

81.0 

1.06 

19.5 

0.224 

10,0 

0,102 

22.2 

0.223 

166.5 

2.13 

42,8 

5*57 

56.8 

59*0 

58.2 

7.7 

5.64 

3.5 

5*65 

38.3 

5.39 

59.6 

6.06 

56.8 

59*0 

58.2 

47*2 

7.16 

51.0 

52.6 

3.9 

7.44 

49*0 

50.6 

21.7 

7.61 

52,0 

7.98 

49*0 

50.6 

21.6 

8.7 

41.6 

40.1 

43.9 

38.0 

8.95 

8.2 

9,3 

65*0 

10,1 

40.0 

40,0 

42.3 

8,2 

10.5 

3*8 

10.8 

40.0 

42,3 

64.5 

15.4 

33*5 

27*7 

36.5 

21,4 

15.9 

3*7 

16.7 

33*5 

27,7 

36.5 

64.8 

16.9 

21.4 

17.3 

3.7 

17*9 

8.1 

30.2 

21.5 

17.6 

25.1 

42,3 

29*2 

8.1 

32.3 

3.8 

32.0 

51*6 

24.0 

23.0 

20.3 

26.5 

5*5 

25.0 

26.5 

51.6 

25*6 

23.0 

20,3 

52.0 

23.2 

23*5 

18,3 

26.8 

21,4 

24.0 

8,1 

25*6 

24.0 

18.3 

27.3 

33.3 

39*2 

14*5 

9*3 

18.3 

3.2 

45.0 

3.7 

52,4 

13.5 

9.3 

17*3 

8,2 

51*3 

16.8 

78,0 

40.2 

13.0 

9*3 

Remarks 


Loss  of  light 
ends 


